您所在位置：主题：关于微分方程组 IVPRK
xwzj
[专家分：0] 发布于 2012-08-17 18:09:00
PROGRAM MAIN
USE IMSL
IMPLICIT NONE
INTEGER,PARAMETER :: N=2
INTEGER IDO , ISTEP,J
REAL(KIND=8) PARAM(50),T,TEND,Y(N),DT
REAL(KIND=8),PARAMETER :: TOL=0.001
EXTERNAL FCN
PARAM=0
Y(1)=0.0
Y(2)=1.0
T=0
DT=0.1
102 FORMAT (1X,'T=',F20.10,1X,'Y=',F20.10,1X,'Z=',F20.10)
END DO
USE IMSL
IMPLICIT NONE
INTEGER,PARAMETER :: N=2
INTEGER IDO , ISTEP,J
REAL(KIND=8) PARAM(50),T,TEND,Y(N),DT
REAL(KIND=8),PARAMETER :: TOL=0.001
EXTERNAL FCN
PARAM=0
Y(1)=0.0
Y(2)=1.0
T=0
DT=0.1
IDO=1
TEND= T+DT
CALL DIVPRK(IDO,N,FCN,T,TEND,TOL,PARAM,Y)
102 FORMAT (1X,'T=',F20.10,1X,'Y=',F20.10,1X,'Z=',F20.10)
END DO
END PROGRAM
subroutine FCN(N,T,Y,YPRIME)
implicit none
integer N
REAL(KIND=8) T,Y,YPRIME(N)
YPRIME(1)=3*Y(1)-2*Y(2)
YPRIME(2)=2*Y(1)-Y(2)
RETURN
END SUBROUTINE
为什么总是错误，倒地哪里错了啊，求问高手! 感激不尽啊