DFN=DBLE (FLOAT (N) ) +DBLE (ALPHA) -1 . D+O BESJ240 



FNU=DFN BESJ241 



X02=X#.5 BESJ242 



SX02=X02*X02 BESJ243 



C BESJ244 



C DECISION TREE FOR REGION WHERE SERIES, ASYMPTOTIC EXPANSION FOR X BESJ245 



C TO INFINITY AND ASYMPTOTIC EXPANSION FOR NU TO INFINITY ARE BESJ246 



C APPLIED. BESJ247 



C BESJ248 



IF (SX02. LE. (FNU+1. ) ) GO TO 850 BESJ249 



TA=AMAX1(20. ,FNU) BESJ250 



IF(X.GT.TA) GO TO 880 BESJ251 



IF (X . ST. 12. ) GO TO 860 BESJ252 



X02L=AL0G(X02) BESJ253 



NS=SX02-FNU BESJ254 



GO TO 852 BESJ255 



850 FN=FNU BESJ256 



FNP1=FN+1. BESJ257 



X02L=AL0G(X02) BESJ258 



IS=KT BESJ259 



IF(X.LE.0.5) GO TO 134 BESJ260 



NS=0 BESJ261 



852 DFN=DFN+DBLE(FLOAT(NS) ) BESJ262 



FN=DFN BESJ263 



FNP1=FN+1. BESJ264 



IS=KT BESJ265 



IF(N-l+NS.GT.O) IS=3 BESJ266 



GO TO 134 BESJ267 



860 NS=AMAX1(36.-FNU,0. ) BESJ268 



DFN=DFN+DBLE(FLOAT(NS) ) BESJ269 



FN=DFN BESJ270 



IS=KT BESJ271 



IF(N-l+NS.GT.O) IS=3 BESJ272 



GO TO 130 BESJ273 



880 CONTINUE BESJ274 



RTX=SQRT(X) BESJ275 



TAU=RTWO#RTX BESJ276 



TA=TAU+FNULIM(KT) BESJ277 



IF(FNU.LE.TA) GO TO 500 BESJ278 



129 FN=FNU BESJ279 

 IS=KT BESJ280 



C BESJ281 



C UNIFORM ASYMPTOTIC EXPANSION FOR NU TO INFINITY BESJ282 



C BESJ283 



130 CONTINUE BESJ284 

 XX=X/FN BESJ285 

 W2=1.-XX#XX BESJ286 

 ABW2=ABS(H2) BESJ287 

 RA=SQRT(ABW2i BESJ28S 

 IF(ABW2.GT. 0.2775) GO TO 200 BESJ289 



C BESJ290 



C CASES NEAR X = FN, ABS ( 1 . - ( X /FN) **2 ) . LE. 0. 2775 BESJ291 



C COEFFICIENTS OF ASYMPTOTIC EXPANSION BY SERIES BESJ292 



C BESJ293 



C BESJ294 



Figure A3. (Sheet 8 of 25) 

 A12 



