\$offlisting * * Equation counts * Total E G L N X C B * 5 5 0 0 0 0 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 12 12 0 0 0 0 0 0 * FX 0 * * Nonzero counts * Total const NL DLL * 37 26 11 0 * * Solve m using NLP minimizing objvar; Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,objvar; Equations e1,e2,e3,e4,e5; e1.. x1 + 2*x2 + 2*x3 + x6 + x10 =E= 2; e2.. x4 + 2*x5 + x6 + x7 =E= 1; e3.. x3 + x7 + x8 + 2*x9 + x10 =E= 1; e4.. -((-6.05576803624071 + log(x1/x11))*x1 + (-17.1307680362407 + log(x2/x11)) *x2 + (-34.0207680362407 + log(x3/x11))*x3 + (-5.88076803624071 + log(x4/ x11))*x4 + (-24.6877680362407 + log(x5/x11))*x5 + (-14.9527680362407 + log(x6/x11))*x6 + (-24.0667680362407 + log(x7/x11))*x7 + (- 10.6747680362407 + log(x8/x11))*x8 + (-26.6287680362407 + log(x9/x11))*x9 + (-22.1447680362407 + log(x10/x11))*x10) + objvar =E= 0; e5.. - x1 - x2 - x3 - x4 - x5 - x6 - x7 - x8 - x9 - x10 + x11 =E= 0; * set non-default bounds x1.lo = 0.001; x2.lo = 0.001; x3.lo = 0.001; x4.lo = 0.001; x5.lo = 0.001; x6.lo = 0.001; x7.lo = 0.001; x8.lo = 0.001; x9.lo = 0.001; x10.lo = 0.001; x11.lo = 0.01; Model m / all /; m.limrow=0; m.limcol=0; m.tolproj=0.0; \$if NOT '%gams.u1%' == '' \$include '%gams.u1%' \$if not set NLP \$set NLP NLP Solve m using %NLP% minimizing objvar;