#3_3.pl
#loop over drop radius from 1 micron to 100 microns in steps of 2 microns
#use each radius to compute cross section
#if radius <20 microns, use (pi)*(r^2)*a*(1-exp(-c*r))
#if radius >= 20 microns, use 2*pi*(r^2)

#VARIABLE DECLARATION
#--------------------
my ($r, $rad, $twopi, $sig, $pi, $a0, $c0);
my (@radii, @rads, @sigs);
#--------------------

use strict;

@radii = (1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 91, 93, 95, 97, 99);
@rads = ();
@sigs = ();

$pi = 3.14159;
$twopi = (2)*(3.14159);
$a0 = 1.18;
$c0 = (0.28e6);

foreach $r (@radii) 
	{
	$rad = (1e-6)*$r;
	push (@rads, $rad);
	if ($r>20)
		{
		$sig = $twopi*($rad**2);
		push (@sigs, $sig);
		}
	else
		{
		$sig = $pi*($rad**2)*$a0*(1-exp((-$c0)*$rad));
		push (@sigs, $sig);
		}
	}
	
print "@rads \n";
print "@sigs \n";