GLOBAL World [ GLOBAL World Home | GLOBALLib | Contact ]

ex6_2_10.gms:

References:

Point:


* NLP written by GAMS Convert at 07/19/01 13:39:43 * * Equation counts * Total E G L N X * 4 4 0 0 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 7 7 0 0 0 0 0 0 * FX 0 0 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 13 7 6 0 * * Solve m using NLP minimizing objvar; Variables objvar,x2,x3,x4,x5,x6,x7; Equations e1,e2,e3,e4; e1.. - ((15.3261663216011*x2 + 23.2043471859416*x4 + 6.69678129464404*x6)*log( 2.1055*x2 + 3.1878*x4 + 0.92*x6) - 2.46348749603266*x2 - 4.33155441248417* x4 - 0.626542690017204*x6 + 6.4661663216011*x2*log(x2/(2.1055*x2 + 3.1878* x4 + 0.92*x6)) + 12.2043471859416*x4*log(x4/(2.1055*x2 + 3.1878*x4 + 0.92* x6)) + 0.696781294644034*x6*log(x6/(2.1055*x2 + 3.1878*x4 + 0.92*x6)) + 9.86*x2*log(x2/(1.972*x2 + 2.4*x4 + 1.4*x6)) + 12*x4*log(x4/(1.972*x2 + 2.4*x4 + 1.4*x6)) + 7*x6*log(x6/(1.972*x2 + 2.4*x4 + 1.4*x6)) + (1.972*x2 + 2.4*x4 + 1.4*x6)*log(1.972*x2 + 2.4*x4 + 1.4*x6) + 1.972*x2*log(x2/( 1.972*x2 + 0.283910843616504*x4 + 3.02002220174195*x6)) + 2.4*x4*log(x4/( 1.45991339466884*x2 + 2.4*x4 + 0.415073537580851*x6)) + 1.4*x6*log(x6/( 0.602183324335333*x2 + 0.115623371371275*x4 + 1.4*x6)) + (15.3261663216011 *x3 + 23.2043471859416*x5 + 6.69678129464404*x7)*log(2.1055*x3 + 3.1878*x5 + 0.92*x7) - 2.46348749603266*x3 - 4.33155441248417*x5 - 0.626542690017204*x7 + 6.4661663216011*x3*log(x3/(2.1055*x3 + 3.1878*x5 + 0.92*x7)) + 12.2043471859416*x5*log(x5/(2.1055*x3 + 3.1878*x5 + 0.92*x7)) + 0.696781294644034*x7*log(x7/(2.1055*x3 + 3.1878*x5 + 0.92*x7)) + 9.86* x3*log(x3/(1.972*x3 + 2.4*x5 + 1.4*x7)) + 12*x5*log(x5/(1.972*x3 + 2.4*x5 + 1.4*x7)) + 7*x7*log(x7/(1.972*x3 + 2.4*x5 + 1.4*x7)) + (1.972*x3 + 2.4* x5 + 1.4*x7)*log(1.972*x3 + 2.4*x5 + 1.4*x7) + 1.972*x3*log(x3/(1.972*x3 + 0.283910843616504*x5 + 3.02002220174195*x7)) + 2.4*x5*log(x5/( 1.45991339466884*x3 + 2.4*x5 + 0.415073537580851*x7)) + 1.4*x7*log(x7/( 0.602183324335333*x3 + 0.115623371371275*x5 + 1.4*x7)) - 17.2981663216011* x2*log(x2) - 25.6043471859416*x4*log(x4) - 8.09678129464404*x6*log(x6) - 17.2981663216011*x3*log(x3) - 25.6043471859416*x5*log(x5) - 8.09678129464404*x7*log(x7)) + objvar =E= 0; e2.. x2 + x3 =E= 0.2; e3.. x4 + x5 =E= 0.4; e4.. x6 + x7 =E= 0.4; * set non default bounds x2.lo = 1E-7; x2.up = 0.2; x3.lo = 1E-7; x3.up = 0.2; x4.lo = 1E-7; x4.up = 0.4; x5.lo = 1E-7; x5.up = 0.4; x6.lo = 1E-7; x6.up = 0.4; x7.lo = 1E-7; x7.up = 0.4; * set non default levels x2.l = 0.19863; x3.l = 0.00137; x4.l = 0.00428; x5.l = 0.39572; x6.l = 0.39922; x7.l = 0.00078; * 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;