File: DINGLE.FC of Tape: Various/ETH/prog1
(.FC,.FO Focal format converted to 'WRITE' listing)
C-DINGLE.FC MAY-75 01.01 C DINGLE FIT !!! F O S 8 !!! 01.02 C ALTES TD MIT KT, GEAENDERT 01.03 C VARIANZ DER FIT WERTE BERUECKSICHTIGT 01.04 C CHI^2 SIGNIFIKANZ AUSSAGE 01.05 E 01.06 G 21.1 01.10 T "FREQ:";S I=0 01.20 S I=I+1;A NU(I);I (-NU(I))1.2;S NF=I-1;G 1.3 01.30 A !?HI,DK? 01.40 T !" TEMP NU NUEX AMP TAU DAMP DTAU A0" 01.50 S J=J+1;A !,X(J);I (-X(J))1.6;S N=J-1;G 1.65 01.60 F I=1,NF;T !,%3,:7, NU(I);S L=NF*J+I;D 10 01.62 G 1.5 01.65 F I=1,NF;D 15 01.70 F K=1,NF;D 6;D 9 01.72 O I TTY:,E;Q 02.01 C SUMS 02.10 S A=0;S B=0;S C=0;S D=0;S E=0 02.20 F I=1,N;S P(I)=1/DY(I)^2 02.40 F I=1,N;S A=A+P(I)*X(I);S B=B+P(I)*Y(I);S C=C+P(I)*X(I)*X(I) 02.50 F I=1,N;S D=D+P(I)*X(I)*Y(I);S E=E+P(I)*Y(I)*Y(I) 02.55 S NP=0;F I=1,N;S NP=NP+P(I) 02.60 S G=NP*C-A*A 02.70 S AA=(B*C-A*D)/G;S BB=(NP*D-A*B)/G 02.80 S MU=BB/KK;S TD=AA/BB 03.01 C FEHLER MIT CHI-QUADRAT 03.05 I (2-N)3.1;S DT=0;S MB=0;S QA=0;R 03.10 F I=1,N;S V(I)=(Y(I)-X(I)*BB-AA) 03.20 S Q=0;S QA=0;F I=1,N;S Q=Q+V(I)*V(I)*P(I);S QA=QA+(V(I)/Y(I))^2 03.30 S M=FSQT(Q/(N-2));S MA=M*FSQT(C/G);S MB=M*FSQT(NP/G) 03.40 S D1=MA/AA;S D2=MB/BB;S DT=FABS[TD*FSQT(D1*D1+D2*D2)] 03.50 S H=NP*E-B*B;S RR=FSQT(BB^2*G/H) 04.01 C FEHLER AUS VARIANZ DER FIT WERTE 04.10 S SA=C/G;S SB=NP/G;S SC=A/G 04.20 S DT(1)=FABS[TD*FSQT(SA/AA^2+SB/BB^2-2*SC/(AA*BB))] 06.01 C TAUFIT 06.10 F I=1,N;S M=NF*I+K;S Y(I)=FX(64+M);S DY(I)=FX(192+M) 06.20 S KK=DK/.838;D 2;D 3;D 4 09.01 C PRINT 09.10 T !!,"NU",%5.02,NE(K) 09.20 T !," TD =",%4.03,TD," +/-",DT,:30,"F =",NE(K)/DK/51.2 09.30 T :50,"CHI^2 =",M,!,:18,DT(1),:30,"P =",512*DK/NE(K) 09.40 T :50,"RMSD =",100*FSQT(QA/N)," %" 09.50 T !," MU =",MU," +/-",MB/KK,:30,"A0 =",A0(K) 09.60 T :50,"RR =",%5.04,RR, 09.70 T !,:18,%4.03,SB/KK,:30,"MU*P =",MU*512*DK/NE(K),! 10.01 C FIT WERTE SPEICHERN 10.05 A NE;S Z=FX(32+L,NE) 10.10 A AF; 10.20 A TA;S Z=FX(64+L,TA) 10.21 D 15.4;S Z=FX(L,AM) 10.22 S Z1=FEXP(-TA*PI/(256*DK*1E-7*HI)) 10.24 S Z=FX(96+L,AM*FSQT(HI*1E-8)/(X(J)*Z1)) 10.30 A DU;S Z=FX(128+L,DU*AM/AF) 10.40 A DU;S Z=FX(192+L,DU) 10.50 T %5.03,FX(96+L) 12.10 T !!!"FREQUENZ",%5.02,NE(K) 12.20 T " PERIODE",512*DK/NE(K) 12.30 T " K.TERM. (A0)",A0(K),"(+/-",KV(K),")",! 14.01 C TAU AUS <A0> MIT KONST. AM 14.10 S TA=FLOG(A0(I)*X(J)/[AM*FSQT(HI*1E-8)])*256*DK*1E-7*HI/PI 15.01 C TAU MIT <A0> MODIFIZIEREN 15.05 S NE(I)=0 15.10 S KT(I)=0;S KG=0;S KV(I)=0;F J=1,N;S L=NF*J+I;D 16 15.15 S NE(I)=NE(I)/N 15.20 S A0(I)=KT(I)/KG;S KV(I)=FSQT((KV(I)/KG-A0(I)^2)/(N-1)) 15.30 F J=1,N;S L=NF*J+I;D 15.8;D 14;D 15.5 15.35 R 15.40 S AM=AF*TA*PI/(1-FEXP(-TA*PI)) 15.50 S Z=FX(64+L,TA) 15.80 S AM=FX(L);S T1=FX(64+L);S AF=AM*(1-FEXP(-T1*PI))/(T1*PI) 16.01 C A0 MITTELWERT <A0> 16.05 S KP=FX(128+L)/FX(L)+FX(192+L)/FX(64+L);S KP=1/KP 16.10 S KT(I)=KT(I)+FX(96+L)*KP 16.20 S KV(I)=KV(I)+FX(96+L)^2*KP 16.30 S KG=KG+KP 16.40 S NE(I)=NE(I)+FX(32+L) 21.01 C GEFITTETE PARAMS AUF FILE "TDTEMP" 21.10 A "MESS#:",M 21.20 A "RUN'S VON:",D1," BIS:",D2,! 21.30 O O DKA1:TDTEMP 21.40 S D=D1;D 23;D 24;D 25 21.50 F D=D1+1,D2;D 23;D 25 21.60 S X=0;T X 21.70 O C 21.80 E 21.90 O I DKA1:TDTEMP,E;G 1.1 23.01 C DATEN VON FILE 23.10 O I DKA1:M(M)D(D) 23.15 F I=1,6;S MD(I)=FIN() 23.20 A D,N;F I=1,N;A X 23.30 A X;I (X-20)23.35,23.4,23.35 23.35 T "FEHLER: X=",X;D 23.8;Q 23.40 A DK,HR,HI,HF,MI,PS,PT,RT,TM,TV,SH,MP;F I=1,8;A X 23.50 F I=1,4;A NU(I),PH(I),AF(I),TA(I) 23.60 F I=1,4;A DN(I),DP(I),DA(I),DT(I) 23.70 O O TTY:,E 23.80 F I=1,6;S Z=FOUT(MD(I)) 23.90 I (-DA(1))23.95 23.92 T " SCHLECHT";S SW=1 23.95 T "$"!;O O " 24.01 C TDTEMP 1.TEIL 24.20 S NN=0;T %3.0 24.30 S NN=NN+1;T FITR(NU(NN));I (-NU(NN))24.3 24.40 S NN=NN-1;T %5.03,HI,DK 25.01 C TDTEMP 2.TEIL 25.05 I (SW),25.1;S SW=0;R 25.10 T %5.04,TM 25.20 F I=1,NN;T NU(I),AF(I)*PS*100,TA(I),DA(I)*PS*100,DT(I)