! ニュートン法によって解を求める。
DEF f(x)=\(x^{2}\)*(x-a)-b
DEF g(x)=2*x*(x-a)+\(x^{2}\) ! g(x)=f'(x)
INPUT a,b
input c
LET x=c
LET E=EPS(1)*10
DO
LET x1=x-f(x)/g(x)
PRINT x1
IF ABS((x1-x)/x)＜E THEN EXIT DO
LET x=x1
LOOP
END

ａ＝５３１、ｂ＝９．９９×１０＾６
ｃ＝１０００を入力すると、
５６２．５６５９３２９５２４２４が求まる。
したがって、
近似解５６３となる。

Solve \(x^{2}\)(x-531)=9.99*1\(0^{6}\) by Newton's method.
Put a=9.99*1\(0^{6}\), f(x)=\(x^{2}\)(x-531)-a.
The root of f(x)=0 is bigger than 531.
Then x0=531.
x1=x0-f(x0)/f'(x0)=531-(-a)/53\(1^{2}\)=531+35.4=566.4
x2=x1-f(x1)/f'(x1)=...