diet2.gms
* LP written by GAMS Convert at 03/28/07 20:22:51
*
* Equation counts
* Total E G L N X C
* 13 1 0 12 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 9 9 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 103 103 0 0
*
* Solve m using LP minimizing objvar;
Variables x1,x2,x3,x4,x5,x6,x7,x8,objvar;
Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13;
e1.. - 60*x1 - 8*x2 - 8*x3 - 40*x4 - 15*x5 - 70*x6 - 25*x7 - 60*x8 =L= -700;
e2.. - 10*x1 - 20*x2 - 15*x3 - 35*x4 - 15*x5 - 15*x6 - 25*x7 - 15*x8 =L= -700;
e3.. - 15*x1 - 20*x2 - 10*x3 - 10*x4 - 15*x5 - 15*x6 - 15*x7 - 10*x8 =L= -700;
e4.. - 20*x1 - 10*x3 - 40*x4 - 35*x5 - 30*x6 - 50*x7 - 20*x8 =L= -700;
e5.. - 938*x1 - 2180*x2 - 945*x3 - 278*x4 - 1182*x5 - 896*x6 - 1329*x7
- 1397*x8 =L= 0;
e6.. - 295*x1 - 770*x2 - 440*x3 - 430*x4 - 315*x5 - 400*x6 - 370*x7 - 450*x8
=L= -16000;
e7.. 60*x1 + 8*x2 + 8*x3 + 40*x4 + 15*x5 + 70*x6 + 25*x7 + 60*x8 =L= 20000;
e8.. 10*x1 + 20*x2 + 15*x3 + 35*x4 + 15*x5 + 15*x6 + 25*x7 + 15*x8
=L= 20000;
e9.. 15*x1 + 20*x2 + 10*x3 + 10*x4 + 15*x5 + 15*x6 + 15*x7 + 10*x8
=L= 20000;
e10.. 20*x1 + 10*x3 + 40*x4 + 35*x5 + 30*x6 + 50*x7 + 20*x8 =L= 20000;
e11.. 938*x1 + 2180*x2 + 945*x3 + 278*x4 + 1182*x5 + 896*x6 + 1329*x7
+ 1397*x8 =L= 48200;
e12.. 295*x1 + 770*x2 + 440*x3 + 430*x4 + 315*x5 + 400*x6 + 370*x7 + 450*x8
=L= 24000;
e13.. - 3.19*x1 - 2.59*x2 - 2.29*x3 - 2.89*x4 - 1.89*x5 - 1.99*x6 - 1.99*x7
- 2.49*x8 + objvar =E= 0;
* set non default bounds
x1.lo = 2; x1.up = 10;
x2.lo = 2; x2.up = 10;
x3.lo = 2; x3.up = 10;
x4.lo = 2; x4.up = 10;
x5.lo = 2; x5.up = 10;
x6.lo = 2; x6.up = 10;
x7.lo = 2; x7.up = 10;
x8.lo = 2; x8.up = 10;
* set non default levels
x1.l = 2;
x2.l = 2;
x3.l = 2;
x4.l = 2;
x5.l = 2;
x6.l = 2;
x7.l = 2;
x8.l = 2;
* set non default marginals
Model m / all /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include %gams.u1%
Solve m using LP minimizing objvar;