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

B. OC-ODE
IMPLICIT NONE
INTEGER IUSER(*)
DOUBLEPRECISION T,X(*),USER(*)
X(1) = 0.25D0
X(2) = 0.8D0
X(3) = 1.0D0
X(4) = 0.0D0
RETURN
END
c-------------------------------------------------
c Initial Estimate CONTROL VARIABLES
c-------------------------------------------------
SUBROUTINE INESTU( T, U, IBOOR, IUSER, USER )
IMPLICIT NONE
INTEGER IBOOR,IUSER(*)
DOUBLEPRECISION T,U(*),USER(*)
U(1)= 0.0D0
U(2)= 0.0D0
U(3)= 0.0D0
RETURN
END
c-------------------------------------------------
c Mass Matrix
c-------------------------------------------------
SUBROUTINE MASS( NX,T,X,XP,U,P,MMASS,IUSER,USER )
IMPLICIT NONE
INTEGER NX,IUSER(*)
DOUBLE PRECISION T,X(*),XP(*),U(*),P(*),
+ MMASS(NX,*),USER(*)
MMASS(1,1)=1.0d0
MMASS(2,1)=1.0d0
MMASS(3,1)=1.0d0
MMASS(4,1)=1.0d0
RETURN
END
c-------------------------------------------------
c Iteration Matrix
c-------------------------------------------------
SUBROUTINE ITMAT( T,X,XP,U,P,JAC1,JAC2,CJ,NX,
+ NSTAB,IUSER,USER )
IMPLICIT NONE
INTEGER NX,NSTAB,IUSER(*)
DOUBLE PRECISION T,X(*),XP(*),U(*),P(*),JAC1(NX,*),
142