本帖最后由 kyra 于 2024-5-16 10:28 编辑
朋友们,我想解方程x³=8的方程,用牛顿法,结果x=4.5,是哪里的步骤有错误呢
[Fortran] 纯文本查看 复制代码 program main
X=1.
K=0
100 K=K+1
FX=X*X*X-8
C 这儿输入公式
DFX=2*X*X
CALL NEWTON (X,FX,DFX,TOL)
IF (TOL.LT.1.E-4) GOTO 200
GOTO 100
200 WRITE(*,*)X
STOP
END program main
SUBROUTINE NEWTON(X,FX,DFX,TOL)
DX=-FX/DFX
TOLER=ABS(DX)
X=X+DX
RETURN
END
|