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

(Source file text) 

---

	SUBROUTINE RFIT(IMAX,M1,N1,X)
	IMAX1=IMAX+1
	IMAX2=2*IMAX-1
	IMAX3=IMAX+2
	DIMENSION A(10,11)
	DIMENSION X(10)
	DIMENSION SUM(19)
	DIMENSION OG(19)
	DIMENSION S(80,2)
	REWIND 5
C	WRITE(4,70)
70	FORMAT(1X,/,' NR DES MITTLEREN DATENPAARES (X,I3):')
	I=1
C	READ(4,31) M1
C	IF (M1-1) 203,203,202)
C203	M=M1
C	N=M1+2
C	GOTO 204
C202	M=M1-1
C	N=M1+1
C204	CONTINUE
	M=M1
	N=N1

31	FORMAT(1X,I3)

1	READ(5,74) IEND,TT,RR
74	FORMAT(1X,I5,2E16.9)

	IF (IEND) 61,51,61
61	S(I,1) = RR
	S(I,2) = 1./TT
	I=I+1
	GOTO 1
51	DO 36 I=1,IMAX2
36	SUM(I)=0.
	DO 37 I=1,IMAX
37	A(I,IMAX1)=0.
	DO 38 I= M,N
	HELP=ALOG(S(I,1))
	OG(1)=1./(HELP*HELP)
	SUM(1)=SUM(1) + OG(1)
	DO 39 K=2,IMAX2
	J=K-1
	OG(K)=OG(J)*HELP
	SUM(K) = SUM(K)+OG(K)
	IF(K-IMAX3) 42,54,54
42	A(J,IMAX1) =A(J,IMAX1) + S(I,2)*OG(K)
54	GOTO 39
39	CONTINUE
38	CONTINUE
C	GAUSSALGORITHMUS
	DO 53 I=1,IMAX
	DO 53 K=1,IMAX
	J=I+K-1
53	A(I,K)=SUM(J)
	N=IMAX
	N1=IMAX-1
	N2=IMAX1
	DO 3 M=1,N1
	M1=M+1
	IF(A(M,M)) 4,5,4
4	DO 3 I=M1,N
	QUO=A(I,M)/A(M,M)
	DO 3 K=M1,N2
3	A(I,K)=A(I,K)-QUO*A(M,K)
	IF(A(N,N)) 6,5,6
6	X(N)=A(N,N2)/A(N,N)
	DO 8 L=1,N1
	I=N1-L+1
	SU=0.
	I1=I+1
	DO 7 K=I1,N
7	SU=SU+A(I,K)*X(K)
8	X(I)=(A(I,N2)-SU)/A(I,I)
	WRITE(4,30)
30	FORMAT(X,//,X,14HFITPARAMETER ://)
	DO 77 I=1,N
	J=-2+I
	WRITE(4,78) J,X(I)
78	FORMAT(X,2HA(I2,4H) = E20.8)
77	CONTINUE
	GOTO 9
5	WRITE(4,40)
40	FORMAT(20X,32HHAUPTDIAGONALENGLIED GLEICH NULL)
	REWIND 5
9	RETURN
	END

---