hs99exp.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:33:46
*
* Equation counts
* Total E G L N X C
* 22 22 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 32 32 0 0 0 0 0 0
* FX 3 3 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 72 50 22 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,objvar;
Positive Variables x1,x2,x3,x4,x5,x9,x13,x17,x21,x25;
Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19
,e20,e21,e22;
e1.. 1250*cos(x1) + x2 - x6 =E= 0;
e2.. 15625*sin(x1) + x3 + 25*x4 - x7 =E= 10000;
e3.. 1250*sin(x1) + x4 - x8 =E= 800;
e4.. 1250*cos(x5) + x6 - x10 =E= 0;
e5.. 15625*sin(x5) + x7 + 25*x8 - x11 =E= 10000;
e6.. 1250*sin(x5) + x8 - x12 =E= 800;
e7.. 3750*cos(x9) + x10 - x14 =E= 0;
e8.. 93750*sin(x9) + x11 + 50*x12 - x15 =E= 40000;
e9.. 3750*sin(x9) + x12 - x16 =E= 1600;
e10.. 3750*cos(x13) + x14 - x18 =E= 0;
e11.. 93750*sin(x13) + x15 + 50*x16 - x19 =E= 40000;
e12.. 3750*sin(x13) + x16 - x20 =E= 1600;
e13.. 3750*cos(x17) + x18 - x22 =E= 0;
e14.. 93750*sin(x17) + x19 + 50*x20 - x23 =E= 40000;
e15.. 3750*sin(x17) + x20 - x24 =E= 1600;
e16.. 9000*cos(x21) + x22 - x26 =E= 0;
e17.. 405000*sin(x21) + x23 + 90*x24 - x27 =E= 129600;
e18.. 9000*sin(x21) + x24 - x28 =E= 2880;
e19.. 9000*cos(x25) + x26 - x29 =E= 0;
e20.. 405000*sin(x25) + x27 + 90*x28 - x30 =E= 100000;
e21.. 9000*sin(x25) + x28 - x31 =E= 1000;
e22.. x29*x29 + objvar =E= 0;
* set non default bounds
x1.up = 1.58;
x2.fx = 0;
x3.fx = 0;
x4.fx = 0;
x5.up = 1.58;
x9.up = 1.58;
x13.up = 1.58;
x17.up = 1.58;
x21.up = 1.58;
x25.up = 1.58;
* set non default levels
x1.l = 0.5;
x5.l = 0.5;
x9.l = 0.5;
x13.l = 0.5;
x17.l = 0.5;
x21.l = 0.5;
x25.l = 0.5;
* 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>