hs070.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:33:42
*
* Equation counts
* Total E G L N X C
* 22 21 1 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 25 25 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 139 21 118 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,objvar;
Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19
,e20,e21,e22;
e1.. (1 - x3)*x4 + x3 =G= 0;
e2.. - ((1 - x3)*x4 + x3) + x5 =E= 0;
e3.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(4.3418554699846 - 4.3418554699846*x2)*exp(x2 - 0.0130582397492818*
x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*(1 - x3)*sqrt(
0.159154570919277*x1)*(1 - 0.0130582397492818*x5*x1/x4)) + x6 =E= 0;
e4.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(2.03927037699056 - 2.03927037699056*x2)*exp(x2 - 0.130582397492818
*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*(1 - x3)*sqrt(
0.159154570919277*x1)*(1 - 0.130582397492818*x5*x1/x4)) + x7 =E= 0;
e5.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(1.34612319643061 - 1.34612319643061*x2)*exp(x2 - 0.261164794985636
*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*(1 - x3)*sqrt(
0.159154570919277*x1)*(1 - 0.261164794985636*x5*x1/x4)) + x8 =E= 0;
e6.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(0.940658088322449 - 0.940658088322449*x2)*exp(x2 -
0.391747192478454*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 0.391747192478454*x5*x1/x4)) + x9
=E= 0;
e7.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(0.652976015870668 - 0.652976015870668*x2)*exp(x2 -
0.522329589971272*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 0.522329589971272*x5*x1/x4))
+ x10 =E= 0;
e8.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(0.429832464556459 - 0.429832464556459*x2)*exp(x2 -
0.65291198746409*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*(
1 - x3)*sqrt(0.159154570919277*x1)*(1 - 0.65291198746409*x5*x1/x4)) + x11
=E= 0;
e9.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277*
x2)*exp(0.247510907762504 - 0.247510907762504*x2)*exp(x2 -
0.783494384956908*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 0.783494384956908*x5*x1/x4))
+ x12 =E= 0;
e10.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp(0.0933602279352459 - 0.0933602279352459*x2)*exp(x2 -
0.914076782449726*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)
*(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 0.914076782449726*x5*x1/x4))
+ x13 =E= 0;
e11.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.0401711646892769) + 0.0401711646892769*x2)*exp(x2 -
1.04465917994254*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.04465917994254*x5*x1/x4))
+ x14 =E= 0;
e12.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.15795420034566) + 0.15795420034566*x2)*exp(x2 -
1.17524157743536*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.17524157743536*x5*x1/x4))
+ x15 =E= 0;
e13.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.263314716003487) + 0.263314716003487*x2)*exp(x2 -
1.30582397492818*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.30582397492818*x5*x1/x4))
+ x16 =E= 0;
e14.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.358624895807811) + 0.358624895807811*x2)*exp(x2 -
1.436406372421*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*(1
- x3)*sqrt(0.159154570919277*x1)*(1 - 1.436406372421*x5*x1/x4)) + x17
=E= 0;
e15.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.445636272797441) + 0.445636272797441*x2)*exp(x2 -
1.56698876991382*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.56698876991382*x5*x1/x4))
+ x18 =E= 0;
e16.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.525678980470977) + 0.525678980470977*x2)*exp(x2 -
1.69757116740663*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.69757116740663*x5*x1/x4))
+ x19 =E= 0;
e17.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.599786952624699) + 0.599786952624699*x2)*exp(x2 -
1.82815356489945*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.82815356489945*x5*x1/x4))
+ x20 =E= 0;
e18.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.668779824111651) + 0.668779824111651*x2)*exp(x2 -
1.95873596239227*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 1.95873596239227*x5*x1/x4))
+ x21 =E= 0;
e19.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.733318345249222) + 0.733318345249222*x2)*exp(x2 -
2.08931835988509*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 2.08931835988509*x5*x1/x4))
+ x22 =E= 0;
e20.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.793942967065657) + 0.793942967065657*x2)*exp(x2 -
2.21990075737791*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 2.21990075737791*x5*x1/x4))
+ x23 =E= 0;
e21.. - ((1 + 0.0833333333333333/x2)*exp(log(x3*x5)*x2)*sqrt(0.159154570919277
*x2)*exp((-0.851101380905606) + 0.851101380905606*x2)*exp(x2 -
2.35048315487072*x5*x2) + (1 + 0.0833333333333333/x1)*exp(log(x5/x4)*x1)*
(1 - x3)*sqrt(0.159154570919277*x1)*(1 - 2.35048315487072*x5*x1/x4))
+ x24 =E= 0;
e22.. - (sqr((-0.00189) + x6) + sqr((-0.1038) + x7) + sqr((-0.268) + x8) +
sqr((-0.506) + x9) + sqr((-0.577) + x10) + sqr((-0.604) + x11) + sqr((-
0.725) + x12) + sqr((-0.898) + x13) + sqr((-0.947) + x14) + sqr((-0.845)
+ x15) + sqr((-0.702) + x16) + sqr((-0.528) + x17) + sqr((-0.385) + x18)
+ sqr((-0.257) + x19) + sqr((-0.159) + x20) + sqr((-0.0869) + x21) +
sqr((-0.0453) + x22) + sqr((-0.01509) + x23) + sqr((-0.00189) + x24))
+ objvar =E= 0;
* set non default bounds
x1.lo = 1E-5; x1.up = 100;
x2.lo = 1E-5; x2.up = 100;
x3.lo = 1E-5; x3.up = 1;
x4.lo = 1E-5; x4.up = 100;
x5.lo = 1E-5;
* set non default levels
x1.l = 2;
x2.l = 4;
x3.l = 0.04;
x4.l = 2;
* 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>