@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C-FOCAL , 8/68

01.05 S MR=100;S ACC=0.00001
01.10 A !!"ORDER",N;I (N-2) 1.1
01.20 T !!"COEFFICIENTS",!,"X^N  REAL  IMAG"
01.25 F I=0,N;T ! %2,N-I,"  ";A CR(I),CI(I)
01.30 T !!"   ITERATIONS          ROOT
01.31 T "                   FUNCTION VALUE"
01.32 T !"          "
01.33 T "     REAL PART   IMAG. PART
01.34 D 1.33
01.40 S RN=1;S XR=1.27;S XI=0.73
01.50 S GR=CR(0);S GI=CI(0);F I=1,N;D 7
01.60 S DR=N*CR(0);S DI=N*CI(0);F I=1,N-1;D 8
01.70 S R1=GR;S R2=DR;S I1=GI;S I2=DI;D 6
01.75 S YR=XR-RE;S YI=XI-IM;S ER=FSQT(YR*YR+YI*YI)
01.80 S ER=FABS((YR-XR)/ER)+FABS((YI-XI)/ER)
01.85 I (ER-AC) 1.9;S XR=YR;S XI=YI;S RN=RN+1;I (RN-MR) 1.5;G 1.99
01.90 T !,"   "%3,RN,"    ",%08.05,YR,"  ",YI,"    ",GR,"  ",GI
01.92 I (AC-FABS(YI)) 2.1;
01.95 S N=N-1;D 9;I (1-N) 1.4,10.1;Q
01.99 T !!"CONVERGENCE NOT ACHIEVED AFTER",%3,RN," ITERATIONS";D 1.9; Q

02.10 S RN=1;S XR=YR;S XI=-YI;S N=N-1;D 9;I (1-N) 1.5,10.1; Q

05.10 S RE=R1*R2-I1*I2;S IM=I1*R2+I2*R1;R

06.10 S I2=-I2;D 5;S I2=-I2
06.20 S TM=R2*R2+I2*I2;S RE=RE/TM;S IM=IM/TM

07.10 S R1=GR;S R2=XR;S I1=GI;S I2=XI;D 5
07.20 S GR=RE+CR(I);S GI=IM+CI(I);R

08.10 S R1=DR;S R2=XR;S I1=DI;S I2=XI;D 5
08.20 S DR=RE+(N-I)*CR(I);S DI=IM+(N-I)*CI(I);R

09.10 F I=1,N+1;D 9.5;R
09.50 S R1=YR;S R2=CR(I-1);S I1=YI;S I2=CI(I-1);D 5;D 9.6
09.60 S CR(I)=CR(I)+RE;S CI(I)=CI(I)+IM

10.10 S RN=0;S ER=0
10.20 S R1=-CR(1);S R2=CR(0);S I1=-CI(1);S I2=CI(0); D 6
10.30 S YR=RE;S YI=IM;D 1.9; Q
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@x|