

200 CONTINUE 





TAU=1./RA 





T2=l./W2 





IF(W2.GE.O. ) BO TO 30 



c 





c 



CASES FOR (X/FN) .ST.SQRTd 



c 





C BESJ350 



C A(ZETA) AND B(ZETA) FOR 3=1,2,3,4 BESJ351 



C BESJ352 



KB=KMAX(5-KS) BESJ353 



KLAST=KB-1 BESJ354 



SA=ALFA(KB,KS) BESJ355 



SB=BETA(KB,KSP1) BESJ356 



DO 26 K=1,KLAST BESJ357 



KB=KB-1 BESJ358 



SA=SA*W2+ALFA(KB,KS> BESJ359 



SB=SB*W2+BETA(KB,KSP1) BESJ360 



26 CONTINUE BESJ361 



TA=SA*RDEN BESJ362 



TB=SB*RDEN BESJ363 



ASUM=ASUM+TA BESJ364 



BSUM=BSUM+TB BESJ365 



IF(ABS(TA) .LE.TOL.AND.ABS(TB).LE.RELB) GO TO 152 BESJ366 



25 CONTINUE BESJ367 



152 CONTINUE BESJ368 



BSUM=BSUM/(FN*FN13) BESJ369 



GO TO 400 BESJ370 



C BESJ371 



BESJ372 

 BESJ373 

 BESJ374 

 BESJ375 

 BESJ376 



!775) BESJ377 



BESJ378 



AZ32=ABS(RA-ATAN(RA) ) BESJ379 



ACZ=AZ32*FN BESJ380 



CZ=-ACZ BESJ381 



Z32=1.5*AZ32 BESJ382 



RTZ=Z32**C0N2 BESJ383 



FN13=FN**C0N2 BESJ384 



RTARY=RTZ*FN13 BESJ385 



ARY=-RTARY*RTARY BESJ386 



GO TO 150 BESJ387 



30 CONTINUE BESJ38B 



C BESJ389 



C CASES FOR (X/FN) .LT.SQRT(0. 7225) BESJ390 



C BESJ391 



AZ32=ABS(AL0G( (1.+RA1/XX) -RA) BESJ392 



C BESJ393 



C TEST FOR UNDERFLOW, 1 . E-280 = EXP(-644.), ONE WORD LENGTH BESJ394 



C UP FROM UNDERFLOW LIMIT OF CDC 6600 BESJ395 



C BESJ396 



ACZ=AZ32*FN BESJ397 



CZ=ACZ BESJ398 



IF(ACZ.BT.ELIM2) GO TO 180 BESJ399 



Z32=1.5*AZ32 BESJ400 



RTZ=Z32**C0N2 BESJ401 



FN13=FN#*C0N2 BESJ402 



RTARY=RTZ*FN13 BESJ403 



ARY=RTARY*RTARY BESJ404 



Figure A3. (Sheet 10 of 25) 

 A14 



