my $PI = 2 * atan2(1, 0);

##</code><code>##

my $PI = 2 * atan2(1, 0);
my $A = 1.18;
my $C = 0.28e6;
my $HWM = 2e-5;
my $SCALE = 1e-6;

sub sigma;

my (@rads, @sigs);
for (my $i = 1; $i <= 100; $i += 2)) {
    my $r = $i * 1e-6;
    push @rads, $r;
    push @sigs, sigma($r, $r >= $HWM)
}
print "radii: @rads\n";
print "sigmas: @sigs\n";
    
sub sigma {
    my ($rad, $high) = @_;
    $PI * $rad ** 2 * ($high ? 2 : $A * (1 - exp(-$C * $rad)));
}