File: DRIFT.FC of Tape: Various/ETH/fc1
(.FC,.FO Focal format converted to 'WRITE' listing) 

C-OS/8 FOCAL, 1972

01.01 C PROG. FLUX DRIFT: D R I F T
01.02 C AUSWERTUNG DER FLUX NOISE MESSUNGEN 27.5.75
01.10 D 2;D 3;D 4
01.90 Q

02.01 C SUMMEN
02.05 E
02.10 A ?I0,N?,!
02.20 F I=I0,I0+N-1;S A=A+I;S B=B+FCOM(I);S C=C+I*I
02.30 F I=I0,I0+N-1;S D=D+I*FCOM(I);S E=E+FCOM(I)^2
02.40 S G=N*C-A*A
02.50 S AA=(B*C-A*D)/G;S BB=(N*D-A*B)/G

03.01 C CHI QUADRAT
03.10 S Q=0;F I=I0,I0+N-1;S Q=Q+(FCOM(I)-I*BB-AA)^2
03.20 S VA=FSQT(Q/(N-2))
03.30 S H=N*E-B*B;S RR=FSQT(BB*BB*G/H)

04.01 C DRUCK
04.10 T !,%5.03,"AA=",AA,"   BB=",BB,"   VA=",VA,"   RR=",RR
04.20 T !

05.10 F I=0,255;S Z=FCOM(I,(1/A^2)-1/(A+B*I)^2);D 13.5

13.50 S Z=FX(0,500)
13.60 S Z=FX(1,500)