ハリー法

 投稿者:しばっち  投稿日:2013年 5月19日(日)22時39分44秒
  !' ハリー法(HALLEY法)
PUBLIC NUMERIC T
INPUT  PROMPT "X=":T
LET EPS=1E-8
LET X=T
DO
   LET X=X-H(X)/DH(X)
   PRINT X
LOOP UNTIL ABS(F(X))<EPS
PRINT X;X^3
END

EXTERNAL  FUNCTION H(X)
LET H=F(X)/SQR(DF(X))
END FUNCTION

EXTERNAL  FUNCTION DH(X)
LET DH=SQR(DF(X))-F(X)*DF2(X)/(2*DF(X)*SQR(DF(X)))
END FUNCTION

EXTERNAL  FUNCTION F(X)!'   f(x)
LET  F=X^3-T
END FUNCTION

EXTERNAL  FUNCTION DF(X)!'  f'(x)
LET  DF=3*X^2
END FUNCTION

EXTERNAL  FUNCTION DF2(X)!'  f''(x)
LET DF2=6*X
END FUNCTION
 

戻る