57 IF (ZR .GT. -7.5) GO TO 40 



S3 S = CAy ( i + 1 ) 



59 T = CAYI(I + 1 ) 



60 CON = (G(I+1) * S + GI(I+1) * T) / (S**2 + T**2) 



61 CON = 1 . / C0N**2 



62 DO 36 J = 1 ,20 



63 TEMP = 5**2 + T**2 



64 TEMPI = (LAMBDI * S - LAMBDA * T) / TEMP 



65 TEMP = (LAMBDA * S + LAMBDI * T) / TEMP 



66 ZR = ((1. + TEMP) * ( 1. - TEMP) + TEMPI**2) * C0N(I+1) 



67 R = ZR / -7.5 



68 IF (DABS(R-1.) . LT . 1 . D-3 ) GO TO 41 



69 S = LAMBDA + (S - LAMBDA) / R 



70 36 CONTINUE 



71 41 ZI = -2. * TEMPI ♦ TEMP • C0N(I+1) 



72 ZT( 1+1 ) = ZR 



73 ZTI(I+1 ) = ZI 

 74 



75 40 CONTINUE 



76 CALL HANKEL(ZR,ZI ,0) 



77 ZT(I+1 ) = ZR 



78 7TTfT+ll=7T 



79 RATI0(2*I+1) = EXPONT 



80 IF (I .NE. 0) GO TO 45 



81 A(1 ,3) = H2R * RHO( 1 ) 



82 0(1,3)= H2I * RHO( 1 ) 



83 A( 1 ,4) = H1R * RHO( 1 ) 



84 0( 1 ,4) = HI I * RHO( 1 ) 



85 GO TO 10 



86 45 CONTINUE 



87 A{2+I,3) =-H2R * RH0(I+1) 



88 0(2*1,3) =-H2I * RH0(I+1) 



89 A(2*I,4) =-H1R + RH0(I+1) 



90 0(2*1,4) =-H1 I * RH0(I+1) 



91 A(2*I+1,3) =-H2PR * G(I+1) + H2PI • GI(I+1) 



92 0(2*1+1,3) =-H2PI * G(I+1) - H2PR * GI(I+1) 



93 A(2*I+1,4) =-H1PR * G( 1+1 ) + H1PI ♦ GI(I+1) 



94 0(2*1+1,4) =-H1PI * G(I+1) - H1PR * GI{I+1) 



95 10 CONTINUE 



96 A( 2*N-2 ,4) = 0. 



97 0( 2*N-2 ,4) = 0. 



98 A( 2*N-1 ,4) = 0. 



99 RETURN 

 100 END 



60 



