loadbal.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:34:05
*
* Equation counts
* Total E G L N X C
* 32 12 20 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 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 123 92 31 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,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;
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;
e1.. - 20*x1 - 80*x2 =G= -999.99;
e2.. - 80*x1 - 20*x2 =G= -999.99;
e3.. - 20*x3 - 80*x4 =G= -999.99;
e4.. - 80*x3 - 20*x4 =G= -999.99;
e5.. - 20*x5 - 80*x6 =G= -999.99;
e6.. - 80*x5 - 20*x6 =G= -999.99;
e7.. - 20*x9 - 80*x10 =G= -999.99;
e8.. - 80*x9 - 20*x10 =G= -999.99;
e9.. - 20*x11 - 80*x12 =G= -999.99;
e10.. - 80*x11 - 20*x12 =G= -999.99;
e11.. - 20*x15 - 80*x16 =G= -999.99;
e12.. - 80*x15 - 20*x16 =G= -999.99;
e13.. - 20*x17 - 80*x18 =G= -999.99;
e14.. - 80*x17 - 20*x18 =G= -999.99;
e15.. - 20*x19 - 80*x20 =G= -999.99;
e16.. - 80*x19 - 20*x20 =G= -999.99;
e17.. - 20*x7 - 80*x8 =G= -9999.99;
e18.. - 80*x7 - 20*x8 =G= -9999.99;
e19.. - 20*x13 - 80*x14 =G= -9999.99;
e20.. - 80*x13 - 20*x14 =G= -9999.99;
e21.. x1 - x2 - x21 =E= -95;
e22.. x3 - x4 - x22 =E= -95;
e23.. x5 - x6 - x23 =E= -19;
e24.. - x1 + x2 - x3 + x4 - x5 + x6 - x7 + x8 - x24 =E= -70;
e25.. x7 - x8 - x9 + x10 - x11 + x12 - x13 + x14 - x25 =E= -70;
e26.. x9 - x10 - x26 =E= -19;
e27.. x11 - x12 - x27 =E= -19;
e28.. x13 - x14 - x15 + x16 - x17 + x18 - x19 + x20 - x28 =E= -70;
e29.. x15 - x16 - x29 =E= -19;
e30.. x17 - x18 - x30 =E= -19;
e31.. x19 - x20 - x31 =E= -19;
e32.. - (0.00972762645914397*x21/(100 - x21) + 0.00972762645914397*x22/(100 -
x22) + 0.00972762645914397*x23/(20 - x23) + 0.00972762645914397*x24/(100
- x24) + 0.00972762645914397*x25/(100 - x25) + 0.00972762645914397*x26/(
20 - x26) + 0.00972762645914397*x27/(20 - x27) + 0.00972762645914397*x28/
(100 - x28) + 0.00972762645914397*x29/(20 - x29) + 0.00972762645914397*
x30/(20 - x30) + 0.00972762645914397*x31/(20 - x31) + 0.155642023346303*
x1/(1000 - 80*x1 - 20*x2) + 0.0389105058365759*x1/(1000 - 20*x1 - 80*x2)
+ 0.155642023346303*x3/(1000 - 80*x3 - 20*x4) + 0.0389105058365759*x3/(
1000 - 20*x3 - 80*x4) + 0.155642023346303*x5/(1000 - 80*x5 - 20*x6) +
0.0389105058365759*x5/(1000 - 20*x5 - 80*x6) + 0.155642023346303*x2/(1000
- 20*x1 - 80*x2) + 0.0389105058365759*x2/(1000 - 80*x1 - 20*x2) +
0.155642023346303*x4/(1000 - 20*x3 - 80*x4) + 0.0389105058365759*x4/(1000
- 80*x3 - 20*x4) + 0.155642023346303*x6/(1000 - 20*x5 - 80*x6) +
0.0389105058365759*x6/(1000 - 80*x5 - 20*x6) + 0.155642023346303*x15/(
1000 - 80*x15 - 20*x16) + 0.0389105058365759*x15/(1000 - 20*x15 - 80*x16)
+ 0.155642023346303*x17/(1000 - 80*x17 - 20*x18) + 0.0389105058365759*
x17/(1000 - 20*x17 - 80*x18) + 0.155642023346303*x19/(1000 - 80*x19 - 20*
x20) + 0.0389105058365759*x19/(1000 - 20*x19 - 80*x20) +
0.155642023346303*x16/(1000 - 20*x15 - 80*x16) + 0.0389105058365759*x16/(
1000 - 80*x15 - 20*x16) + 0.155642023346303*x18/(1000 - 20*x17 - 80*x18)
+ 0.0389105058365759*x18/(1000 - 80*x17 - 20*x18) + 0.155642023346303*
x20/(1000 - 20*x19 - 80*x20) + 0.0389105058365759*x20/(1000 - 80*x19 - 20
*x20) + 0.155642023346303*x9/(1000 - 80*x9 - 20*x10) + 0.0389105058365759
*x9/(1000 - 20*x9 - 80*x10) + 0.155642023346303*x10/(1000 - 20*x9 - 80*
x10) + 0.0389105058365759*x10/(1000 - 80*x9 - 20*x10) + 0.155642023346303
*x11/(1000 - 80*x11 - 20*x12) + 0.0389105058365759*x11/(1000 - 20*x11 -
80*x12) + 0.155642023346303*x12/(1000 - 20*x11 - 80*x12) +
0.0389105058365759*x12/(1000 - 80*x11 - 20*x12) + 0.155642023346303*x8/(
10000 - 20*x7 - 80*x8) + 0.0389105058365759*x8/(10000 - 80*x7 - 20*x8) +
0.155642023346303*x7/(10000 - 80*x7 - 20*x8) + 0.0389105058365759*x7/(
10000 - 20*x7 - 80*x8) + 0.155642023346303*x13/(10000 - 80*x13 - 20*x14)
+ 0.0389105058365759*x13/(10000 - 20*x13 - 80*x14) + 0.155642023346303*
x14/(10000 - 20*x13 - 80*x14) + 0.0389105058365759*x14/(10000 - 80*x13 -
20*x14)) + objvar =E= 0;
* set non default bounds
x21.up = 99.99;
x22.up = 99.99;
x23.up = 19.99;
x24.up = 99.99;
x25.up = 99.99;
x26.up = 19.99;
x27.up = 19.99;
x28.up = 99.99;
x29.up = 19.99;
x30.up = 19.99;
x31.up = 19.99;
* set non default levels
x21.l = 95;
x22.l = 95;
x23.l = 19;
x24.l = 70;
x25.l = 70;
x26.l = 19;
x27.l = 19;
x28.l = 70;
x29.l = 19;
x30.l = 19;
x31.l = 19;
* 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>