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

(Source file text) 

---

	SUBROUTINE TEMPFI(IMAX,ICENT,R,T)
	COMMON 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(50,2)
	M1=ICENT
	IF (M1-1) 203,203,202)
203	M=M1
	N=M1+2
	GOTO 204
202	M=M1-1
	N=M1+1
204	CONTINUE

	DO 61 I=M,N
	S(I,1) = R(I)
61	S(I,2) = 1./T(I)
	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

---