150 CONTINUE BESJ405 

 PHI=SQRT( (RTZ+RTZ)*TAU> BESJ406 

 TB=1. BESJ407 

 ASUM=1. BESJ408 

 TFN=TAU/FN BESJ409 

 UP0L(2)=(C(1,1)*T2+C(2,1) )#TFN BESJ410 

 RCZ=C0N1/CZ BESJ411 

 CRZ32=C0N548*RCZ BESJ412 

 BSUM=UP0L(2)+CRZ32 BESJ413 

 RELB=TOL*ABS(BSUM) BESJ414 

 AP=TFN BESJ415 

 KS=0 BESJ416 

 KP1=2 BESJ417 

 RZDEN=RCZ BESJ418 

 DO 155 LR=2,8,2 BESJ419 



C BESJ420 



C COMPUTE TWO U POLYNOMIALS FOR NEXT A(ZETA) AND B(ZETA) BESJ421 



C BESJ422 



LRP1=LR+1 BESJ423 



DO 101 K=LR,LRP1 BESJ424 



KS=KS+1 BESJ425 



KP1=KP1+1 BESJ426 



S1=C(1,K) BESJ427 



DO 102 J=2,KP1 BESJ428 



S1=S1*T2+C(J,K) BESJ429 



102 CONTINUE BESJ430 



AP=AP*TFN BESJ431 



UP0L(KP1)=AP*S1 BESJ432 



CR(KS)=BR(KS)*RZDEN BESJ433 



RZDEN=RZDEN#RCZ BESJ434 



DR(KS)=AR(KS)*RZDEN BESJ435 



101 CONTINUE BESJ436 



SUMA=UP0L(LRP1) BESJ437 



SUMB=UP0L(LR+2)+UP0L(LRPl)*CRZ32 BESJ438 



JU=LRP1 BESJ439 



DO 151 JR=1,LR BESJ440 



JU=JU-1 BESJ441 



SUMA=SUMA+CR(JR)*UPOL(JU) BESJ442 



SUMB=SUMB+DR(JR)*UPOL(JU) BESJ443 



151 CONTINUE BESJ444 

 TB=-TB BESJ445 

 IF(W2.GT.O.) TB=ABS(TB) BESJ446 

 ASUM=ASUM+SUMA*TB BESJ447 

 BSUM=BSUM+SUMB*TB BESJ448 

 IF(ABS(SUMA).LE.TOL.AND.ABS(SUMB).LE.RELB) GO TO 165 BESJ449 



155 CONTINUE BESJ450 



165 TB=RTARY BESJ451 



IF(W2.GT.0.) TB=-TB BESJ452 



BSUM=BSUM/TB BESJ453 



C BESJ454 



400 CONTINUE BESJ455 



CALL JAIRY(ARY,RTARY,ACZ,AI ,DAI) BESJ456 



TEMP(IS)=PHI*(AI*ASUM+DAI*BSUM)/FN13 BESJ457 



GO TO (401,202,650), IS BESJ458 



402 TEMPU)=TEMP<3) BESJ459 



Figure A3. (Sheet 11 of 25) 

 A15 



