171 IF (RZR .GT. 4.4) GO TO 120 



172 TEM1 = 7. - .2632 * RZR**2 



173 IF (TEMP .GT. TEMI ) GO TO 120 



174 GO TO 53 



175 51 IF (RZR .LT. -9.) GO TO 120 



176 TEMI = 4.4 + .1375 * RZR 



177 IF (TEMP .GT. TEMI) GO TO 120 



178 53 FLPS = 1 



179 ST0RE1 = ZR*ZR-ZI*ZI 



180 ST0RE2 = 2.*ZR*ZI 



181 XR = ST0RE1*ZR -ST0RE2*ZI 



182 XI = ST0RE1*ZI +ST0RE2*ZR 



183 DO 55 MLS=1 , LA2 



184 IF (ZMSO - ZMLA2(MLS)) 62,62,55 



185 55 CONTINUE 



186 62 FR = AG 



187 FI = 0.0 



188 XPR(l) = XR 



189 XPI(1 ) = XI 



190 DO 65 M = 1 ,MLS 



191 FR=FR + A(P/!)*XPR(M) 



192 FI=FI+A(M)*XPI (M) 



193 XPR(M+1 )=XR*XPR(M)-XI*XPI(M) 



194 XPI(M+1 )=XI*XPR(M)+XR*XPI (M) 



195 65 CONTINUE 



196 GR=BO 



197 GI=0.0 



198 DO 72 M = 1 ,MLS 



199 GR=GR+B(M)+XPR(M) 



200 GI=GH-B(M)+XPI(IVI) 



201 72 CONTINUE 



202 X =ZRtGR-ZI*GI 



203 GI=ZR*GI+ZI*GR 



204 GR=X 



205 SR=-C1+ (GI-FI-FI ) 



206 SI=C1 ♦(GR-FR-FR) 



207 H2R=GR-SR 



208 H2I=GI-SI 



209 GO TO 317 



210 120 FLPS = 



211 ZM = DSQRT(ZMSO) 



212 ZRT2M = DSORT(ZM) 



213 IF (ZR .LT. O.DO) GO TO 125 



214 ZRT2R = DSORT (0.5D0 * (ZR + ZM)) 



215 ZRT2I = ZI / (ZRT2R + ZRT2R) 



216 Z32R = ZR*ZRT2R - ZI*ZRT2I 



217 Z32I = ZR+ZRT2I + ZI*ZRT2R 



218 GO TO 130 



219 125 ZRT2I = DSORT (0.5D0 » (ZM - ZR)) 



220 IF (ZI .LT. O.DO) ZRT2I = -ZRT2I 



221 ZRT2R = ZI / (ZRT2I + ZRT2I) 



222 Z32R = ZR*ZRT2R - ZI+ZRT2I 



223 Z32I = ZR*ZRT2I + ZI*ZRT2R 



224 ZM1R = DABS(Z32I) 



225 IF (ZM1R .LT. TLIM) GO TO 130 



226 R = (TLIM / ZM1R) 



227 Z32R = Z32R * R 



68 



