#!/usr/bin/perl
#

$G = 6.67E-11;
$Msun = 1.99E30;
$AU = 1.50E11;


# $vi = 19580.0;
$vi = 21860.0;

$ri = 5.27*$AU;
$rf = 32.83*$AU;


$ke_i = 0.5*$vi*$vi;
$gpe_i = 0.0 - $G*$Msun/$ri;
$tot_i = $ke_i + $gpe_i;

printf "  ke_i %12.5e  GPE_i  %12.5e   tot_i  %12.5e \n",
        $ke_i, $gpe_i, $tot_i;

$gpe_f = 0.0 + $G*$Msun/$rf;
$ke_f = $tot_i + $gpe_f;

printf "  ke_f %12.5e  GPE_f  %12.5e   tot_f  %12.5e \n",
               $ke_f, $gpe_f, $tot_i;



$vf = sqrt(2.0*$ke_f);

printf " vf is %12.5e \n", $vf;


exit 0;