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

B.1. Programm 1. Berechnung mit stückweise konstanten Steuerungen
+ *U(2)*6.4D0*250.0D0*X(2)/1.2D+6)+(460.0D0*U(3)*
+ 6.4D0*250.0D0*X(3)/1.3D+6)-500.0D0*250.0D0*(U(1)+U(2)
+ +U(3))/1.0D+6)*exp(-0.06D0*T)
RETURN
END
c-------------------------------------------------
c Nonlinear Constraints
c-------------------------------------------------
SUBROUTINE NLCSTR( T, X, U, P, G, IUSER, USER )
IMPLICIT NONE
INTEGER IUSER(*)
DOUBLEPRECISION T,X(*),U(*),P(*),G(*),USER(*)
DOUBLE PRECISION L
L = 1900.0D0
G(1) = -U(1)-U(2)-U(3)+L
RETURN
END
c-------------------------------------------------
c JACOBIAN OF Nonlinear Constraints
c-------------------------------------------------
SUBROUTINE JACNLC( T, X, U, P, NG, GJAC, IUSER, USER )
IMPLICIT NONE
INTEGER NG,IUSER(*)
DOUBLEPRECISION T,X(*),U(*),P(*),GJAC(NG,*),USER(*)
RETURN
END
c-------------------------------------------------
c Boundary Conditions
c-------------------------------------------------
SUBROUTINE BDCOND( T0,TF,X0,XF,U0,UF,P,PSI,IUSER,USER )
IMPLICIT NONE
INTEGER IUSER(*)
DOUBLEPRECISION T0,TF,X0(*),XF(*),U0(*),UF(*),P(*),
+ PSI(*),USER(*)
PSI(1) = XF(1)
PSI(2) = XF(2)
PSI(3) = XF(3)
RETURN
END
c-------------------------------------------------
c Initial Estimate STATE VARIABLES
c-------------------------------------------------
SUBROUTINE INESTX( T, X, IUSER, USER )
141