File: TMF.FT of Disk: Disks/MyPDP/m8-backup-rka1-rkb1

(Source file text) 

---

C PGM TMF ENTHALPIEFIT
C FIT FUER DATEN AUS MESSUNG DER SPEZ. WAERME UM THERMOMETER ZU EICHEN

	DIMENSION R1(10),R2(10),Q(10),A(10,6),X(6),B(10),S(6),DN(6)
	INTEGER FRAGE
	READ(4,10)R1,R2,Q
 10	FORMAT(3F12.0)
	WRITE(3,20)
 20	FORMAT(1H1,1X,'EINGEGEBENE WERTE:'/
     $		1X,'    R1          R2          Q')
	DO 30 I=1,10
 30	WRITE(4,40)R1(I),R2(I),Q(I)
 40	FORMAT(1X,3F12.2)
	S(1)=0.25
	S(2)=0.01
	S(3)=0.5
	S(4)=-0.1
	S(5)=0.05
	S(6)=0.0005
	WRITE(3,45)
 45	FORMAT(1X,'NAEHERUNGSWERTE:')
	DO 48 I=1,6
	WRITE(3,50)I,S(I)
 50	FORMAT(1X,'X(',I1,') =',E14.8)
 48 	CONTINUE
	DO 60 I=1,10
	LN1=ALOG(R1(I))
	LN2=ALOG(R2(I))
	T1=1./(S(3)/LN1+S(4)+S(5)*LN1+S(6)*LN1*LN1)
	T2=1./(S(3)/LN2+S(4)+S(5)*LN2+S(6)*LN2*LN2)
	B(I)=1.-S(1)/(2.*Q(I))*(T2*T2-T1*T1)
	B(I)=B(I)-S(2)/(4.*Q(I))*(T2**4-T1**4)
	A(I,1)=.5/Q(I)*(T2*T2-T1*T1)
	A(I,2)=.25/Q(I)*(T2**4-T1**4)
	A(I,3)=-S(1)/Q(I)*(T2**3/LN2-T1**3/LN1)
	A(I,3)=-S(2)/Q(I)*(T2**5/LN2-T1**5/LN1)+A(I,3)
	A(I,4)=-S(1)/Q(I)*(T2**3-T1**3)
	A(I,4)=-S(2)/Q(I)*(T2**5-T1**5)+A(I,4)
	A(I,5)=-S(1)/Q(I)*(T2**3*LN2-T1**3*LN1)
	A(I,5)=-S(2)/Q(I)*(T2**5*LN2-T1**5*LNI)+A(I,5)
	A(I,6)=-S(1)/Q(I)*(T2**3*LN2**2-T1**3*LN1**2)
	A(I,6)=-S(2)/Q(I)*(T2**5*LN2**2-T1**5*LN1**2)+A(I,6)
 60	CONTINUE
	CALL HLS(A,B,6,1,IER,AUX,IPIV,ERS,X)
	WRITE(3,70)
 70	FORMAT(///,1X'BERECHNETE WERTE UND FEHLER IN PROZENT')
	DO 75 I=1,6
	DN(I)=(X(I)-S(I))*1
	WRITE(3,77)I,S(I),DN(I)
 77	FORMAT(1X,I2,E14.8,E6.3)
 75 	CONTINUE
	WRITE(4,80)
 80	FORMAT(1X,'NEUER FIT, FALLS JA 1 EINGEBEN SONST 0:',T70,$)
	READ(4,90) FRAGE
 90	FORMAT(I1)
	DO 95 I=1,6
 95	S(I)=X(I)
	IF (FRAGE.EQ.1) GOTO 48
	END

---