djtl.gms
<![CDATA[
* DNLP written by GAMS Convert at 10/06/06 11:31:19
*
* 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
* 3 3 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 3 1 2 0
*
* Solve m using DNLP minimizing objvar;
Variables x1,x2,objvar;
Equations e1;
e1.. - 0.001*(POWER((-10) + x1,3) + POWER((-20) + x2,3) + 10000000000*(1 -
max(0,201 + (-sqr((-5) + x1)) - sqr((-5) + x2))/(1e-10 + max(0,201 + (-
sqr((-5) + x1)) - sqr((-5) + x2))))*sqr(200 + (-sqr((-5) + x1)) - sqr((-5)
+ x2)) - max(0,201 + (-sqr((-5) + x1)) - sqr((-5) + x2))/(1e-10 + max(0,
201 + (-sqr((-5) + x1)) - sqr((-5) + x2)))*log(1e-10 + abs(201 + (-sqr((-5
) + x1)) - sqr((-5) + x2))) + 10000000000*(1 - max(0,(-99) + sqr((-5) + x1
) + sqr((-5) + x2))/(1e-10 + max(0,(-99) + sqr((-5) + x1) + sqr((-5) + x2)
)))*sqr((-100) + sqr((-5) + x1) + sqr((-5) + x2)) - max(0,(-99) + sqr((-5)
+ x1) + sqr((-5) + x2))/(1e-10 + max(0,(-99) + sqr((-5) + x1) + sqr((-5)
+ x2)))*log(1e-10 + abs((-99) + sqr((-5) + x1) + sqr((-5) + x2))) +
10000000000*(1 - max(0,1 + sqr((-5) + x2) + sqr((-6) + x1))/(1e-10 + max(0
,1 + sqr((-5) + x2) + sqr((-6) + x1))))*sqr(sqr((-5) + x2) + sqr((-6) + x1
)) - max(0,1 + sqr((-5) + x2) + sqr((-6) + x1))/(1e-10 + max(0,1 + sqr((-5
) + x2) + sqr((-6) + x1)))*log(1e-10 + abs(1 + sqr((-5) + x2) + sqr((-6)
+ x1))) + 10000000000*(1 - max(0,83.81 + (-sqr((-5) + x2)) - sqr((-6) +
x1))/(1e-10 + max(0,83.81 + (-sqr((-5) + x2)) - sqr((-6) + x1))))*sqr(
82.81 + (-sqr((-5) + x2)) - sqr((-6) + x1)) - max(0,83.81 + (-sqr((-5) +
x2)) - sqr((-6) + x1))/(1e-10 + max(0,83.81 + (-sqr((-5) + x2)) - sqr((-6)
+ x1)))*log(1e-10 + abs(83.81 + (-sqr((-5) + x2)) - sqr((-6) + x1))) +
10000000000*(1 - max(0,101 - x1)/(1e-10 + max(0,101 - x1)))*sqr(100 - x1)
- max(0,101 - x1)/(1e-10 + max(0,101 - x1))*log(1e-10 + abs(101 - x1)) +
10000000000*(1 - max(0,(-12) + x1)/(1e-10 + max(0,(-12) + x1)))*sqr((-13)
+ x1) - max(0,(-12) + x1)/(1e-10 + max(0,(-12) + x1))*log(1e-10 + abs((-
12) + x1)) + 10000000000*(1 - max(0,101 - x2)/(1e-10 + max(0,101 - x2)))*
sqr(100 - x2) - max(0,101 - x2)/(1e-10 + max(0,101 - x2))*log(1e-10 + abs(
101 - x2)) + 10000000000*(1 - max(0,1 + x2)/(1e-10 + max(0,1 + x2)))*sqr(
x2) - max(0,1 + x2)/(1e-10 + max(0,1 + x2))*log(1e-10 + abs(1 + x2)))
+ objvar =E= 0;
* set non default bounds
* set non default levels
x1.l = 15;
x2.l = -1;
* set non default marginals
Model m / all /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include '%gams.u1%'
Solve m using DNLP minimizing objvar;
]]></plaintext></body>