C-FOCAL,1969 01.10 T !! "REAL ROOTS OF POLYNOMIAL - PART 1 01.30 A ! "ORDER OF POLYNOMIAL ",N 01.40 T !! "TYPE POLYNOMIAL COEFFS",! 01.45 F I=1,N+1;A AM(I) 01.50 D 1.1 01.55 T !! "POLYNOMIAL COEFFS",! 01.60 F I=1,N+1;T %,AM(I),! 01.70 S ST=0 01.80 F I=1,N+1;S ST=ST+FABS(AM(I)) 02.05 S XS=-ST/2;S T=0 02.10 S X2=XS;S X1=0 02.15 D 3 02.20 I (Y1*Y2)4.1,3.25,2.25 02.25 S X2=X2+XS;S X1=X1+XS 02.30 I (FABS(X2)-ST)2.15,2.15,2.4 02.40 I (X2)2.5,2.5,2.55 02.50 S XS=-XS;G 2.1 02.55 S T=T+1;I (T-3)2.57,2.6,2.6 02.57 S XS=-XS/10;G 2.1 02.60 T !! "FUNCTION VALUE HAS NOT CROSSED REAL AXIS DURING 02.65 T ! "PRELIMINARY SEARCH - ROOTS MAY BE COMPLEX",!;Q 03.10 S X=X1;D 7 03.15 S Y1=Z;S X=X2;D 7 03.20 S Y2=Z;R 03.25 T !! "ONE ROOT IS " 03.30 IF (Y2)3.5,3.4,3.5 03.40 T X2;S XR=X2;G 6.1 03.50 T X1;S XR=X1;G 6.1 04.10 S X=X2;D 7 04.12 S Y=Z 04.13 S X=X1 04.15 D 7 04.20 IF (FABS(X2-X)-.1)5.1,5.1,4.25 04.25 IF (Y*Z)4.3,4.6,4.4 04.30 S X1=X;G 4.5 04.40 S X2=X 04.50 S X=(X1+X2)/2;G 4.15 04.60 D 3.25;T X;S XR=X;G 6.1 05.10 S XC=X 05.15 F I=1,N;S AP(I)=AM(I)*(N-I+1) 05.20 S X=XC;D 7 05.25 S YA=Z 05.30 S YB=AP(1);F I=1,N-1;S YB=YB*X+AP(I+1) 05.40 S XN=XC-YA/YB 05.50 I (FABS(XN-XC)-.00001)5.7,5.7,5.6 05.60 S XC=XN;G 5.2 05.70 D 3.25;T XN;S XR=XN;Q 07.10 S Z=AM(1) 07.20 F I=1,N;S Z=Z*X+AM(I+1) 07.30 R *o