      SUBROUTINE FF(SRHO,OMEGA)                                         00000100
      COMMON/RATIO/ITT,RAT(500,3),ERR(500,3)                            00000200
      DIMENSION BETA(3),SRHO(3),OMEGA(3),B(9),A(3,3),FSUM(3)            00000300
      DIMENSION LW(9),M(9),F(9)                                         00000400
      DATA F/.577,.0387,.00104,.3965,.00753,.154,.212,.111,.0524/       00000500
      WRITE(6,2000)                                                     00000600
 2000 FORMAT(' -----------------------------')                          00000700
      N=3                                                               00000800
      DO 10 I=1,3                                                       00000900
      FSUM(I)=0.                                                        00001000
      DO 10 J=1,3                                                       00001100
      A(I,J)=F(3*(I-1)+J)*OMEGA(J)                                      00001200
10    FSUM(I)=FSUM(I)+A(I,J)*SRHO(J)                                    00001300
      DO 30 I=1,3                                                       00001400
      DO 30 J=1,3                                                       00001500
      L=3*(I-1)+J                                                       00001600
      SUN=0.                                                            00001700
      DO 20 K=1,3                                                       00001800
      LL=3*(I-1)+K                                                      00001900
   20 SUN=SUN+A(J,K)*F(LL)                                              00002000
   30 B(L)=SUN                                                          00002100
      CALL MINV(B,N,D,LW,M)                                             00002200
      DO 40 I=1,3                                                       00002300
      BETA(I)=0.                                                        00002400
      DO 40 J=1,3                                                       00002500
      L=3*(J-1)+I                                                       00002600
   40 BETA(I)=BETA(I)+B(L)*FSUM(J)                                      00002700
      WRITE(6,2001)BETA                                                 00002800
      B1=SQRT(B(1))                                                     00002900
      B5=SQRT(B(5))                                                     00003000
      B9=SQRT(B(9))                                                     00003100
      WRITE(6,2002)B1,B5,B9                                             00003200
      WRITE(6,2003)                                                     00003300
 2003 FORMAT(7X,' AL/SI ',5X,' MG/SI ',5X,' MG/AL ')                    00003400
      RAT(ITT,1)=BETA(2)/BETA(1)                                        00003500
      RAT(ITT,2)=BETA(3)/BETA(1)                                        00003600
      RAT(ITT,3)=BETA(3)/BETA(2)                                        00003700
      WRITE(6,2004)(RAT(ITT,J),J=1,3)                                   00003800
      ERR(ITT,1)=SQRT((B1/BETA(1))**2+(B5/BETA(2))**2)*RAT(ITT,1)       00003900
      ERR(ITT,2)=SQRT((B1/BETA(1))**2+(B9/BETA(3))**2)*RAT(ITT,2)       00004000
      ERR(ITT,3)=SQRT((B5/BETA(2))**2+(B9/BETA(3))**2)*RAT(ITT,3)       00004100
      WRITE(6,2005)(ERR(ITT,J),J=1,3)                                   00004200
      WRITE(6,2000)                                                     00004300
      WRITE(6,2006)                                                     00004400
      B2=B(2)**2/(B(1)*B(5))                                            00004500
      B3=B(3)**2/(B(1)*B(9))                                            00004600
      B6=B(6)**2/(B(5)*B(9))                                            00004700
      WRITE(6,2007)B2                                                   00004800
      WRITE(6,2007)B3,B6                                                00004900
 2004 FORMAT(1X,'RATIOS',1X, 3(F8.6,4X))                                00005000
 2005 FORMAT(1X,'ERRORS',1X,3(F8.6,4X))                                 00005100
 2007 FORMAT(1X,3(2X,E12.6))                                            00005200
 2001 FORMAT(' BETAS -- SI,AL,MG --  ',3(2X,E12.6))                     00005300
 2002 FORMAT(' ERRORS-- SI,AL,MG --  ',3(2X,E12.6))                     00005400
 2006 FORMAT(' INTERFERENCE MATRIX')                                    00005500
      RETURN                                                            00005600
      END                                                               00005700
