File: UNSHIF.FC of Tape: Various/ETH/fc2

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

---

C-OS/8 FOCAL, 1972

14.10 I (N-64)14.2,14.6,14.2
14.20 I (N-128)14.3,14.7,14.3
14.30 I (N-256)14.4,14.8,14.4
14.40 T !"UNSHIF ACCEPTS ONLY N=64,128,256!"!;L C DATA
14.60 S A(0)=63;S A(3)=2016;S A(6)=85344
14.65 S A(4)=85344;S A(7)=4064256;S A(8)=.206447E9;G 14.9
14.70 S A(0)=127;S A(3)=8126;S A(6)=690880
14.75 S A(4)=.69088E6;S A(7)=.660644E8;S A(8)=.6738429E10;G 14.9
14.80 S A(0)=255;S A(3)=32640;S A(6)=.555968E7
14.85 S A(4)=.555968E7;S A(7)=.10653696E10;S A(8)=.2177604E12
14.90 S A(1)=A(3);S A(2)=A(4);S A(5)=A(7)

15.10 F J=0,2;S A(J+9)=0
15.20 F I=0,N-1;S Z=FCOM(I);D 15.3
15.25 G 15.6
15.30 F J=0,2;S A(J+9)=A(J+9)+Z*I^J
15.60 D 20
15.70 R

20.13 F K=0,2;S R(K)=K+1;S I=-1
20.14 S MU=1E-6
20.16 F J=0,2;F K=0,2;D 24
20.17 S R(P)=0
20.18 F K=0,3;S A(P+3*K)=A(P+3*K)/MU
20.20 F J=0,2;D 25
20.22 S I=I+1
20.23 I (I-2)20.14,20.26,20.14
20.26 F J=0,2;F K=0,2;D 27
20.40 R

24.05 I (R(J))0,24.3,24.1
24.10 I (FABS(A(J+3*K))-FABS(MU))24.3;
24.20 S MU=A(J+3*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 DU=A(J+3*Q)
25.30 F K=0,3;S A(J+3*K)=A(J+3*K)-A(P+3*K)*DU
25.40 R

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

---