File: LINREG.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1

(Source file text) 

---

	SUBROUTINE LINREG(TC,CC)
C	LINEARE REGRESSION FUER SPEZIELLE FUNKTION
C	F(X)= GAMMA*X+DELTA*X**3.


	DIMENSION TC(100)
	DIMENSION CC(100)
	DIMENSION A(5)
	DO 1000 I=1,5
1000	A(I)=0
	N=0
	I=1
	WRITE(3,160)
100	IF(CC(I)) 10,20,10
10	HELP=CC(I)/TC(I)
	HSQR=TC(I)*TC(I)
	WRITE(3,170) TC(I), HELP,HSQR
	A(1)=A(1)+HSQR*HELP
	A(2)=A(2)+HSQR
	A(3)=A(3)+HELP
	A(4)=A(4)+HSQR*HSQR
	A(5)=A(5)+HELP*HELP
	N=N+1
20	I=I+1
	IF(I-100) 100,100,300


300	AN=FLOAT(N)

	BB=(A(1)-A(2)*A(3)/AN)/(A(4)-A(2)*A(2)/AN)
	AA=A(3)/AN-BB*A(2)/AN
	H1=(A(1)-A(2)*A(3)/AN)
	H2=(A(4)-A(2)*A(2)/AN)*(A(5)-A(3)*A(3)/AN)
	RSQR=H1*H1/H2

	WRITE(4,110)
	WRITE(4,120) N
	WRITE(4,130) AA
	WRITE(4,140) BB
	WRITE(4,150) RSQR

	WRITE(3,110)
	WRITE(3,120) N
	WRITE(3,130) AA
	WRITE(3,140) BB
	WRITE(3,150) RSQR

110	FORMAT(/1X,'RESULTAT:')
120	FORMAT(//1X,'ANZAHL FITPUNKTE :',T50,I3)
130	FORMAT(/1X,'DELTA =',T50,E16.5,' J/GK*K')
140	FORMAT(/1X,'GAMMA=',T50,E16.5,' J/GK**4')
150	FORMAT(/1X,'KORRELATION (RSQR)=',E16.5)
160	FORMAT(/1X,'    T(K)         C/T             T**2 ',//)
170	FORMAT(/1X,E16.4,'  ',E16.4,'  ',E16.4)
	RETURN
	END

---