File: EDYFIT.FC of Tape: Various/ETH/fc1

(.FC,.FO Focal format converted to 'WRITE' listing) 

---

C-OS/8 FOCAL, 1972

01.10 G 15.01

15.01 E
15.20 S C=2*PI*PI*1E-10;S C=1/C;S C=FSQT(C)
15.30 S L=5;S N=4;S SW=-1
15.35 O I EDY
15.40 A X;I (X)0,20.13;A PH;
15.50 S PR=FSQT(90-PH)+271/PH;S PT=X*C
15.60 F J=0,N;S ZY=PT*PT+ZY;S Z(J)=Z(J)+PT*PR^J;S A(J+L*L)=Z(J);D 15.7
15.65 G 15.4
15.70 F K=0,N;S A(J+L*K)=A(J+L*K)+PR^(J+K)
15.80 S Q=0;F K=0,N;S Q=Q+X(K)*Z^K
15.90 Q

16.01 C TESTWERTE
16.05 A "D",D,"NU",NU
16.10 O I EDY
16.15 T %4.03
16.20 A !,Z,PH;T %4.02,"PH",PH,%5.03;S R=2*PI*PI*D*D*NU/Z^2;T :13,"R",R
16.30 S Z=FSQT(90-PH)+271/PH
16.32 S Q=-23872+8950*P-350*P^2+13.018*P^3-.1517*P^4
16.34 D 15.8
16.40 S R1=D*D*NU*1E10/Q^2;
16.50 T :25,"R1",R1,:38,"DR/R",100*(R-R1)/R,"%"
16.60 G 16.2

20.01 E
20.11 S N=L-1;
20.12 F J=0,N;T !;F K=0,L;T "A(",%2,(J+1)*10+K+1,")";A A(J+L*K)
20.13 F K=0,N;S R(K)=K+1;S I=-1
20.14 S M=1E-6;T !
20.16 F J=0,N;F K=0,N;D 24
20.17 S R(P)=0
20.18 F K=0,L;S A(P+L*K)=A(P+L*K)/M
20.20 F J=0,N;D 25
20.22 S I=I+1
20.23 I (I-N)20.14,20.26,20.14
20.26 F J=0,N;F K=0,N;D 27
20.28 F K=0,N;T !%1,"X("K,")=",%10.07,X(K)
20.30 O I TTY:,E
20.40 T !!;Q

24.05 I (R(J))0,24.3,24.1
24.10 I (FABS(A(J+L*K))-FABS(M))24.3;
24.20 S M=A(J+L*K)
24.22 S P=J;S Q=K
24.30 R

25.10 I (J-P)25.2,25.4,25.2
25.20 S D=A(J+L*Q)
25.30 F K=0,L;S A(J+L*K)=A(J+L*K)-A(P+L*K)*D
25.40 R

27.10 I (1E-6-FABS(A(J+L*K)))27.2;R
27.20 S X(K)=A(J+L*L)

---