File PRIN1.RM

Directory of image this file is from
This file as a plain text file

	SUBROUTINE PRIN
        COMMON IBUFF(255),DUMI(85),MANU,IBLK,CP,CV,CF
        COMMON/BLK1/VOLN(250),VOLX(250),FLON(260),FLOX(260),PRSN(250)
        COMMON/BLK2/ACCN(250),ACCX(250),PRSX(250)
        COMMON/BLK3/INEG,IPOS,IDIM,ITMP,IDIM1,IPAS,ISTOP,IREC
        COMMON/BLK4/HOLD(30),K,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 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
C(((((
	SEM=SQRT(SSQ*AMAT(I,I)/AM)
	TV=AMAT(I,I)/SEM
	WRITE(3,999)I,SEM,TV
  999	FORMAT(1X,I2,2E15.6)
C)))))
	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
	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,/))
	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	DUMI(I)=0.
	AM=AM-1.
	X=AM
	DO 140 L=1,4
	DO 130 I=1,4
  130	DUMI(L)=DUMI(L)+AMAT(I,L)*A12(I)
  140	CONTINUE
	RETURN
	END



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search