SUBROUTINE PRIN COMMON IBUFF(255),DUMI(85),MANU,IBLK,CP,CV,CF COMMON/BLK1/VOLN(100),VOLX(100),FLON(110),FLOX(110),PRSN(100) COMMON/BLK2/ACCN(100),ACCX(100),PRSX(100) COMMON/BLK3/INEG,IPOS,IDIM,ITMP,IDIM1,IPAS,ISTOP,IREC COMMON/BLK4/HOLD(30),K,BUFF(1400),X,Y,APOS,J COMMON/BLK5/AMAT(5,5),COEF(4),COV(5,5),COR(5,5),A21(4) COMMON/BLK6/SD(5),A12(4),CEPT,DET,EPS,SX(5),SXY(5),SY DIMENSION T(5,5),CEPT1(4),SE(4),FV(4),V(4),CV2(4,4),CR2(4,4) AM=X A22=APOS TEM=0. DO 10 I=1,4 10 TEM=TEM+COEF(I)*A12(I) DFNU=AM-4.-1. SSQ=(A22-TEM)/DFNU DO 40 I=1,4 DO 30 L=1,4 CV2(I,L)=SSQ*AMAT(I,L) 30 CR2(I,L)=AMAT(I,L)/SQRT(AMAT(I,I)*AMAT(L,L)) 40 CONTINUE CEPT=0 DO 110 I=1,4 110 CEPT=CEPT+COEF(I)*(SX(I)/AM) CEPT=SY/AM-CEPT DO 120 I=1,4 CEPT1(I)=0. 120 T(I,1)=0. AM=AM-1. C .....A11 INVERSE TIMES A12 DO 140 L=1,4 DO 130 I=1,4 130 T(L,1)=T(L,1)+AMAT(I,L)*A12(I) 140 CONTINUE BDD=0. DO 150 I=1,4 150 BDD=BDD+T(I,1)*A21(I) WRITE(3,997)BDD,(A21(I),I=1,4) 997 FORMAT(' .......BDD,A21='5E15.6) CORCOF=SQRT(BDD/A22) BDD=A22-BDD DF=AM-4 C .....SS=BDD..... C .....AVSQ=MS..... AVSQ=BDD/DF C .....RDF=4..... RSS=A22-BDD RMS=RSS/4. F=RMS/AVSQ SEE=SQRT(AVSQ) CRC=SQRT(RSS/A22) DO 160 I=1,4 SE(I)=SEE*AMAT(I,I)**.5 160 FV(I)=(COEF(I)/SE(I))**2 C ......PRINT STAT PARAMS.....///// WRITE(3,165)(T(L,1),L=1,4),A22 165 FORMAT(' B(I,D)=',5E15.6) WRITE(3,170)(COEF(I),I=1,4),(SE(I),I=1,4), 2 (FV(I),I=1,4), 1 (SD(I),I=1,5),DET,CEPT,BDD,CORCOF,DF,AVSQ,BDD WRITE(3,180)RSS,RMS,F,SEE,CRC WRITE(3,190)((COV(I,L),L=1,5),I=1,5),((COR(I,L),L=1,5),I=1,5) WRITE(3,210)((CV2(I,L),L=1,4),I=1,4) WRITE(3,210)((CR2(I,L),L=1,4),I=1,4) 210 FORMAT(///4(4E15.6,/)) 170 FORMAT(' 'T10,'COEFS',/,4E15.6,/,T10,'STAND ERR' 1 ,/,4E15.6,/,T10,'FVALUES',/,4E15.6,/, 2 /,T15,'STAND.DEVS',5E15.6,/,T5,'DETERM', 3 1X,E15.6,' INTERCEPT',1X,E15.6,'B(D,D)=',E15.6, 4 'CORREL COEF',1X,E15.6,//,T5,'DEG FREE=',1X,E15.6, 5 5X,'MEAN SQ=',1X,E15.6,4X,'RDF=4',5X,'SS=',1X,E15.6) 180 FORMAT(//,T5,'RSS='1X,E15.6,5X,'RMS=',1X,E15.6,5X,'F=',1X, 1 E15.6,5X,'EST STAND ERR=',1X,E15.6,//,T5,'MULT COR COEF', 2 1X,E15.6) 190 FORMAT(' COVARIANT MAT',/,5(T5,5E15.6,/),/////, 1 ' CORRELATION MAT',/,5(T5,5E15.6,/)) RETURN END