File: TDMU.FC of Tape: Various/ETH/fc2
(.FC,.FO Focal format converted to 'WRITE' listing)
C-OS/8 FOCAL, 1972 01.01 C DINGLE FIT: TDMU.FC 01.05 E 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 AMP TAU DAMP DTAU" 01.50 S J=J+1;A !,X(J);I (-X(J))1.6;S N=J-1;G 1.7 01.60 F I=1,NF;T !,%3,:7, NU(I);S L=NF*J+I;D 1.65 01.62 G 1.5 01.65 D 10 01.70 F K=1,NF;T !!"FREQUENZ",%4.01,NU(K),!,"AMPFIT:";D 1.75; 01.72 Q 01.75 D 5;D 7;D 9;T "TAUFIT:";D 6;D 9 02.01 C SUMS 02.10 S A=0;S B=0;S C=0;S D=0;S E=0 02.20 S DE=DY(1)^2;F I=1,N;S P(I)=DE/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 UND PRINT 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-1));S MA=M*FSQT(C/G);S MB=M*FSQT(NP/G) 03.40 S D1=MA/AA;S D2=MB/BB;S DT=TD*FSQT(D1*D1+D2*D2) 04.10 S H=NP*E-B*B;S RR=BB*BB*G/H;S RR=FSQT(RR);T %5.04 !,"RR =", RR 05.01 C AMPFIT 05.05 F I=1,N;S M=NF*I+K;D 5.1 05.07 G 5.2 05.10 S Y(I)=-FLOG(FCOM(M)/X(I));S DY(I)=FCOM(128+M)/FCOM(M) 05.20 S KK=1.47E5/HI;D 2;D 3 06.01 C TAUFIT 06.10 F I=1,N;S M=NF*I+K;S Y(I)=FCOM(64+M);S DY(I)=FCOM(192+M) 06.20 S KK=DK/.838;D 2;D 3 07.01 C TD AUS AMPFIT 07.05 S KK=DK/.838;S TM=0;F I=1,N;D 7.1 07.07 G 7.2 07.10 S M=NF*I+K;S TD(I)=FCOM(64+M)/(KK*MU)-X(I);S TM=TM+TD(I) 07.20 S TD=TM/N;S TQ=0 07.30 F I=1,N;S TQ=TQ+(TD-TD(I))^2 07.40 S DT=FSQT(TQ/(N-1)) 09.01 C PRINT 09.10 T " TD="%4.03,TD," (+/-",DT,")" 09.20 T " MU="%4.03,MU," (+/-",MB/KK,")" 09.30 T " RMSD=",%3.01,100*FSQT(QA/(N)),"%" 09.40 D 4;T ! 10.10 A :13,DU;S Z=FCOM(L,DU) 10.20 A :20,DU;S Z=FCOM(64+L,DU) 10.30 A :26,DU;S Z=FCOM(128+L,DU) 10.40 A :32,DU;S Z=FCOM(192+L,DU) 10.50 S Z=FCOM(256+L,PS) 10.60 S Z=FX(0,500)