GLOBAL World [ GLOBAL World Home | GLOBALLib | Contact ]

ex6_2_12.gms:

References:

Point:


* NLP written by GAMS Convert at 07/19/01 13:39:44 * * Equation counts * Total E G L N X * 3 3 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 * 9 5 4 0 * * Solve m using NLP minimizing objvar; Variables objvar,x2,x3,x4,x5; Equations e1,e2,e3; e1.. - (x2*log(x2/(8*x2 + x4)) + x4*log(x4/(8*x2 + x4)) + 0.0696225416798359* x2 + 0.752006*x4 + (8*x2 + 1.6*x4)*log(8*x2 + 1.6*x4) + 5*x2*log(x2/( 5.00000397494442*x2 + 0.480353357956269*x4)) + 3*x2*log(x2/( 8.96062592375197*x2 + 1.13841069150863*x4)) + 1.6*x4*log(x4/( 1.69889877049372*x2 + 1.6*x4)) + x3*log(x3/(8*x3 + x5)) + x5*log(x5/(8*x3 + x5)) + 0.0696225416798359*x3 + 0.752006*x5 + (8*x3 + 1.6*x5)*log(8*x3 + 1.6*x5) + 5*x3*log(x3/(5.00000397494442*x3 + 0.480353357956269*x5)) + 3 *x3*log(x3/(8.96062592375197*x3 + 1.13841069150863*x5)) + 1.6*x5*log(x5/( 1.69889877049372*x3 + 1.6*x5)) - 8*x2*log(x2) - 1.6*x4*log(x4) - 8*x3*log( x3) - 1.6*x5*log(x5)) + objvar =E= 0; e2.. x2 + x3 =E= 0.5; e3.. x4 + x5 =E= 0.5; * set non default bounds x2.lo = 1E-7; x2.up = 0.5; x3.lo = 1E-7; x3.up = 0.5; x4.lo = 1E-7; x4.up = 0.5; x5.lo = 1E-7; x5.up = 0.5; * set non default levels x2.l = 0.4994; x3.l = 0.0006; x4.l = 0.1179; x5.l = 0.3821; * 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;