hydrothermal.gms
<![CDATA[
* NLP written by GAMS Convert at 11/29/04 16:21:11
*
* Equation counts
* Total E G L N X C
* 42 34 8 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 53 53 0 0 0 0 0 0
* FX 5 5 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 127 99 28 0
*
* Solve m using NLP 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,objvar,x49,x50,x51,x52
,x53;
Positive Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x18
,x19,x20,x23,x24,x25,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;
e1.. x17 - x18 - x35 =E= -1.296;
e2.. x18 - x19 - x36 =E= -1.296;
e3.. x19 - x20 - x37 =E= -1.296;
e4.. x20 - x21 - x38 =E= -1.296;
e5.. x22 - x23 + x35 - x39 =E= -0.0864;
e6.. x23 - x24 + x36 - x40 =E= -0.0864;
e7.. x24 - x25 + x37 - x41 =E= -0.0864;
e8.. x25 - x26 + x38 - x42 =E= -0.0864;
e9.. - (7.7167 + 0.2*(0.0957033*x17 - 0.000143410266666667*sqr(x18 - x17) +
0.0957033*x18 - 0.0004302308*x17*x18 + 1.212651e-7*(sqr(x17) + sqr(x18))*(
x17 + x18)))*x35 + x27 =E= 0;
e10.. - (7.7167 + 0.2*(0.0957033*x18 - 0.000143410266666667*sqr(x19 - x18) +
0.0957033*x19 - 0.0004302308*x18*x19 + 1.212651e-7*(sqr(x18) + sqr(x19))*
(x18 + x19)))*x36 + x28 =E= 0;
e11.. - (7.7167 + 0.2*(0.0957033*x19 - 0.000143410266666667*sqr(x20 - x19) +
0.0957033*x20 - 0.0004302308*x19*x20 + 1.212651e-7*(sqr(x19) + sqr(x20))*
(x19 + x20)))*x37 + x29 =E= 0;
e12.. - (7.7167 + 0.2*(0.0957033*x20 - 0.000143410266666667*sqr(x21 - x20) +
0.0957033*x21 - 0.0004302308*x20*x21 + 1.212651e-7*(sqr(x20) + sqr(x21))*
(x20 + x21)))*x38 + x30 =E= 0;
e13.. - (2.28 + 0.2*(0.0584999*x22 - 6.46057666666667e-5*sqr(x23 - x22) +
0.0584999*x23 - 0.0001938173*x22*x23 + 7.739465e-8*(sqr(x22) + sqr(x23))*
(x22 + x23)))*x39 + x31 =E= 0;
e14.. - (2.28 + 0.2*(0.0584999*x23 - 6.46057666666667e-5*sqr(x24 - x23) +
0.0584999*x24 - 0.0001938173*x23*x24 + 7.739465e-8*(sqr(x23) + sqr(x24))*
(x23 + x24)))*x40 + x32 =E= 0;
e15.. - (2.28 + 0.2*(0.0584999*x24 - 6.46057666666667e-5*sqr(x25 - x24) +
0.0584999*x25 - 0.0001938173*x24*x25 + 7.739465e-8*(sqr(x24) + sqr(x25))*
(x24 + x25)))*x41 + x33 =E= 0;
e16.. - (2.28 + 0.2*(0.0584999*x25 - 6.46057666666667e-5*sqr(x26 - x25) +
0.0584999*x26 - 0.0001938173*x25*x26 + 7.739465e-8*(sqr(x25) + sqr(x26))*
(x25 + x26)))*x42 + x34 =E= 0;
e17.. - x27 - x31 + x44 =E= 0;
e18.. - x28 - x32 + x45 =E= 0;
e19.. - x29 - x33 + x46 =E= 0;
e20.. - x30 - x34 + x47 =E= 0;
e21.. x43 + x49 =E= 0;
e22.. x44 + x50 =E= 280;
e23.. x45 + x51 =E= 250;
e24.. x46 + x52 =E= 275;
e25.. x47 + x53 =E= 250;
e26.. - x1 - x9 - x50 =E= -250;
e27.. - x2 - x10 - x51 =E= -250;
e28.. - x3 - x11 - x52 =E= -250;
e29.. - x4 - x12 - x53 =E= -250;
e30.. - x5 - x13 + x50 =E= 100;
e31.. - x6 - x14 + x51 =E= 100;
e32.. - x7 - x15 + x52 =E= 100;
e33.. - x8 - x16 + x53 =E= 100;
e34.. x1 - x44 =G= -45;
e35.. x2 - x45 =G= -45;
e36.. x3 - x46 =G= -45;
e37.. x4 - x47 =G= -45;
e38.. x5 + x44 =G= 35;
e39.. x6 + x45 =G= 35;
e40.. x7 + x46 =G= 35;
e41.. x8 + x47 =G= 35;
e42.. - (sqr(x50) + 10*x50 + sqr(x51) + 10*x51 + sqr(x52) + 10*x52 + sqr(x53)
+ 10*x53) + objvar =E= -44000;
* set non default bounds
x1.up = 45;
x2.up = 45;
x3.up = 45;
x4.up = 45;
x5.up = 65;
x6.up = 65;
x7.up = 65;
x8.up = 65;
x9.up = 105;
x10.up = 105;
x11.up = 105;
x12.up = 105;
x13.up = 85;
x14.up = 85;
x15.up = 85;
x16.up = 85;
x17.fx = 350;
x18.up = 354;
x19.up = 354;
x20.up = 354;
x21.fx = 350;
x22.fx = 158;
x23.up = 160;
x24.up = 160;
x25.up = 160;
x26.fx = 158;
x35.up = 5.184;
x36.up = 5.184;
x37.up = 5.184;
x38.up = 5.184;
x39.up = 3.456;
x40.up = 3.456;
x41.up = 3.456;
x42.up = 3.456;
x43.fx = 0;
x44.up = 75;
x45.up = 75;
x46.up = 75;
x47.up = 75;
x50.lo = 100; x50.up = 250;
x51.lo = 100; x51.up = 250;
x52.lo = 100; x52.up = 250;
x53.lo = 100; x53.up = 250;
* 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 NLP minimizing objvar;
]]></plaintext></body>