KT=1 BESJ460 



401 IS=2 BESJ461 



DFN=DFN-l.D+0 BESJ462 



FN=DFN BESJ463 



GO TO 130 BESJ464 



C BESJ465 



C SERIES FOR (X/2) **2. LE. NU+1 BESJ466 



C BESJ467 



134 CONTINUE BESJ468 



GLN=GAMLN(FNP1) BESJ469 



ARG=FN#X02L-GLN BESJ470 



IF (ARB. LT. -ELIM1 ) GO TO 123 BESJ471 



EARG=EXP(ARG) BESJ472 



300 CONTINUE BESJ473 

 S=l. BESJ474 

 AK=3. BESJ475 

 T2=l. BESJ476 

 T=l. BESJ477 

 S1=FN BESJ478 

 DO 125 K = 1 , 1 7 BESJ479 

 S2=T2+S1 BESJ480 

 T=-T*SX02/S2 BESJ481 

 S = S + T BESJ482 

 IF(ABS(T) .LT.TOL) GO TO 127 BESJ483 

 T2=T2+AK BESJ484 

 AK=AK+2. BESJ485 

 S1=S1+FN BESJ486 



125 CONTINUE BESJ487 



127 CONTINUE BESJ488 



TEMP(IS)=S*EARG BESJ489 



GO TO (301,202,600). IS BESJ490 



301 EARG=EARG«FN/X02 BESJ491 

 DFN=DFN-l.D+0 BESJ492 

 FN=DFN BESJ493 

 IS=2 BESJ494 

 GO TO 300 BESJ495 



C BESJ496 



C SET UNDERFLOW VALUE AND UPDATE PARAMETERS BESJ497 



C BESJ498 



180 Y(NN)=0. BESJ499 



NN=NN-1 BESJ500 



DFN=DFN-l.D+0 BESJ501 



FN=DFN BESJ502 



IF (NN-1) 170,171,130 BESJ 50: 



171 KT=2 BESJ504 

 IS=2 BESJ505 

 GO TO 130 BESJ506 



123 Y(NN)=0. BESJ507 



NN=NN-1 BESJ508 



FNP1=FN BESJ509 



DFN=DFN-l.D+0 BESJ510 



FN=DFN BESJ511 



IF (NN-1 ) 170,172,173 BESJ512 



172 KT=2 BESJ513 

 IS=2 BESJ514 



Figure A3. (Sheet 12 of 25) 

 A16 



