*  QCP written by GAMS Convert at 01/05/04 11:35:41
*  
*  Equation counts
*     Total       E       G       L       N       X       C
*        54      45       0       9       0       0       0
*  
*  Variable counts
*                 x       b       i     s1s     s2s      sc      si
*     Total    cont  binary integer    sos1    sos2   scont    sint
*        47      47       0       0       0       0       0       0
*  FX     0       0       0       0       0       0       0       0
*  
*  Nonzero counts
*     Total   const      NL     DLL
*       195     177      18       0
*
*  Solve m using QCP minimizing objvar;


Variables  objvar,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18
          ,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35
          ,x36,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47;

Equations  e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19
          ,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36
          ,e37,e38,e39,e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,e50,e51,e52,e53
          ,e54;
Scalar MU / 0 /;


e1..    objvar - x29 =E= 0;

e2.. 0.5*x30*x30 + 0.5*x31*x31 =L= 0.5;

e3.. 0.5*x32*x32 + 0.5*x33*x33 =L= 0.5;

e4.. 0.5*x34*x34 + 0.5*x35*x35 =L= 0.5;

e5.. 0.5*x36*x36 + 0.5*x37*x37 =L= 0.5;

e6.. 0.5*x38*x38 + 0.5*x39*x39 =L= 0.5;

e7.. 0.5*x40*x40 + 0.5*x41*x41 =L= 0.5;

e8.. 0.5*x42*x42 + 0.5*x43*x43 =L= 0.5;

e9.. 0.5*x44*x44 + 0.5*x45*x45 =L= 0.5;

e10.. 0.5*x46*x46 + 0.5*x47*x47 =L= 0.5;

e11..    0.166666666666667*x3 - 0.166666666666667*x4 - 0.166666666666667*x6
       - 0.166666666666667*x7 =E= 0;

e12..  - 0.166666666666667*x5 + 0.166666666666667*x7 - 0.166666666666667*x8
       - 0.166666666666667*x10 =E= 0;

e13..    0.166666666666667*x6 - 0.166666666666667*x7 - 0.166666666666667*x9
       - 0.166666666666667*x10 =E= 0;

e14..  - 0.166666666666667*x8 + 0.166666666666667*x10 =E= 0;

e15..    0.166666666666667*x9 - 0.166666666666667*x10 =E= 0;

e16..    0.166666666666667*x2 - 0.166666666666667*x4 - 0.166666666666667*x11
       - 0.166666666666667*x13 =E= 0;

e17..    0.166666666666667*x2 + 0.166666666666667*x4 + 0.166666666666667*x5
       - 0.166666666666667*x7 - 0.166666666666667*x11 + 0.166666666666667*x13
       - 0.166666666666667*x14 - 0.166666666666667*x16 =E= 0;

e18..    0.166666666666667*x3 + 0.166666666666667*x4 - 0.166666666666667*x6
       + 0.166666666666667*x7 + 0.166666666666667*x12 - 0.166666666666667*x13
       - 0.166666666666667*x15 - 0.166666666666667*x16 =E= 0;

e19..    0.166666666666667*x5 + 0.166666666666667*x7 + 0.166666666666667*x8
       - 0.166666666666667*x10 - 0.166666666666667*x14 + 0.166666666666667*x16
       - 0.166666666666667*x17 - 0.166666666666667*x19 =E= 0;

e20..    0.166666666666667*x6 + 0.166666666666667*x7 - 0.166666666666667*x9
       + 0.166666666666667*x10 + 0.166666666666667*x15 - 0.166666666666667*x16
       - 0.166666666666667*x18 - 0.166666666666667*x19 =E= 0;

e21..    0.166666666666667*x8 + 0.166666666666667*x10 - 0.166666666666667*x17
       + 0.166666666666667*x19 =E= 0;

e22..    0.166666666666667*x9 + 0.166666666666667*x10 + 0.166666666666667*x18
       - 0.166666666666667*x19 =E= 0;

