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

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

---

C-FOU   .FC MAY-76

01.10 A ?N?
01.20 S NR=N/2
01.30 A ?TA?
01.40 D 2
01.50 G 1.3

02.10 S TA=TA*2*PI/N
02.20 S ME=FEXP(-NR*TA);S PE=FEXP(NR*TA)
02.30 S TH=(PE-ME)/(PE+ME)
02.40 S CE=NR/TH-1/TA
02.50 T %6.04,CE*PI/1024,!

03.20 S ME=FEXP(-TA);S PE=FEXP(TA)
03.30 S TH=(PE-ME)/(PE+ME)

04.10 A !,?A?
04.15 S X1=A*3
04.20 I (FABS(A),>,.5)4.25;S X1=1/(FSGN(A)-A)
04.25 S D1=1E10
04.30 S TA=X1;D 3;S X2=TH/(1-A*TH)
04.31 S D2=X2-X1
04.32 I (FABS(D1),>,1E-6)4.33;S Q=D2/D1;S P=Q/(1-Q)
04.33 S X2=X2+P*D2
04.40 I (FABS(X2-X1),<,1E-3)4.5;T X2;G 4.1
04.50 S X1=X2;S D1=D2
04.70 G 4.3

05.10 S Z=(1+X)/(1-X)
05.20 S Y=FLOG(Z)/2

---