capacityexpansionexhaustive.gms
* LP written by GAMS Convert at 09/03/08 03:39:34
*
* Equation counts
* Total E G L N X C
* 26 7 9 10 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 49 49 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 133 133 0 0
*
* Solve m using LP minimizing 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,x29,x30,x31,x32,x33,x34,x35,x36
,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48,objvar;
Positive 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,x29,x30,x31,x32,x33,x34
,x35,x36,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48;
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;
e1.. x10 + x11 + x12 + x13 + x14 + x15 =L= 350;
e2.. - x11 - x16 + x17 =E= 0;
e3.. - x13 - x18 + x19 =E= 0;
e4.. - x15 - x20 + x21 =E= 0;
e5.. - x10 + x16 =E= 0;
e6.. - x12 + x18 =E= 0;
e7.. - x14 + x20 =E= 0;
e8.. x1 + x22 + x31 + x40 =G= 25;
e9.. x2 + x23 + x32 + x41 =G= 33;
e10.. x3 + x24 + x33 + x42 =G= 40;
e11.. x4 + x25 + x34 + x43 =G= 55;
e12.. x5 + x26 + x35 + x44 =G= 63;
e13.. x6 + x27 + x36 + x45 =G= 70;
e14.. x7 + x28 + x37 + x46 =G= 85;
e15.. x8 + x29 + x38 + x47 =G= 93;
e16.. x9 + x30 + x39 + x48 =G= 100;
e17.. - x16 + x23 + x26 + x29 =L= 25;
e18.. - x17 + x24 + x27 + x30 =L= 15;
e19.. - x18 + x32 + x35 + x38 =L= 50;
e20.. - x19 + x33 + x36 + x39 =L= 20;
e21.. - x20 + x41 + x44 + x47 =L= 35;
e22.. - x21 + x42 + x45 + x48 =L= 20;
e23.. x22 + x25 + x28 =L= 40;
e24.. x31 + x34 + x37 =L= 80;
e25.. x40 + x43 + x46 =L= 70;
e26.. - 1300*x1 - 1300*x2 - 1300*x3 - 1300*x4 - 1300*x5 - 1300*x6 - 1300*x7
- 1300*x8 - 1300*x9 - 534*x10 - 490*x11 - 464*x12 - 500*x13 - 508*x14
- 507*x15 - 55.1*x22 - 52.2*x23 - 46.4*x24 - 55.1*x25 - 52.2*x26
- 46.4*x27 - 55.1*x28 - 52.2*x29 - 46.4*x30 - 116*x31 - 98.6*x32
- 92.8*x33 - 116*x34 - 98.6*x35 - 92.8*x36 - 116*x37 - 98.6*x38
- 92.8*x39 - 87*x40 - 63.8*x41 - 60.9*x42 - 87*x43 - 63.8*x44 - 60.9*x45
- 87*x46 - 63.8*x47 - 60.9*x48 + objvar =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 LP minimizing objvar;