File: FIT.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1
(Source file text)
C FIT FUER KOHLENSTOFFTHERMOMETER INTEGER IMAX WRITE(4,95) 95 FORMAT(1X,/,' ANZAHL TERME (FORMAT(X,I3)):') READ(4,31) IMAX 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) REWIND 5 WRITE(4,70) 70 FORMAT(1X,/,' NR DES ERSTEN ZU VERWENDENDEN DATENPAARES (X,I3):') I=1 READ(4,31) M WRITE(4,71) 71 FORMAT(1X,/,' NR DES LETZTEN DATENPAARES (FORMAT(X,I3)):') READ(4,31) N 31 FORMAT(1X,I3) 1 READ(5,74) IEND,T,R 74 FORMAT(1X,I5,2E16.9) IF (IEND) 61,51,61 61 S(I,1) = R S(I,2) = 1./T 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) 9 STOP END