R = DCBRT(R) 







ZRT2R 



= ZRrSR 



* 



R 



ZRT2I 



= ZRT2I 



♦ 



R 



2RT2M 



= ZRT2M 



* 



R 



R = R 



* R 







ZM = ZM * R 







ZMSO = 



= Z1VI**2 







228 

 223 

 230 

 231 

 232 

 233 

 234 



235 ZR = ZR * R 



236 ZI = ZI * R 



237 130 ZRT4R = DSQRT (0.5D0 * (ZRT2R + ZRT2M)) 



238 ZRT4I = 0.5D0 * ZRT2I / ZRT4R 



239 ZRTM4R = ZRT4R/ZRT2M 



240 ZRTIVI4I = -ZRT4I/ZRT2IVI 



241 IF (ZR .GT. 0. ) GO TO 210 



242 IF (ZM1R .LT. TLIM) GO TO 210 



243 ABK = ABS(K2) 



244 IF (Z32I .GT. 0.) GO TO 205 



245 K1 = K * EXPONT 



246 K2 = K / EXPONT 



247 Z32I = -TLIM 



248 GO TO 220 



249 905 K3 = K * FXPONT 



250 K1 = K / EXPONT 



251 Z32I = TLIM 



252 GO TO 220 



253 210 K2 = C2 * Z32I 



254 S2 = DEXP(K2) 



255 K2 = K+S2 



256 K1 = K/S2 



257 220 THR = FPI12 - C2 * Z32R 



258 STHR =DSIN(THR) 



259 CTHR =DCOS(THR) 



260 STP = -C3*CTHR +0.5*STHR 



261 CTP = -C3*STHR -0.5*CTHR 



262 TEMP = DABS (Z32R) 



263 TEM1 = DABS (Z32I) 



264 IF (TEMP . LT . TEM1 ) TEMP = TEM1 



265 230 DO 235 ML = 1 , LAS 



266 IF (TEMP . GT . ZMLA5(ML)) GO TO 250 



267 235 CONTINUE 



268 250 CONTINUE 



269 YR = Z32I 



270 YI = -Z32R 



271 CALL CFR (YR, YI, F2R, F2 I , C4 , C5, ML) 



272 CPR = F2R 



273 CPI = F2I 



274 ST0RE3=K2*(ZRTM4R*F2R-ZRTM4I*F2I) 



275 ST0RE4=K2*(ZRTM4I*F2R+ZRTM4R*F2I) 



276 H22R =ST0RE3*CTKR-ST0RE4*STHR 



277 H22I =STORE3*STHR+STORE4*CTHR 



278 IF (ZR) 280,270,270 



279 270 FLQUAD =0 



280 GO TO 300 



281 280 IF (ZI) 290,310,310 



282 290 FLQUAD = 1 



283 300 H2R = H22R 



284 H2I = H22I 



69 



