

T! 



XX=2. +T 



DO 24 1=1, NXM 



DGAM=DGAH*XX 





24 



XX=XX+1. 





25 



GAMLN=ALOG(DGAN 

 RETURN 



c 







c 





X.GT. XLIM1 



c 







GLN121 

 GLN122 



GLN123 

 GLN124 

 GLN125 

 GLN126 

 GLN127 

 GLN128 

 GLN129 



40 RX=1. /X GLN130 

 RXX=RX*RX GLN131 

 IF( (X-XLIM2) . LT. 0. ) GO TO 41 GLN132 

 PX=Q(i)*RXX+Q(2) GLN133 

 GAMLN=PX*RX+(X-.5)*AL0G(X)-X+RTWPIL GLN134 

 RETURN GLN135 



C GLN136 



C X.LT.XLIM2 GLN137 



C GLN138 



41 PX=P(1) GLN139 

 SUM=(X-.5i *ALOG(X)-X GLN140 

 DO 42 1=2,5 GLN141 

 PX=PX*RXX+P ( I ) GLN142 



42 CONTINUE GLN143 

 GAMLN=PX*RX+SUM+RTWPIL GLN144 

 RETURN GLN145 



C GLN146 



C TABLE LOOK UP FOR INTEGER ARGUMENTS LESS THAN OR EQUAL 100. GLN147 



C GLN148 



51 IF (NDX. GT. 100) GO TO 40 GLN149 



GAMLN=GLN(NDX) GLN150 



RETURN GLN151 



END GLN152 



Figure A3. (Sheet 25 of 25) 



A29 



