F17(p,g)Ne18
$T9	NA(av)    (cm3 mole-1/sec)
0.1	2.68E-9		  0.38E-9  
0.2	5.15E-6		  0.75E-6  
0.3	1.97E-4		  0.29E-4  
0.4	2.29E-3		  0.40E-3  
0.5	1.77E-2		  0.49E-2  
0.6	9.29E-2		  3.28E-2  
0.7	3.32E-1		  1.30E-1  
0.8	8.80E-1		  3.61E-1  
0.9	1.88E0		  0.78E0   
1.0	3.43E0		  1.44E0   
1.5	1.97E1		  0.78E1   
2.0	4.62E1		  1.64E1   



Subroutine code in PERL format:
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
use strict;
my $T9=$ARGV[0];
print "#Bardayan2002: F17(p,g)Ne18 rate:\n";
unless($T9) { print "Input T9:"; chop($T9=<STDIN>);}
my $rate=rate($T9);
printf "F17(p,g)Ne18 \t %E \t0 \t #Bardayan2002\n",$rate;


sub rate 
{
my  $T9=$_[0];
my  $T913=$T9**(1/3);
my  $T953=$T9**(5/3);
my $rate=exp(2.6357E1 + 6.7829E-2/$T9-2.0195E1/$T913- 1.0068E1*$T913 + 4.7875E0*$T9
	           - 1.1482E0*$T953-4.1554E-1*log($T9))
           + exp(3.2659E1-9.1645E-1/$T9-1.5977E2/$T913
		   + 1.2833E2*$T913 + 2.6893E0*$T9 - 2.0447E0*$T953
		   - 9.0961E1*log($T9));
}
#=============================================================================================