

S2=TM*S2-S1 



BESJ625 





TM=TM+TX 



BESJ626 





S1=S 



BESJ627 



520 



CONTINUE 



BESJ628 



C 





BESJ629 



C 



FORWARD RECUR FROM INDEX ALPHA TO ALPHA+N-1 



BESJ630 



C 





3ESJ631 





Y(1)=S1 



BESJ632 





Y(2)=S2 



BESJ633 





IF(NN.EQ.2) RETURN 



BESJ634 





DO 530 1=3, NN 



BESJ635 





Y(I)=TM*Y(I-l)-Y(I-2) 



BESJ636 





TM=TM+TX 



BESJ637 



530 



CONTINUE 



BESJ638 





RETURN 



BESJ639 



535 



Y(1)=S2 



BESJ640 





RETURN 



BESJ641 



C 





BESJ642 



c 



BACKWARD RECURSION WITH NORMALIZATION BY 



BESJ643 



c 



ASYMPTOTIC EXPANSION FOR NU TO INFINITY OR POWER SERIES. 



BESJ644 



c 





BESJ645 



600 



CONTINUE 



BESJ646 



C 



COMPUTATION OF LAST ORDER FOR SERIES NORMALIZATION 



BESJ647 





KM=AMAX1 (3.-FN,0. ) 



BESJ648 





TFN=FN+FLOAT(KM) 



BESJ649 





TA= (GLN+TFN-O. 9189385332-0. 0833333333/TFN) / (TFN+O. 5) 



BESJ650 





TA=X02L-TA 



BESJ651 





TB=-(1.-1.5/TFN)/TFN 



BESJ652 





IN=CE/(-TA+SQRT(TA*TA-CE#TB) )+1.5 



BESJ653 





IN=IN+KM 



BESJ654 





GO TO 603 



BESJ655 



650 



CONTINUE 



BESJ656 



c 



COMPUTATION OF LAST ORDER FOR ASYMPTOTIC EXPANSION NORMALIZATION 



BESJ657 





BLN=AZ32+RA 



BESJ658 





IF(ARY.GT.30. ) GO TO 675 



BESJ659 





RDEN=(PP(4)»ARY+PP(3i )*ARY+1. 



BESJ660 





RZDEN=PP(1)+PP(2)*ARY 



BESJ661 





TA=RZDEN/RDEN 



BESJ662 





I F ( W 2 . L T . . 10) GO TO 651 



BESJ663 





TB=GLN/RTARY 



BESJ664 





GO TO 677 



BESJ665 



651 



TB= (1.25992 1049 + 0. 1679894730 *W2) /FN 13 



BESJ666 





GO TO 677 



BESJ667 



675 



CGNTINUE 



BESJ668 





TA=C0N1*TCE/ACZ 



BESJ669 





TA=( (0.0493827 160*TA-0. 1 1 1 1 1 1 1 1 1 1 ) *TA + 0. 6666666667) *TA*ARY 



BESJ670 





IF(W2.LT.0. 10) GO TO 651 



BESJ671 





TB=GLN/RTARY 



BESJ672 



677 



IN=TA/TB+1.5 



BESJ673 





IFdN.GT. INLIM) GO TO 402 



BESJ674 



603 



DX=FLOAT(IN) 



BESJ675 





DTM=DFN+DX 



BESJ676 





DX = X 



BESJ677 





TRX=2.D+0/DX 



BESJ678 





DTM=DTM*TRX 



BESJ679 



Figure A3. (Sheet 15 of 25) 



A19 



