C-FOCAL,1969

11.05 T !! "ROOTS OF POLYNOMIAL - PART II
11.10 T !! "POLYNOMIAL DIVISION FOR REAL ROOT REDUCTION
11.20 A ! "ORDER OF POLYNOMIAL  ",N
11.30 S M=1;S B(1)=1
11.40 T !! "TYPE IN POLYNOMIAL COEFFS",!
11.45 F I=1,N+1;A A(I)
11.50 A ! "TYPE IN VALUE OF REAL ROOT  ",RR
11.55 S B(2)=-RR
11.60 T !! "POLYNOMIAL COEFFS",!
11.65 F I=1,N+1;T %,A(I), !
11.70 T !! "VALUE OF REAL ROOT  ",RR,!

12.10 S LL=N+1
12.20 F I=1,LL;S G(I)=A(I)
12.30 S LL=N-M+1
12.40 F I=1,LL;D 13
12.50 F I=1,M;S R(I)=G(N-M+I+1)
12.60 S L=M-1
12.70 IF (R(1)) 12.9,12.75,12.9
12.75 IF (L) 12.9,12.9,12.77
12.77 F I=1,L;S R(I)=R(I+1)
12.80 S L=L-1;G 12.7
12.90 G 15.1

13.10 S Q(I)=G(I)/B(1)
13.20 F J=1,M;D 14
13.30 R

14.10 S W=Q(I)*B(J+1)
14.20 S G(I+J)=G(I+J)-W
14.30 IF (FABS(G(I+J))-.0001*FABS(W)) 14.4,14.4,14.5
14.40 S G(I+J)=0
14.50 R

15.10 T !! "RESULTS OF DIVISION
15.20 T !! "REDUCED POLYNOMIAL",!
15.25 F I=1,LL;T %,Q(I), !
15.40 Q
*