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

(Source file text) 

---

C	FITPROGRAMM

	LOGICAL GOON
	DIMENSION A(3,4)
	DIMENSION X(40)
	DIMENSION Y(40)
	DIMENSION XX(3)
	N=0
	DO 100 I=1,3
	DO 100 K=1,4
100	A(I,K) = 0
	DO 200 I=1,3
200	X(I)=0
300	WRITE(4,1000)
	READ(4,1100) GOON
	IF(.NOT.GOON) GOTO 400
	N=N+1
	WRITE(4,1200)
	READ(4,1300) X(N)
	WRITE(4,1400)
	READ(4,1300) Y(N)
	GOTO 300
400	DO 500 I=1,N
	H=X(I)*X(I)
	A(1,4)=A(1,4) +Y(I)
	A(2,4)=A(2,4)+X(I)*Y(I)
	A(3,4)=A(3,4)+H*X(I)*Y(I)
500	CONTINUE
	DO 600 I=1,N
	X2=X(I)*X(I)
	X4=X2*X2
	A(1,2)=A(1,2)+X2
	A(1,3)=A(1,3)+X4
	A(2,1) =A(2,1)+ X(I)
	A(2,2)=A(2,2)+X2*X(I)
	A(2,3)=A(2,3)+X4*X(I)
	A(3,1)=A(2,2)
	A(3,2)=A(2,3)
	A(3,3)=A(3,3)+X4*X2*X(I)
600	CONTINUE
	A(1,1)=FLOAT(N)
	N=3
	N2=N+1

	N1=N-1
	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	XX(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)*XX(K)
8	XX(I)=(A(I,N2)-SU)/A(I,I)
	WRITE(4,1600)
	DO 90 I=1,N
90	WRITE(4,3000) I,XX(I)
	GOTO 9
5	WRITE(4,4000)

1000	FORMAT(/1X,'PUNKTEPAAR? (T/F):',T40,$)
1100	FORMAT(L1)
1200	FORMAT(/1X,'X= ',T8,$)
1300	FORMAT(E16.8)
1400 	FORMAT(/1X,'Y= ',T8,$)
1600	FORMAT(//1X,'FITPARAMETER :',//)
3000	FORMAT(1X,' X(',I2,')  =',E18.5)
4000	FORMAT(/1X,30HHAUPTDIAGONALGLIED GLEICH NULL)
9	STOP
	END
	

6

---