kowosb.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:33:49
*
* 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
* 5 5 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 5 1 4 0
*
* Solve m using NLP minimizing objvar;
Variables x1,x2,x3,x4,objvar;
Equations e1;
e1.. - (sqr(0.1957 - x1*(16 + 4*x2)/(16 + 4*x3 + x4)) + sqr(0.1947 - x1*(4 + 2
*x2)/(4 + 2*x3 + x4)) + sqr(0.1735 - x1*(1 + x2)/(1 + x3 + x4)) + sqr(0.16
- x1*(0.25 + 0.5*x2)/(0.25 + 0.5*x3 + x4)) + sqr(0.0844 - x1*(0.0625 +
0.25*x2)/(0.0625 + 0.25*x3 + x4)) + sqr(0.0627 - x1*(0.027889 + 0.167*x2)/
(0.027889 + 0.167*x3 + x4)) + sqr(0.0456 - x1*(0.015625 + 0.125*x2)/(
0.015625 + 0.125*x3 + x4)) + sqr(0.0342 - x1*(0.01 + 0.1*x2)/(0.01 + 0.1*
x3 + x4)) + sqr(0.0323 - x1*(0.00693889 + 0.0833*x2)/(0.00693889 + 0.0833*
x3 + x4)) + sqr(0.0235 - x1*(0.00509796 + 0.0714*x2)/(0.00509796 + 0.0714*
x3 + x4)) + sqr(0.0246 - x1*(0.00390625 + 0.0625*x2)/(0.00390625 + 0.0625*
x3 + x4))) + objvar =E= 0;
* set non default bounds
* set non default levels
x1.l = 0.25;
x2.l = 0.39;
x3.l = 0.415;
x4.l = 0.39;
* 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>