File: FOULEG.FC of Tape: Various/ETH/fc-files

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

---

C-FOULEG.FC NOV-77

05.01 C FOUR
05.10 A ?SI?,?SF?;S X=(SF-SI)/16;I (X-FITR(X)-.5)5.14;S X=X+.5
05.14 S V=FITR(X)/16;S SF=SI+V*256;S SI=SI+V*128
05.16 T %3,"'"V*32"'"
05.20 D 13.9;D 12
05.30 S Z=FOUR(SI,V);S Z=FX(1,760)

12.01 C FOURSCALE
12.10 S X=194;S U=SI/100;S Y=FITR(U)
12.12 I (-Y)12.14;S Z=FDIS(X+2,-96);G 12.16
12.14 S LG=12.2;G 12.52
12.16 D 12.8;I (-Y)12.18;S Z=FDIS(X+1,-96);G 12.8
12.18 D 12.9;G 12.8
12.20 S X=201;S U=SI/100;S Y=FITR(U);S LG=12.3;G 12.5
12.30 S X=212;S U=SF/100;S Y=FITR(U);S LG=12.34;G 12.5
12.34 S LG=2.01;R
12.50 I (-Y)12.52,12.54
12.52 D 12.9;G 12.56
12.54 S Z=FDIS(X,-96);S X=X+1
12.56 D 12.8;D 12.9;D 12.8;D 12.9;G LG
12.80 S U=(U-Y)*10+1E-4;S Y=FITR(U)
12.90 S Z=FDIS(X,Y-80);S X=X+1
12.92 G 12.2

13.01 C DIS SWITCH
13.40 S Z=FDIS(0,0)
13.50 S Z=FX(0,500)
13.60 S Z=FX(ST,500)
13.70 S Z=FDIS(Z,FABS(FDIS(Z)));S Z=FDIS(Y,-FABS(FDIS(Y)))
13.80 S Z=4;S Y=5;D 13.7
13.90 S Z=5;S Y=4;D 13.7

23.01 C LEGENDRE
23.05 S A(1)=0;S A(2)=0;S A(3)=0;S A(4)=0
23.10 F I=0,N-1;D 24
23.15 S A(2)=A(2)*3;S A(3)=A(3)*5;S A(4)=A(4)*7
23.20 T !,%6.04,"A0=",A(1),"A1=",A(2),"A2=",A(3),"A3=",A(4)

24.01 C SUB OF 23
24.10 S U=FCOM(I)/(N-1)
24.20 S Y=I/(N-1)
24.30 I (-FABS(I*(I-N+1)))24.6;S U=U/2
24.60 S A(1)=A(1)+U
24.70 S A(2)=A(2)+(2*Y-1)*U
24.80 S A(3)=A(3)+(6*Y^2-6*Y+1)*U
24.90 S A(4)=A(4)+(20*Y^3-30*Y^2+12*Y-1)*U

---