5^ 



1=1- e V* V ° M^(y)M^y)(Z 2 + yJdZ 

 7 o J o Vz2+2>T 



In this form, the integrand behaves properly, (there is no longer 

 a singularity at y = y ) and the integral converges rapidly. 



The integral was actually computed by using the form in [5]. The 

 numerical integration was performed once using Simpson's rule and a second 

 time using the trapezoidal rule— for checking purposes. The interval 4Z was 

 taken as 0.1 and the range extended from 0.0 to approximately J>. 5. The M 

 functions were computed from previous tables by using 4-point Lagrangian in- 

 terpolation. The exponential function was computed from tables and the use 

 of the approximation e~ x =1 - x + x 2 /2, x < 0.01 . The algebraic function 

 in the integrand was computed in straight form and fashion. All of this work 

 was done on the IBM electronic calculator (type 60U) and the auxiliary IBM 

 punch card equipment. All the IBM operations were checked— independently 

 wherever possible. 



