*  NLP written by GAMS Convert at 03/28/07 20:24:11
*  
*  Equation counts
*      Total        E        G        L        N        X        C
*         13       13        0        0        0        0        0
*  
*  Variable counts
*                   x        b        i      s1s      s2s       sc       si
*      Total     cont   binary  integer     sos1     sos2    scont     sint
*         29       29        0        0        0        0        0        0
*  FX      0        0        0        0        0        0        0        0
*  
*  Nonzero counts
*      Total    const       NL      DLL
*         71       59       12        0
*
*  Solve m using NLP maximizing objvar;


Variables  x1,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,objvar;

Positive Variables  x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12;

Equations  e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13;


e1..    0.005*x1 + 0.00714285714285714*x5 - x9 =E= 0;

e2..    0.005*x2 + 0.00714285714285714*x6 - x10 =E= 0;

e3..    0.005*x3 + 0.00714285714285714*x7 - x11 =E= 0;

e4..    0.005*x4 + 0.00714285714285714*x8 - x12 =E= 0;

e5..    x2 - x14 + x21 - x22 =E= 0;

e6..    x3 - x15 + x22 - x23 =E= 0;

e7..    x4 - x16 + x23 - x24 =E= 0;

e8..    x6 - x18 + x25 - x26 =E= 0;

e9..    x7 - x19 + x26 - x27 =E= 0;

e10..    x8 - x20 + x27 - x28 =E= 0;

e11..    x1 - x13 - x21 =E= -10;

e12..    x5 - x17 - x25 =E= 0;

e13..  - (-2.5*(sqrt(x21*x21) - x21) - 1.25*(sqrt(x21*x21) + x21) + (-2.5*(
      sqrt(x22*x22) - x22)) - 1.25*(sqrt(x22*x22) + x22) + (-2.5*(sqrt(x23*x23)
       - x23)) - 1.25*(sqrt(x23*x23) + x23) + (-2.5*(sqrt(x24*x24) - x24)) - 
      1.25*(sqrt(x24*x24) + x24) + (-3*(sqrt(x25*x25) - x25)) - 1.5*(sqrt(x25*
      x25) + x25) + (-3*(sqrt(x26*x26) - x26)) - 1.5*(sqrt(x26*x26) + x26) + (-
      3*(sqrt(x27*x27) - x27)) - 1.5*(sqrt(x27*x27) + x27) + (-3*(sqrt(x28*x28)
       - x28)) - 1.5*(sqrt(x28*x28) + x28) - (1550*(sqrt((-35 + x9)*(-35 + x9))
       + x9) + 1500*(sqrt((-35 + x10)*(-35 + x10)) + x10) + 1850*(sqrt((-30 + 
      x11)*(-30 + x11)) + x11) + 1550*(sqrt((-35 + x12)*(-35 + x12)) + x12)))
       + 10*x1 + 10*x2 + 10*x3 + 10*x4 + 11*x5 + 11*x6 + 11*x7 + 11*x8 - 25*x13
       - 26*x14 - 27*x15 - 27*x16 - 30*x17 - 35*x18 - 37*x19 - 39*x20 + objvar
       =E= 216500;

* set non default bounds

x9.up = 42; 
x10.up = 42; 
x11.up = 40; 
x12.up = 42; 
x13.lo = 3000; x13.up = 6000; 
x14.lo = 3000; x14.up = 6000; 
x15.lo = 3000; x15.up = 4000; 
x16.lo = 3000; x16.up = 6500; 
x17.lo = 2000; x17.up = 4000; 
x18.lo = 2000; x18.up = 2500; 
x19.lo = 2000; x19.up = 3500; 
x20.lo = 2000; x20.up = 4200; 

* set non default levels

x13.l = 3000; 
x14.l = 3000; 
x15.l = 3000; 
x16.l = 3000; 
x17.l = 2000; 
x18.l = 2000; 
x19.l = 2000; 
x20.l = 2000; 

* set non default marginals


Model m / all /;

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

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

Solve m using NLP maximizing objvar;