| 本帖最后由 kyra 于 2024-5-16 10:28 编辑 
 朋友们,我想解方程x³=8的方程,用牛顿法,结果x=4.5,是哪里的步骤有错误呢
 
 [Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode 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
 
 
 |