e23..    0.166666666666667*x11 - 0.166666666666667*x13 - 0.166666666666667*x20
       - 0.166666666666667*x22 =E= 0;

e24..    0.166666666666667*x11 + 0.166666666666667*x13 + 0.166666666666667*x14
       - 0.166666666666667*x16 - 0.166666666666667*x20 + 0.166666666666667*x22
       - 0.166666666666667*x23 - 0.166666666666667*x25 =E= 0;

e25..    0.166666666666667*x12 + 0.166666666666667*x13 - 0.166666666666667*x15
       + 0.166666666666667*x16 + 0.166666666666667*x21 - 0.166666666666667*x22
       - 0.166666666666667*x24 - 0.166666666666667*x25 =E= 0;

e26..    0.166666666666667*x14 + 0.166666666666667*x16 + 0.166666666666667*x17
       - 0.166666666666667*x19 - 0.166666666666667*x23 + 0.166666666666667*x25
       - 0.166666666666667*x26 - 0.166666666666667*x28 =E= 0;

e27..    0.166666666666667*x15 + 0.166666666666667*x16 - 0.166666666666667*x18
       + 0.166666666666667*x19 + 0.166666666666667*x24 - 0.166666666666667*x25
       - 0.166666666666667*x27 - 0.166666666666667*x28 =E= 0;

e28..    0.166666666666667*x17 + 0.166666666666667*x19 - 0.166666666666667*x26
       + 0.166666666666667*x28 =E= 0;

e29..    0.166666666666667*x18 + 0.166666666666667*x19 + 0.166666666666667*x27
       - 0.166666666666667*x28 =E= 0;

e30..    0.166666666666667*x20 - 0.166666666666667*x22 - 0.166666666666667*x29
       =E= 0;

e31..    0.166666666666667*x20 + 0.166666666666667*x22 + 0.166666666666667*x23
       - 0.166666666666667*x25 - 0.333333333333333*x29 =E= 0;

e32..    0.166666666666667*x21 + 0.166666666666667*x22 - 0.166666666666667*x24
       + 0.166666666666667*x25 =E= 0;

e33..    0.166666666666667*x23 + 0.166666666666667*x25 + 0.166666666666667*x26
       - 0.166666666666667*x28 - 0.333333333333333*x29 =E= 0;

e34..    0.166666666666667*x24 + 0.166666666666667*x25 - 0.166666666666667*x27
       + 0.166666666666667*x28 =E= 0;

e35..    0.166666666666667*x26 + 0.166666666666667*x28 - 0.166666666666667*x29
       =E= 0;

e36..    0.166666666666667*x27 + 0.166666666666667*x28 =E= 0;

e37..    0.5*x2 - 0.5*x3 - x30 =E= 0;

e38..    x4 - x31 =E= 0;

e39..    0.5*x5 - 0.5*x6 - x32 =E= 0;

e40..    x7 - x33 =E= 0;

e41..    0.5*x8 - 0.5*x9 - x34 =E= 0;

e42..    x10 - x35 =E= 0;

e43..    0.5*x11 - 0.5*x12 - x36 =E= 0;

e44..    x13 - x37 =E= 0;

e45..    0.5*x14 - 0.5*x15 - x38 =E= 0;

e46..    x16 - x39 =E= 0;

e47..    0.5*x17 - 0.5*x18 - x40 =E= 0;

e48..    x19 - x41 =E= 0;

e49..    0.5*x20 - 0.5*x21 - x42 =E= 0;

e50..    x22 - x43 =E= 0;

e51..    0.5*x23 - 0.5*x24 - x44 =E= 0;

e52..    x25 - x45 =E= 0;

e53..    0.5*x26 - 0.5*x27 - x46 =E= 0;

e54..    x28 - x47 =E= 0;


* set non default bounds


* set non default levels


* set non default marginals


Model m / all /;

m.limrow=0; m.limcol=0;

$if NOT '%gams.u1%' == '' $include '%gams.u1%'

Solve m using QCP minimizing objvar;