disser1.pdf (2006 KB) - Ernst-Moritz-Arndt-Universität Greifswald

B. OC-ODE
c Call to OC-ODE
c OPEN( INFO(9),FILE=’OUT’,STATUS=’UNKNOWN’)
CALL OCODE( T, XL, XU, UL, UU, P, G, BC,
+ TOL, TAUU, TAUX, LIW, LRW, IRES,
+ IREALTIME, NREALTIME, HREALTIME,
+ IADJOINT, RWADJ, LRWADJ, IWADJ, LIWADJ, .FALSE.,
+ MERIT,IUPDATE,LENACTIVE,ACTIVE,IPARAM,PARAM,
+ DIM,INFO,IWORK,RWORK,SOL,NVAR,IUSER,USER)
PRINT*,’Ausgabe der Loesung: NVAR=’,NVAR
WRITE(*,’(E30.16)’) (SOL(I),I=1,NVAR)
c CLOSE(INFO(9))
c READ(*,*)
END
c-------------------------------------------------
c Objective Function
c-------------------------------------------------
SUBROUTINE OBJ( X0, XF, TF, P, V, IUSER, USER )
IMPLICIT NONE
INTEGER IUSER(*)
DOUBLEPRECISION X0(*),XF(*),TF,P(*),V,USER(*)
V = -XF(4)
RETURN
END
c-------------------------------------------------
c Differential Equation
c-------------------------------------------------
SUBROUTINE DAE( T, X, XP, U, P, F, IFLAG, IUSER, USER )
IMPLICIT NONE
INTEGER IFLAG,IUSER(*)
DOUBLEPRECISION T,X(*),XP(*),U(*),P(*),F(*),USER(*)
F(1) = 0.4D0*X(1)*(1-X(1))-U(1)*1.5D0*25.0D0*X(1)/1.0
+ D+5-0.02D0*X(1)*X(2)/1.2D0-0.02D0*X(1)*X(3)/1.3D0
F(2) = 0.6D0*X(2)*(1-X(2)/1.2D0)-U(2)*6.4D0*250.0D0*X(2)
+ /1.2D+6-0.0125D0*X(1)*X(2)/1.2D0 - 0.01D0*X(2)*X(3)
+ /1.56D0
F(3) = 0.6D0*X(3)*(1-X(3)/1.3D0)-U(3)*6.4D0*250.0D0*X(3)
+ /1.3D+6-0.0125D0*X(1)*X(3)/1.3D0 - 0.01D0*X(2)*X(3)
+ /1.56D0
F(4) = ((1130.0D0*U(1)*1.5D0*250.0D0*X(1)/1.0D+6)+(270.0D0
140