Catalyst Mixing COPS 2.0

   {   Catmix.mpl   }

   {   GAMS Model Library, http://www.gams.com/modlib/libhtml/catmix.htm }

   {   Catalyst Mixing COPS 2.0,  NLP,  Size: 100x153  }



TITLE
    Catmix;

Option
    Modeltype=nonlinear
    Parsertype=extended

DATA
    Dnh := 50;

INDEX
    nh := 0..Dnh;
    i := nh;

DATA
    Tf := 1;
    X1_0 := 1;
    X2_0 := 0;
    Alpha := 0;
    H := (Tf/Dnh);

VARIABLES
    U[nh];
    X1[nh];
    X2[nh];

MODEL
    MIN Obj = -1 + X1[nh:=Dnh] + X2[nh:=Dnh]
              + Alpha * H * SUM(nh: (U[nh+1] - U[nh-1])^2);

SUBJECT TO

    Ode1[nh>0]: X1
              =
                X1[nh-1] + (H/2)*(U[nh-1]*(10*X2[nh-1]-X1[nh-1])+U*(10*X2-X1));

    Ode2[nh>0]: X2
              =
                X2[nh-1] + (H/2)*(U[nh-1]*(U[nh-1]*(X1[nh-1]-10*X2[nh-1])
              - (1-U[nh-1])*X2[nh-1]+U*(X1-10*X2)
              - (1-U)*X2));

BOUNDS
   0 <= U <= 1;
   -100 <= X1 <= 100;
   -100 <= X2 <= 100;
   X1[nh=0] = X1_0;
   X2[nh=0] = X2_0;

END