C ZETA AND TRUNCATION FOR A(ZETA) AND B(ZETA) SERIES BESJ295 



C BESJ296 



C KMAX IS TRUNCATION INDEX FOR A(ZETA) AND B(ZETA) SER I ES=MAX (2 , SA) BESJ297 



C BESJ298 



SA=0. BESJ299 



IF(ABW2.EQ.O. ) GO TO 21 BESJ300 



SA=T0LS/AL0G(ABW2) BESJ301 



21 SB=SA BESJ302 

 DO 22 1=1,5 BESJ303 

 KMAX(I)=AMAX1 (SA,2. ) BESJ304 

 SA=SA+SB BESJ305 



22 CONTINUE BESJ306 

 KB=KMAX(5) BESJ307 

 KLAST=KB-1 BESJ30B 

 SA=GAMA(KB) BESJ309 

 DO 24 K=1,KLAST BESJ310 

 KB=KB-1 BESJ311 

 SA=SA*W2+GAMA(KB) BESJ312 



24 CONTINUE BESJ313 



Z=W2#SA BESJ314 



AZ=ABS(Z) BESJ315 



RTZ=SQRT(AZ) BESJ316 



FN13=FN**C0N2 BESJ317 



RTARY=RTZ*FN13 BESJ31B 



ARY=-RTARY*RTARY BESJ319 



AZ32=AZ*RTZ#C0N1 BESJ320 



ACZ=FN*AZ32 BESJ321 



IFfZ.LE.O. ) GO TO 27 BESJ322 



C BESJ323 



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



C UP FROM UNDERFLOW LIMIT OF CDC 6600 BESJ325 



C BESJ326 



IFIACZ.GT.ELIM2) GO TO 180 BESJ327 



ARY=-ARY BESJ32B 



27 PHI = SQRT(SQRT(SA + SA + SA + SA) ) BESJ329 



C BESJ330 



C B(ZETA) FOR S=0 BESJ331 



C BESJ332 



KB=KMAX(5) BESJ333 



KLAST=KB-1 BESJ334 



SB=BETA(KB,1) BESJ335 



DO 23 K=1,KLAST BESJ336 



KB=KB-1 BESJ337 



SB=SB*W2+BETA(KB,1) BESJ338 



23 CONTINUE BESJ339 

 K S P 1 = 1 BESJ340 

 FN2=FN*FN BESJ341 

 RFN2=1./FN2 BESJ342 

 RDEN=1. BESJ343 

 ASUM=1. BESJ344 

 RELB=TOL*ABS(SB) BESJ345 

 BSUM=SB BESJ346 

 DO 25 KS=1 ,4 BESJ347 

 KSP1=KSP1+1 BESJ348 

 RDEN=RDEN*RFN2 BESJ349 



Figure A3. (Sheet 9 of 25) 

 A13 



