30 - PARAMETER OPTIMIZATION OF DYNAMIC SYSTEMS
----------------------------------------------


      CALL EIUD30(N,NE,X,YA,EW,EY,FMIN,XMAX,TA,TAMAX,NEXT,IEXT,IERR)

         N - NUMBER OF VARIABLES (OUTPUT)
         NE - NUMBER OF DIFFERENTIAL EQUATIONS (OUTPUT)
         X(N) - INITIAL VECTOR OF VARIABLES (OUTPUT)
         YA(NE) - INITIAL SOLUTION VECTOR (OUTPUT)
         EW(NE) - WEIGHT COEFFICIENTS IN THE TERMINAL FUNCTION (OUTPUT)
         EY(NE) - REQUIRED VALUES IN THE TERMINAL FUNCTION (OUTPUT)
         FMIN - LOWER BOUND OF OBJECTIVE FUNCTION (OUTPUT)
         XMAX - UPPER BOUND OF STEPSIZE (OUTPUT)
         TA - INITIAL TIME (OUTPUT)
         TAMAX - FINAL TIME (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT): THERE IS 4 PROBLEMS
                (NEXT FROM 1 TO 4)
         IEXT - TYPE OF EXTREMUM (0-MINIMUM, 1-MAXIMUM)
         IERR - CONFLICT IN INPUT DATA IF IERR <> 0 (OUTPUT)


      CALL EEFU30(N,NE,KE,TA,YA,X,FE,NEXT)

         N - NUMBER OF VARIABLES (INPUT)
         NE - NUMBER OF DIFFERENTIAL EQUATIONS (INPUT)
         KE - INDEX OF DIFFERENTIAL EQUATION (INPUT)
         TA - VALUE OF THE TIME (INPUT)
         YA(NE) - SOLUTION VECTOR (INPUT)
         X(N) - ACTUAL VECTOR OF VARIABLES (INPUT)
         FE - VALUE OF THE RIGHT HAND SIDE FUNCTION (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT)


      CALL EEGU30(N,NE,KE,TA,YA,X,DE,GE,NEXT)

         N - NUMBER OF VARIABLES (INPUT)
         NE - NUMBER OF DIFFERENTIAL EQUATIONS (INPUT)
         KE - INDEX OF DIFFERENTIAL EQUATION (INPUT)
         TA - VALUE OF THE TIME (INPUT)
         YA(NE) - INITIAL CONDITIONS VECTOR (INPUT)
         X(N) - ACTUAL VECTOR OF VARIABLES (INPUT)
         DE(N) - DERIVATIVES OF THE RIGHT HAND SIDE FUNCTION (OUTPUT)
         GE(N) - GRADIENT OF THE RIGHT HAND SIDE FUNCTION (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT)


      CALL EYFU30(N,KE,X,FE,NEXT)

         N - NUMBER OF VARIABLES (INPUT)
         KE - INDEX OF DIFFERENTIAL EQUATION (INPUT)
         X(N) - ACTUAL VECTOR OF VARIABLES (INPUT)
         FE - VALUE OF THE INITIAL CONDITION (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT)


      CALL EYGU30(N,KE,X,GE,NEXT)

         N - NUMBER OF VARIABLES (INPUT)
         KE - INDEX OF DIFFERENTIAL EQUATION (INPUT)
         X(N) - ACTUAL VECTOR OF VARIABLES (INPUT)
         GE(N) - GRADIENT OF THE INITIAL CONDITION (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT)


RECOMMENDED INPUT VALUE N=5, NE=5.
