EXTERNAL FS
X0=1.5
EPS=1.0E-06
CALL DNEWT(X0,EPS,FS,L)
IF (L.NE.0) THEN
WRITE(*,10) X0
END IF
10FORMAT(1X,'X=',E15.6)
END
SUBROUTINE FS(X,F,DY)
F=X*X*(X-1.0)-1.0
DY=3.0*X*X-2.0*X
RETURN
END
SUBROUTINE DNEWT(X,EPS,FS,L)
L=60
CALL FS(X,F,DY)
10IF (ABS(DY)+1.0.EQ.1.0) THEN
L=0
WRITE(*,20)
RETURN
END IF
20FORMAT(1X,' ERR')
X1=X-F/DY
CALL FS(X1,F,DY)
IF ((ABS(X1-X).GE.EPS).OR.(ABS(F).GE.EPS)) THEN
L=L-1
X=X1
IF (L.EQ.0) RETURN
GOTO 10
END IF
X=X1
RETURN
END
----------------------------------------------------------------------------------------------------------------------------------牛顿法---------哥负责的帮你调试好了.题中有个例子直接可以解答!没有一点问题.不过最好用Fortran77..90的没事实验.有什么问题可以问.请采纳哦这么多时间哥的
