程序如下,连接的时候出现如图片所示的错误,请各位前辈赐教!
program NLEQUAT
real x0,x1,XEBC,FEBC,XR,FR,F0,F1
INTEGER IPARA,ICDIS
ITMAX=50
XEBC=1.0E-6
FEBC=1.0E-8
ICDIS=1
PRINT*,"enter IPARA x0 x1"
read 10,IPARA,x0,x1
10 format(I1,2f10.0)
IF(IPARA.EQ.1) THEN
READ 20,ICDIS,ITMAX,XEBC,FEBC
20 FORMAT(2I2,2F10.0)
ENDIF
F0=FUN(X0)
F1=FUN(X1)
IF ((F0*F1).GE.0.) STOP 'ERROR IN ROOT REGEON(X0,X1)'
CALL DHRT
WRITE(*,70)ITER,XR,FR
70 FORMAT(//"NUMBER OF ITERATION:",I15/"EQUATION ROOT",E15.4/"FUNCTION VALUE F(x)=",E15.4)
END
SUBROUTINE DHRT
ITER=0
20 ITER=ITER+1
XR=(X0+X1)/2
FR=FUN(XR)
IF(ICDIS.EQ.1) THEN
WRITE(*,30)ITER,X0,X1,XR,FR
30 FORMAT(1X,I6,3F12.6,F15.6)
ENDIF
IF(ABS(FR).LT.FEBC)RETURN
IF(FR*F0.LT.0.)THEN
X1=XR
F1=FR
ELSE
X0=XR
F0=FR
ENDIF
IF (ITER.GT.ITMAX)RETURN
IF(ABC(X0-X1).GT.XEBC) GOTO 20
RETURN
END
FUNCTION FUN(X)
REAL FUN
FUN=EXP(-X/4.)*(2.-X)-1.0
RETURN
END
![毕业论文遇到棘手问题,求前辈们指点!]()
构建结果.jpg |