Paviani.gms
<![CDATA[
* NLP written by GAMS Convert at 11/11/08 09:34:02
*
* Equation counts
* Total E G L N X C
* 1 1 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 11 11 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 11 1 10 0
*
* Solve m using NLP minimizing objvar;
*
Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,objvar;
Equations e1;
e1.. - (sqr(log((-2) + x1)) + sqr(log(10 - x1)) + sqr(log((-2) + x2)) + sqr(
log(10 - x2)) + sqr(log((-2) + x3)) + sqr(log(10 - x3)) + sqr(log((-2) +
x4)) + sqr(log(10 - x4)) + sqr(log((-2) + x5)) + sqr(log(10 - x5)) + sqr(
log((-2) + x6)) + sqr(log(10 - x6)) + sqr(log((-2) + x7)) + sqr(log(10 -
x7)) + sqr(log((-2) + x8)) + sqr(log(10 - x8)) + sqr(log((-2) + x9)) +
sqr(log(10 - x9)) + sqr(log((-2) + x10)) + sqr(log(10 - x10)) - exp(0.2*
log(x1*x2*x3*x4*x5*x6*x7*x8*x9*x10))) + objvar =E= 0;
* set non default bounds
x1.lo = 2.000001; x1.up = 9.999999;
x2.lo = 2.000001; x2.up = 9.999999;
x3.lo = 2.000001; x3.up = 9.999999;
x4.lo = 2.000001; x4.up = 9.999999;
x5.lo = 2.000001; x5.up = 9.999999;
x6.lo = 2.000001; x6.up = 9.999999;
x7.lo = 2.000001; x7.up = 9.999999;
x8.lo = 2.000001; x8.up = 9.999999;
x9.lo = 2.000001; x9.up = 9.999999;
x10.lo = 2.000001; x10.up = 9.999999;
* set non default levels
x1.l = 2.000001;
x2.l = 2.000001;
x3.l = 2.000001;
x4.l = 2.000001;
x5.l = 2.000001;
x6.l = 2.000001;
x7.l = 2.000001;
x8.l = 2.000001;
x9.l = 2.000001;
x10.l = 2.000001;
* set non default marginals
$if not set startsol $goto cont
$goto %startsol%
$label sol1
x1.l = 9.351 ;
x2.l = 9.351 ;
x3.l = 9.351 ;
x4.l = 9.351 ;
x5.l = 9.351 ;
x6.l = 9.351 ;
x7.l = 9.351 ;
x8.l = 9.351 ;
x9.l = 9.351 ;
x10.l = 9.351 ;
objvar.l = -45.778 ;
$goto cont
$label cont
Model m / all /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include '%gams.u1%'
Solve m using NLP minimizing objvar;
]]></plaintext></body>