File PRIN.MA (MACREL macro assembler source file)

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(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



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