$sig = ($r>20)
         ? $twopi*($rad**2)
         : $pi*($rad**2)*$a0*(1-exp((-$c0)*$rad));
    push (@rads, $rad);
    push (@sigs, $sig);