1 SUBROUTINE RCOEF (K3) 



2 IMPLICIT DOUBLE PRECISION (A-H,0-Z) 



3 COMMON/ INPUT/ Z(12), N, OMEGA, V, VI, GCU(12), GS0(12), 



4 1 CAY(12), LAMBDA, LAMBDI , G(12) 



5 2,RH0(12), GI(12), G S0I(12), CAYI(12) 



6 DIMENSION RR(12), RI(12), RA(12), RT(12), CYS0(12), CYS0I(12) 



7 COMMON/RFFL/ AF(12,200), AG(12,200), BF(12,200), BG(12,200), 



8 2 EIGEN(350). EIGEN I ( 350 ) , BR( 25 , 4) , BI(25,4), 08(12), CBI(12), 



9 3 CAYS(3(12), CAYSQI(12), NN 



10 NM = N - 1 



11 I = K3 



12 IF (I .GT. NM) I = NM 



13 110 J = I + I 



14 K = d + 1 



15 IF (NN .NE. 1 ) GO TO 102 



16 L = I + 1 



17 TEMP = CB(I)**2 + CBI(I)**2 



18 CY = OMEGA * CB(I) / TEMP 



19 CYI = -OMEGA * CBI(I) / TEMP 



20 CYSO(I) = CY**2 - CYI**2 



21 CYSQ! ; I ; - CY + CYI 



22 CYSOI(I) = CYSOKD + CYSOI(I) 



23 102 EL SO = C YSOd) " EIGEN(NN)**2 + EIGENI (NN ) ♦♦2 



24 ELSOI = C YSOI(I) - 2. DO * EIGEN(NN) * EIGENI(NN) 



25 TEMP = ELSO + DSORT (ELS0**2 + ELS0I**2) 



26 IF (TEMP .LE. . DO ) GO TO 107 



27 EL = DSORT (TEMP * . 5D0 ) 



28 ELI = ELSOI / (EL + EL) 



29 103 A = AF( I ,NN)*BR( J,2) - AG ( I , NN ) * BI ( d , 2 ) 



30 ♦ + BF( I ,NN)*BR(d, 1 ) - BG( I ,NN)*BI(d, 1 ) 



31 B = AF( I ,NN)*BI(d,2) + AG ( I , NN ) * BR ( d , 2 ) 



32 * + BF( I ,NN)*BI(d,1 ) + BG( I ,NN)*BR(d, 1 ) 



33 E = AF( I ,NN)*BR(K,2) - AG( I , NN ) *BI ( K , 2 ) 



34 * + BF(I ,NN)*BR(K, 1 ) - BG( I , NN ) * BI ( K , 1 ) 



35 F = AF( I ,NN)*BI(K,2) + AG( I , NN ) *BR ( K , 2 ) 



36 * + BF(I ,NN)*BI(K,1 ) + BG( I , NN ) * BR ( K , 1 ) 



37 C = (F * EL - E + ELI) / (ELSO + ELSOI) 



38 D = -(E * EL + F * ELI) / (ELSO + ELSOI) 



39 TEMP = (A + C)**2 + (B + D)**2 



40 RR(I) = (A**2 - C**2 + B**2 - D**2) / TEMP 



41 RI(I) = -2. DO * (A * D - B * C) / TEMP 



42 10 FORMAT (10D13.5) 



43 RA(I) = 



44 IF (CB(I) .GT. V) GO TO 104 



45 RX = CB( I) / V 



46 RA(I) = ACOS(RX) * 57.296 



47 104 RT(I) = RR(I)**2 + RI(I)**2 



48 RT(I) = 1 .DO / RT( I) 



49 RI(I) = -DATAN2 (RI(I), RR(I)) * 57.296D0 



50 RR(I) = -4.34294D0 * DLOG (RT(I)) 



51 IF (K3 .NE.1 ) GO TO 108 



52 1=1+1 



53 IF (I . LT. N) GO TO 110 



54 105 CONTINUE 



55 PRINT 106, (RR(I), I = 1,NM) 



56 PRINT 106, (RI(I), I = 1,NM) 



63 



