3体問題

日常の数学・物理    ーー それでも地球はまわっている



ニュートン法-逆数計算


逆数は、ニュートン法で求められる。 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

十進BACIC

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