my ($r, $rad, $twopi, $sig, $pi, $a0, $c0);

##</code><code>##

    if ($r>20)
        {
        my $sig = $twopi*($rad**2);
        push (@sigs, $sig);
        }
    else
        {
        my $sig = $pi*($rad**2)*$a0*(1-exp((-$c0)*$rad));
        push (@sigs, $sig);
        }