|
|
!' 減速ニュートン法
LET X=3 !'初期値
LET EPS=1E-8
DO
LET U=2
DO
LET U=U/2
LET Y=X-U*F(X)/DF(X) !'減速ニュートン法
LOOP UNTIL ABS(F(Y))<(1-U/2)*ABS(F(X))
LET X=Y
LOOP UNTIL ABS(F(X))<EPS
PRINT X;F(X)
LET X=3 !'初期値
LET I=0
DO
LET X=X-F(X)/DF(X) !'ニュートン法
LET I=I+1
IF I>100 THEN
PRINT"収束しません"
STOP
END IF
LOOP UNTIL ABS(F(X))<EPS
PRINT X;F(X)
END
EXTERNAL FUNCTION F(X)
LET F=3*ATN(X-1)+X/4
END FUNCTION
EXTERNAL FUNCTION DF(X)
LET DF=3/((X-1)^2+1)+1/4
END FUNCTION
|
|