173 IF <SX02. LE. FNP1 ) GO TO 133 BESJ515 



BO TO 130 BESJ516 



133 ARG=ARG-X02L+AL0G(FNP1> BESJ517 



IF(ARG.LT.-ELIMl) GO TO 123 BESJ51B 



GO TO 134 BESJ519 



170 NZ=N-NN BESJ520 



RETURN BESJ521 



C BESJ522 



C BACKWARD RECURSION SECTION BESJ523 



C BESJ524 



202 CONTINUE BESJ525 

 NZ=N-NN BESJ526 

 IF(KT.EQ.2) GO TO 250 BESJ527 



203 CONTINUE BESJ528 

 C BACKWARD RECUR FROM INDEX ALPHA+NN-1 TO ALPHA BESJ529 



Y(NN>=TEMP(1) BESJ530 



Y(NN-1)=TEMP(2> BESJ531 



IF(NN.EQ.2) RETURN BESJ532 



DX=X BESJ533 



TRX=2.D+0/DX BESJ534 



DTM=DFN*TRX BESJ535 



TM=DTM BESJ536 



K=NN+1 BESJ537 



DO 230 1=3, NN BESJ53B 



K=K-1 BESJ539 



Y(K-2)=TM*Y(K-1)-Y(K) BESJ540 



DTM=DTM-TRX BESJ541 



TM=DTM BESJ542 



230 CONTINUE BESJ543 



RETURN BESJ544 



250 Y(1)=TEMP(2) BESJ545 



RETURN BESJ546 



C BESJ547 



C ASYMPTOTIC EXPANSION FOR X TO INFINITY WITH FORWARD RECURSION IN BESJ548 



C OSCILLATORY REGION X.GT.MAX(20, NU), PROVIDED THE LAST MEMBER BESJ549 



C OF THE SEQUENCE IS ALSO IN THE REGION. BESJ550 



C BESJ551 



500 CONTINUE BESJ552 



IN=ALPHA-TAU+2. BESJ553 



IF(IN.LE.O) GO TO 502 BESJ554 



I NP 1 = I N+ 1 BESJ555 



DALPHA=ALPHA-FL0AT(INP1) BESJ556 



KT=1 BESJ557 



GO TO 511 BESJ558 



502 DALPHA=ALPHA BESJ559 

 IN=0 BESJ560 



511 IS=KT BESJ561 



512 ARG=X-PIDT*DALPHA-PDF BESJ562 

 SA=SIN(ARG) BESJ563 

 SB=COS(ARG) BESJ564 

 RA=RTTP/RTX BESJ565 

 ETX=8.*X BESJ566 



503 DX=DALPHA BESJ567 

 DX=DX+DX BESJ568 

 DTM=DX*DX BESJ569 



Figure A3. (Sheet 13 of 25) 

 A17 



