質問＜１０６９＞２００３／１／１３
from=クロ
「三次方程式の解き方について」

Ｔ２＾2(Ｔ２－５３１）＝９．９９×１０＾6

Ｔ２＝５６３

の解き方をおしえてください。

お返事２００３／１／１３
from=武田

! ニュートン法によって解を求める。
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

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

お便り２００３／１／１３
from=juin

Solve x^2(x-531)=9.99*10^6 by Newton's method.
Put a=9.99*10^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)/531^2=531+35.4=566.4
x2=x1-f(x1)/f'(x1)=...