launch.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:33:50
*
* Equation counts
* Total E G L N X C
* 30 7 9 14 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 27 27 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 144 80 64 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,objvar;
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;
e1.. 2*x1 - x2 =E= 0;
e2.. - x8 + 0.6*x9 =E= 0;
e3.. - x15 + 0.7*x16 =E= 0;
e4.. - x4 + 5*x6 =E= 0;
e5.. - x11 + 5*x13 =E= 0;
e6.. x18 - x20 =E= 0;
e7.. 12*x2 - x5 =L= 0;
e8.. 17*x2 - x5 =G= 0;
e9.. 10*x9 - x12 =L= 0;
e10.. 13*x9 - x12 =G= 0;
e11.. 7*x16 - x19 =L= 0;
e12.. 10*x16 - x19 =G= 0;
e13.. 1.2*x2 - x4 + 1.2*x5 + 1.2*x9 + 1.2*x12 + 1.2*x16 + 1.2*x19 + 1.2*x22
=L= -24;
e14.. 1.4*x2 - x4 + 1.4*x5 + 1.4*x9 + 1.4*x12 + 1.4*x16 + 1.4*x19 + 1.4*x22
=G= -28;
e15.. 0.6*x9 - x11 + 0.6*x12 + 0.6*x16 + 0.6*x19 + 0.6*x22 =L= -12;
e16.. 0.75*x9 - x11 + 0.75*x12 + 0.75*x16 + 0.75*x19 + 0.75*x22 =G= -15;
e17.. 0.7*x16 - x18 + 0.7*x19 + 0.7*x22 =L= -14;
e18.. 0.9*x16 - x18 + 0.9*x19 + 0.9*x22 =G= -18;
e19.. x3*x22 + 20*x3 - x22 - x2 - x9 - x12 - x16 - x19 - x26 =L= 20;
e20.. x10*x22 + 20*x10 - x22 - x9 - x16 - x19 - x26 =L= 20;
e21.. x17*x22 + 20*x17 - x22 - x16 - x26 =L= 20;
e22.. - x4*x23 + 240*x5 =L= 0;
e23.. - x4*x23 + 290*x5 =G= 0;
e24.. - x11*x24 + 240*x12 =L= 0;
e25.. - x11*x24 + 290*x12 =G= 0;
e26.. - x18*x25 + 340*x19 =L= 0;
e27.. - x18*x25 + 375*x19 =G= 0;
e28.. - (-32*x4*x23*log(x3)/x5 - 32*x11*x24*log(x10)/x12 - 32*x18*x25*log(x17)
/x19) =L= -35000;
e29.. (-32*x4*x23*log(x3)/x5) - 32*x11*x24*log(x10)/x12 - 32*x18*x25*log(x17)/
x19 =L= 50000;
e30.. - (x26*x26 + 1e-8*(5272.77*exp(1.2781*log(x1))*exp(-0.1959*log(x2))*exp(
2.4242*log(x3))*exp(0.38745*log(x4))*exp(0.9904*log(x5)) + 160.909*exp(-
0.146*log(0.001*x6)) + 282.874*exp(0.648*log(0.001*x6)) + 0.64570846*exp(
0.3322*log(x1))*exp(-1.5935*log(x3))*exp(0.2363*log(x5))*exp(0.1079*log(
x7)) + 31.136196*exp(0.736*log(0.001*x6)) + 12.092112*exp(-0.229*log(
0.001*x6)) + 31.136196*exp(0.736*log(0.001*x13)) + 12.092112*exp(-0.229*
log(0.001*x13)) + 0.0002587*x6) + 1e-8*(5272.77*exp(1.2781*log(x8))*exp(-
0.1959*log(x9))*exp(2.4242*log(x10))*exp(0.38745*log(x11))*exp(0.9904*
log(x12)) + 160.909*exp(-0.146*log(0.001*x13)) + 282.874*exp(0.648*log(
0.001*x13)) + 0.64570846*exp(0.3322*log(x8))*exp(-1.5935*log(x10))*exp(
0.2363*log(x12))*exp(0.1079*log(x14)) + 0.0002587*x13) + 1e-8*(5272.77*
exp(1.2781*log(x15))*exp(-0.1959*log(x16))*exp(2.4242*log(x17))*exp(
0.38745*log(x18))*exp(0.9904*log(x19)) + 181.806*exp(0.539*log(0.001*x20)
) + 232.57*exp(0.772*log(0.001*x20)) + 0.49783215*exp(0.3322*log(x15))*
exp(-1.5935*log(x17))*exp(0.2363*log(x19))*exp(0.1079*log(x21)) -
0.22424514*exp(-1.33*log(0.01*x20)) + 20.708238*exp(0.498*log(0.01*x20))
+ 0.001958*x20) + 2.55000001785e-8*exp(0.46*log(0.003*x5 + 0.003*x12 +
0.003*x19))) - 4.7040096E-7*x22 + objvar =E= -5.64017E-6;
* set non default bounds
x1.lo = 1E-8; x1.up = 10000;
x2.lo = 1E-8; x2.up = 10000;
x3.lo = 0.25; x3.up = 0.3;
x4.lo = 1E-8; x4.up = 10000;
x5.lo = 1E-8; x5.up = 10000;
x6.lo = 1E-8; x6.up = 10000;
x7.lo = 125; x7.up = 150;
x8.lo = 1E-8; x8.up = 10000;
x9.lo = 1E-8; x9.up = 10000;
x10.lo = 0.24; x10.up = 0.29;
x11.lo = 1E-8; x11.up = 10000;
x12.lo = 1E-8; x12.up = 10000;
x13.lo = 1E-8; x13.up = 10000;
x14.lo = 75; x14.up = 100;
x15.lo = 1E-8; x15.up = 10000;
x16.lo = 1E-8; x16.up = 10000;
x17.lo = 0.16; x17.up = 0.21;
x18.lo = 1E-8; x18.up = 10000;
x19.lo = 1E-8; x19.up = 10000;
x20.lo = 1E-8; x20.up = 10000;
x21.lo = 50; x21.up = 70;
x22.lo = 2.5; x22.up = 4;
x23.lo = 1E-8; x23.up = 10000;
x24.lo = 1E-8; x24.up = 10000;
x25.lo = 1E-8; x25.up = 10000;
* set non default levels
x1.l = 68;
x2.l = 136;
x3.l = 0.29988744;
x4.l = 3733;
x5.l = 2177;
x6.l = 746.6;
x8.l = 28.2;
x9.l = 47;
x10.l = 0.28939109;
x11.l = 480;
x12.l = 566;
x13.l = 96;
x15.l = 11.2;
x16.l = 16;
x17.l = 0.20980926;
x18.l = 129;
x19.l = 145;
x20.l = 129;
x23.l = 155;
x24.l = 314;
x25.l = 403;
* 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>