C NZ - ERROR INDICATOR BESJ135 



C NZ=0 NORMAL RETURN - COMPUTATION COMPLETED BESJ136 



C NZ=-1 X IS LESS THAN 0.0 BESJ137 



C NZ=-2 ALPHA IS LESS THAN 0.0 BESJ13B 



C NZ=-3 N IS LESS THAN 1 BESJ139 



C NZ.GT.O LAST NZ COMPONENTS OF Y SET TO 0.0 BESJ140 



C BECAUSE OF UNDERFLOW BESJ141 



C BESJ142 



C ERROR CONDITIONS BESJ143 



C BESJ144 



C IMPROPER INPUT ARGUMENTS - A FATAL ERROR BESJ145 



C UNDERFLOW - A NON-FATAL ERROR (NZ.BT.O) BESJ146 



C BESJ147 



c -BESJ148 



C BESJ149 



DOUBLE PRECISION DX , TRX , DTM.DFN BESJ150 



DIMENSION Y(N) BESJ151 



DIMENSION C(ll, 10) , ALFA (26,4) , BETA (26,5) BESJ152 



DIMENSION CI (88) ,C2 (22) BESJ153 



DIMENSION AK52) ,A2(52) ,B1(52) ,B2<52) ,B3(26) BESJ154 



DIMENSION GAMA(26) , TEMP (3) ,KMAX (5) ,AR(B) ,BR(10) ,UP0L(10) BESJ 155 



DIMENSION FNULIM(2) ,PP(4) BESJ156 



DIMENSION CR(10) ,DR(10) BESJ157 



C BESJ15B 



EQUIVALENCE (C (1 , 1 ) ,C1 ( 1 ) ) BESJ159 



EQUIVALENCE (C ( 1 ,9) ,C2 ( 1 ) ) BESJ160 



EQUIVALENCE ( ALFA ( 1 , 1 ) , Al ( 1 ) ) BESJ161 



EQUIVALENCE ( ALFA ( 1 , 3) , A2 ( 1 ) ) BESJ162 



EQUIVALENCE (BETA ( 1 , 1 ) ,B1 ( 1) ) BESJ163 



EQUIVALENCE (BETA ( 1 , 3) , B2 ( li i BESJ 164 



EQUIVALENCE ( BETA ( 1 , 5) , B3 ( 1 ) ) BESJ165 



C BESJ166 



DATA ELIM1,ELIM2,T0L / 667. , 644. , l.E-15 /BESJ167 



C BESJ168 



DATA PP(1) /8.7290915393555E+00/, PP (2) 12. 6569393226503E-01 / , BESJ169 



1 PP(3)/1.2457357686559E-01/, PP (4) /7. 7013374743039E-04/ BESJ170 



C BESJ 171 



C T0LS=LN(l.E-3) BESJ172 



DATA TOLS / -6 . 907755278982 1 E + 00 / BESJ173 



C BESJ174 

 DATA UPOL(l) ,C0Nl,CON2,CON3,CON548 /' 1 . OOOOOOOOOOOOOE + 00, BESJ 175 



1 6.6666666666667E-01, 3. 3333333333333E-01 , 1 . 4 1 42 1 3562373 1 E + 00 , BESJ176 



2 1.0416666666667E-01/ BESJ177 

 C BESJ178 



DATA RTWO,PDF,RTTP,PIDT / 1 . 3483997249265E + 00, BESJ179 



1 7.8539816339745E-01, 7. 9788456080286E-01 . 1 . 5707963267949E+00/ BESJ180 



C BESJ181 



DATA FNULIM ( 1 ) /' 100. / , FNUL IM (2) /60. / BESJ182 



C BESJ183 



C CE=-ALOG(TOL) . TCE=-0. 75#AL0G (TOD BESJ184 



DATA CE , TCE / 3. 453877639491 1E+01 , 2. 5904082296 183E + 01 / BESJ185 



C BESJ186 



DATA INLIM / 150 / BESJ 187 



C BESJ188 



DATA AR ( 1 ) /8.3550347222222E-02/, AR (2) / 1 . 2822657455633E-01 / , BESJ139A 



Figure A3. (Sheet 3 of 25) 

 A7 



