Halley's Method

 投稿者:しばっち  投稿日:2015年10月28日(水)22時04分34秒
 
!'Halley's Method
INPUT T
LET X=T
DO
   LET X=X-2*F(X)*DF1(X)/(2*DF1(X)^2-F(X)*DF2(X))
   PRINT X
LOOP UNTIL ABS(F(X))<1E-13
PRINT X;X^3

FUNCTION F(X)
   LET F=X*X*X-T
END FUNCTION

FUNCTION DF1(X)
   LET DF1=3*X^2
END FUNCTION

FUNCTION DF2(X)
   LET DF2=6*X
END FUNCTION
END
 

戻る