#!/usr/bin/perl
#
#

$pi = 3.14159;
$a = 1.0;

$start_r = 0.0;
$end_r = 100.0;
$dr = 0.05;

for ($r = $start_r; $r < $end_r; $r += $dr) {

  $radfunc = ( sqrt(1.0/$pi)*($a**(-1.5)) * exp( -$r/$a ) );
  $dvol = $r*$r*$dr;

  $radfunc_squared = $radfunc*$radfunc;
  $prob = $radfunc_squared*$dvol;

  printf "%8.3f %10.5e %10.5e %10.5e %10.5e \n",
     $r, $radfunc, $dvol, $radfunc_squared, $prob;

}


exit 0;
