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