ex8_5_3.gms:
References:
- Floudas, C A, Pardalos, P M, Adjiman, C S, Esposito, W R, Gumus, Z H, Harding, S T, Klepeis, J L, Meyer, C A, and Schweiger, C A, Handbook of Test Problems in Local and Global Optimization. Kluwer Academic Publishers, 1999.
- Michelsen, M L, The Isothermal Flash Problem - Part i. Stability. Fluid Phase Equilibria 9 (1982), 1.
- Original source: Global Model of Chapter 8 ex8.5.3.gms from Floudas e.a. Test Problems
Point:
p1
Best known point: p1 with value -0.0041
* NLP written by GAMS Convert at 07/19/01 13:40:14
*
* Equation counts
* Total E G L N X
* 5 5 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 6 6 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 17 7 10 0
*
* Solve m using NLP minimizing objvar;
Variables objvar,x2,x3,x4,x5,x6;
Equations e1,e2,e3,e4,e5;
e1.. - (x2*log(x2) + x3*log(x3) - log(x4 - x6) + x4 - x5*log(1 + x6/x4)/x6 +
5.0464317551216*x2 + 0.366877055769689*x3) + objvar =E= -1;
e2.. POWER(x4,3) - sqr(x4) + (-sqr(x6) - x6 + x5)*x4 - x5*x6 =E= 0;
e3.. - (1.04633*x2*x2 + 0.579822*x2*x3 + 0.579822*x3*x2 + 0.379615*x3*x3) + x5
=E= 0;
e4.. - 0.0771517*x2 - 0.0765784*x3 + x6 =E= 0;
e5.. x2 + x3 =E= 1;
* set non default bounds
* set non default levels
x2.l = 0.5;
x3.l = 0.5;
x4.l = 2;
x5.l = 1;
x6.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 NLP minimizing objvar;