350 CONTINUE AIRY274 



T=10./C-l. AIRY275 



TT=T+T AIRY276 



J=N4 AIRY277 



F1=A(J) AIRY278 



E1=B(J) AIRY279 



F2=0. AIRY280 



E2=0. AIRY281 



DO 310 1=1 ,M4 AIRY282 



J=J-1 AIRY283 



TEMP1=F1 AIRY284 



TEMP2=E1 AIRY285 



F1=TT*F1-F2+A<J) AIRY286 



E1=TT*E1-E2+B(J) AIRY287 



F2=TEMP1 AIRY28B 



E2=TEMP2 AIRY289 



310 CONTINUE AIRY290 



TEMP1=T*F1-F2+A(1) AIRY291 



TEMP2=T#Ei-E2+B(l) AIRY292 



RTRX=SQRT(RX) AIRY293 



CV=C-FPI12 AIRY294 



CCV=COS(CV) AIRY295 



SCV=SIN(CV) AIRY296 



AI=(TEMP1«CCV-TEMP2*SCV)/RTRX AIRY297 



J=N4D AIRY298 



F1=DA(J) AIRY299 



E1=DB(J) AIRY300 



F2=0. AIRY301 



E2=0. AIRY302 



DO 311 1=1, M4D AIRY303 



J=J-1 AIRY304 



TEMP1=F1 AIRY305 



TEMP2=E1 AIRY306 



F1=TT*F1-F2+DA(J) AIRY307 



E1=TT*E1-E2+DB(J) AIRY308 



F2=TEMP1 AIRY309 



E2=TEMP2 AIRY310 



311 CONTINUE AIRY311 



TEMP1=T#F1-F2+DA(1) AIRY312 



TEMP2=T*E1-E2+DB(1) AIRY313 



E1=CCV*C0N5+.5*SCV AIRY314 



E2=SCV*C0N5-.5*CCV AIRY315 



DAI=(TEMPI*E1-TEMP2*E2)*RTRX AIRY316 



RETURN AIRY317 



END AIRY318 



FUNCTION GAMLN(X) 6LN1 



c GLN2 



C WRITTEN BY D. E. AMOS, SEPTEMBER, 1977. GLN3 



C REFERENCES GLN4 



C * SAND-77-1518 GLN5 

 C * COMPUTER APPROXIMATIONS BY J.F.HART, ET.AL., SIAM SERIES IN GLN6 



C APPLIED MATHEMATICS, WILEY, 1968, P. 135-136. GLN7 



C * NBS HANDBOOK OF MATHEMATICAL FUNCTIONS, AMS 55, BY GLN3 



C M. ABRAMOWITZ AND I. A. STEGUN, DECEMBER. 1955, P. 257. GLN9 



C ABSTRACT GLN10 



Figure A3. (Sheet 22 of 25) 

 A26 



