逆数は、ニュートン法で求められる。 Aを求めたい数の逆数とすると(= 1/Aを求めたいとする)
f(x)=1/x-A=0
f’(x)=-1/x^2
したがって、逆数のニュートン法の漸化式は
Xn+1 :=Xn-(f(Xn)/f’(Xn))
=Xn-(1/Xn-A)/(-1/Xn^2)
=Xn+(Xn-A・Xn^2)
=2・Xn-A・Xn^2
OPTION ARITHMETIC DECIMAL_HIGH
LET A=3 ! 1/3の場合
LET x=0.4 ! 初期値
FOR J=1 TO 5
LET x=2*x-A*x^2
NEXT J
PRINT "1/3=";x
END
操作を5回施した結果: 1/3= .33333333333333333333331901677568