5 ! Find integer coefficient less than M for a polynome of degree 4
6 ! with a root near K (use UTIL/1 binary)
10 T=TIME  @ CLEAR @ K=PI  @ M=20 @ E=INF  @ N=M+1 @ O=N*K @ P=N/K @ Q=(P+M)/K
20 T0=K @ FOR A=1 TO M @ AWRITE 0,0 @ DISP A,T0 @ IF ABS (T0)*K*K-N>O*(K+1) THEN 90
25 B0=MAX (-M,INT (-Q-T0+.5)) @ B1=MIN (M,INT (Q-T0+.5)) @ IF B0>B1 THEN 90
27    ! AWRITE 0,40 @ DISP B0,B1
30 T1=(T0+B0)*K @ FOR B=B0 TO B1 @ AWRITE 1,0 @ DISP B,T1 @ IF ABS (T1)*K-N>O THEN 80
40 C0=MAX (-M,INT (-P-T1+.5)) @ C1=MIN (M,INT (P-T1+.5)) @ IF C0>C1 THEN 80
42    ! AWRITE 1,40 @ DISP C0,C1
50 T2=(T1+C0)*K @ FOR C=C0 TO C1 @ D=-INT (T2+.5) @ IF ABS (D)>M THEN 70
60 IF ABS (T2+D)<E THEN F=A @ G=B @ H=C @ I=D @ E=ABS (T2+D) @ AWRITE 5,0 @ DISP A;B;C;D;E
70 T2=T2+K @ NEXT C
80 T1=T1+K @ NEXT B
90 T0=T0+K @ NEXT A
100 AWRITE 6,0 @ DISP TIME -T @ END
