|
投稿者：しばっち
投稿日：2013年 5月19日(日)22時30分29秒
|
|
|
|
!' ステッフェッセン法(STEFFENSEN法)
PUBLIC NUMERIC T
INPUT PROMPT "X=":T
LET EPS=1E-8
LET X=T
DO
LET X=X-F(X)/G(X)
PRINT X
LOOP UNTIL ABS(F(X))<EPS
PRINT X;X*X
END
EXTERNAL FUNCTION F(X)
LET F=X*X-T
END FUNCTION
EXTERNAL FUNCTION G(X) !' f'(x)の代用
LET G=(F(F(X)+X)-F(X))/F(X)
END FUNCTION
|
|
|