投稿者:しばっち
投稿日: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
|
|
|