ex14_2_9.gms:
Reference:
- 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.
- Original source: Global Model of Chapter 14 ex14.2.9.gms from Floudas e.a. Test Problems
Point:
p1
Best known point: p1 with value 0.0000
* NLP written by GAMS Convert at 07/19/01 13:40:31
*
* Equation counts
* Total E G L N X
* 6 2 0 4 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
* 20 8 12 0
*
* Solve m using NLP minimizing objvar;
Variables x1,x2,x3,objvar,x5;
Positive Variables x5;
Equations e1,e2,e3,e4,e5,e6;
e1.. objvar - x5 =E= 0;
e2.. 8.86*log(2.1055*x1 + 4.0456*x2) - 7.888*log(1.972*x1 + 3.236*x2) - (
2.1105532*x2 - 0.922208999999999*x1)/(2.1055*x1 + 4.0456*x2) - (0.848*log(
1.52337552625369*x1 + 3.236*x2) + 1.124*log(1.17581829697036*x1 +
0.197740576646344*x2)) - ((1.29182244626313*x1 + 1.29182244626313*x2)/(
1.52337552625369*x1 + 3.236*x2) + 3.29049113670798*x2/(1.52337552625369*x1
+ 3.236*x2) + 0.347329619985842*x2/(1.52337552625369*x1 + 3.236*x2) +
1.32161976579469*x1/(1.17581829697036*x1 + 0.197740576646344*x2)) -
3803.98/(231.47 + x3) - x5 =L= -13.1111702786953;
e3.. 15.18*log(2.1055*x1 + 4.0456*x2) - 12.944*log(1.972*x1 + 3.236*x2) - (
4.05530944*x2 - 1.7719728*x1)/(2.1055*x1 + 4.0456*x2) - (0.848*log(
1.52337552625369*x1 + 3.236*x2) + 2.16*log(1.52337552625369*x1 + 3.236*x2)
+ 0.228*log(1.52337552625369*x1 + 3.236*x2)) - ((2.744128*x1 + 2.744128*
x2)/(1.52337552625369*x1 + 3.236*x2) + 6.98976*x2/(1.52337552625369*x1 +
3.236*x2) + 0.737808*x2/(1.52337552625369*x1 + 3.236*x2) +
0.222260408150491*x1/(1.17581829697036*x1 + 0.197740576646344*x2)) -
2735.58621973158/(226.276 + x3) - x5 =L= -11.2003192377536;
e4.. 7.888*log(1.972*x1 + 3.236*x2) - 8.86*log(2.1055*x1 + 4.0456*x2) + (
2.1105532*x2 - 0.922208999999999*x1)/(2.1055*x1 + 4.0456*x2) + 0.848*log(
1.52337552625369*x1 + 3.236*x2) + 1.124*log(1.17581829697036*x1 +
0.197740576646344*x2) + (1.29182244626313*x1 + 1.29182244626313*x2)/(
1.52337552625369*x1 + 3.236*x2) + 3.29049113670798*x2/(1.52337552625369*x1
+ 3.236*x2) + 0.347329619985842*x2/(1.52337552625369*x1 + 3.236*x2) +
1.32161976579469*x1/(1.17581829697036*x1 + 0.197740576646344*x2) + 3803.98
/(231.47 + x3) - x5 =L= 13.1111702786953;
e5.. 12.944*log(1.972*x1 + 3.236*x2) - 15.18*log(2.1055*x1 + 4.0456*x2) + (
4.05530944*x2 - 1.7719728*x1)/(2.1055*x1 + 4.0456*x2) + 0.848*log(
1.52337552625369*x1 + 3.236*x2) + 2.16*log(1.52337552625369*x1 + 3.236*x2)
+ 0.228*log(1.52337552625369*x1 + 3.236*x2) + (2.744128*x1 + 2.744128*x2)
/(1.52337552625369*x1 + 3.236*x2) + 6.98976*x2/(1.52337552625369*x1 +
3.236*x2) + 0.737808*x2/(1.52337552625369*x1 + 3.236*x2) +
0.222260408150491*x1/(1.17581829697036*x1 + 0.197740576646344*x2) +
2735.58621973158/(226.276 + x3) - x5 =L= 11.2003192377536;
e6.. x1 + x2 =E= 1;
* set non default bounds
x1.lo = 1E-6; x1.up = 1;
x2.lo = 1E-6; x2.up = 1;
x3.lo = 40; x3.up = 90;
* set non default levels
x1.l = 0.371;
x2.l = 0.629;
x3.l = 60.632;
* 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;