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

(Source file text) 

---

C	PROGRAMM FUER THERMOMETERFIT
	DIMENSION A(10)
	LOGICAL GOON
199	WRITE(4,200)
200	FORMAT(/1X,'FIT GEWUENSCHT? (T/F)',T70,$)
	READ(4,210) GOON
210	FORMAT(L1)
	IF (.NOT.GOON) GOTO 99
	WRITE(4,100)
100	FORMAT(2X,' ANZAHL TERME IN DER REIHE',T70,$)
	READ(4,110) IM
110	FORMAT(I2)
	WRITE(4,120)
120	FORMAT(/1X,'ERSTES WERTEPAAR (NR)',T70,$)
	READ(4,110) I1
	WRITE(4,140)
140	FORMAT(/1X,'LETZTES PUNKTEPAAR (NR)',T70,$)
	READ(4,110) I2

	CALL RFIT(IM,I1,I2,A)
	GOTO 199
	WRITE(4,300)
300	FORMAT(/1X,'TABELLE ?',T70,$)
	READ(4,210) GOON
	IF (.NOT.GOON) GOTO 199
	WRITE( 4,400)
400	FORMAT(/1X,'SCHRITTWEITE IN R ?',T70,$)
	READ(4,500) RD
500	FORMAT(E16.8)
	WRITE(4,600)
600	FORMAT(/1X,'ANFANGSWERT FUER R?',T70,$)
	READ(4,500) RA
	WRITE(4,700)
700	FORMAT(/1X,'ABBRUCH BEI DELTA T < ?',T70,$)
	READ(4,500) DELTA

	T1=0
198	H1=ALOG(RA)
	H2=1.
	SUM =A(1)/H1 + A(2)
	DO 800 I=3,IM
	H2=H2*H1
	SUM =SUM+A(I)*H2
800	CONTINUE
	TE=1./SUM
	WRITE(4,900) TE,RA
900	FORMAT(1X,2E16.8)
	IF(ABS(TE-T1).LT.DELTA) GOTO 199
	RA=RA+RD
	T1=TE
	GOTO 198
99	STOP
	END

---