質問<1063>
「「3次方程式」」
日付 2003/1/12
質問者 クロ


どなたか教えてください。

T\(2^{2}\)(T2-454)=2.97×1\(0^{8}\)
T2=858
の答えのとき方をおしえてください。
よろしくお願いします。

お返事(武田)
日付 2003/1/13
回答者 武田


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

a=454、b=2.97×10^8
c=1000を入力すると、
857.712870633521が求まる。
したがって、
近似解858となる。