GLOBAL World [ GLOBAL World Home | GLOBALLib | Contact ]

prolog.gms:

Reference:

Point:


* NLP written by GAMS Convert at 07/30/01 17:04:26 * * Equation counts * Total E G L N X * 23 1 0 22 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 21 21 0 0 0 0 0 0 * FX 0 0 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 129 115 14 0 * * Solve m using NLP minimizing objvar; Variables objvar,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18 ,x19,x20,x21; Positive Variables x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20 ,x21; Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19 ,e20,e21,e22,e23; e1.. x5 + x6 - 0.94*x11 - 0.94*x12 - 0.94*x13 + 0.244*x17 + 0.244*x18 + 0.244*x19 =L= 0; e2.. 0.064*x11 + 0.064*x12 + 0.064*x13 - 0.58*x14 - 0.58*x15 - 0.58*x16 + 0.172*x17 + 0.172*x18 + 0.172*x19 =L= 0; e3.. x7 + x8 + 0.048*x11 + 0.048*x12 + 0.048*x13 + 0.247*x14 + 0.247*x15 + 0.247*x16 - 0.916*x17 - 0.916*x18 - 0.916*x19 =L= 0; e4.. x11 + 1.2*x12 + 0.8*x13 + 2*x14 + 1.8*x15 + 2.4*x16 + 3*x17 + 2.7*x18 + 3.2*x19 =L= 3712; e5.. 2*x11 + 1.8*x12 + 2.2*x13 + 3*x14 + 3.5*x15 + 2.3*x16 + 3*x17 + 3.2*x18 + 2.7*x19 =L= 5000; e6.. 356.474947137148*x2 + 53.7083537310174*x4 + x5 - 0.564264890180399*x20 =L= 352; e7.. 339.983422262764*x2 + 43.5418249774113*x4 + x6 - 0.405939876920766*x21 =L= 430; e8.. 106.946746119538*x2 + 145.018955433089*x4 + x7 - 0.507117039797071*x20 =L= 222; e9.. 173.929713444361*x2 + 203.031384299627*x4 + x8 - 0.578889145413521*x21 =L= 292; e10.. x5*x2 + x7*x4 - x20 =L= 0; e11.. x6*x2 + x8*x4 - x21 =L= 0; e12.. - 3340.8*x9 - 500*x10 + x20 =L= 0; e13.. - 371.2*x9 - 4500*x10 + x21 =L= 0; e14.. 0.94*x2 - 0.064*x3 - 0.048*x4 - x9 - 2*x10 =L= 0; e15.. 0.94*x2 - 0.064*x3 - 0.048*x4 - 1.2*x9 - 1.8*x10 =L= 0; e16.. 0.94*x2 - 0.064*x3 - 0.048*x4 - 0.8*x9 - 2.2*x10 =L= 0; e17.. 0.58*x3 - 0.247*x4 - 2*x9 - 3*x10 =L= 0; e18.. 0.58*x3 - 0.247*x4 - 1.8*x9 - 3.5*x10 =L= 0; e19.. 0.58*x3 - 0.247*x4 - 2.4*x9 - 2.3*x10 =L= 0; e20.. - 0.244*x2 - 0.172*x3 + 0.916*x4 - 3*x9 - 3*x10 =L= 0; e21.. - 0.244*x2 - 0.172*x3 + 0.916*x4 - 2.7*x9 - 3.2*x10 =L= 0; e22.. - 0.244*x2 - 0.172*x3 + 0.916*x4 - 3.2*x9 - 2.7*x10 =L= 0; e23.. - (x5*x2 + x6*x2 + x7*x4 + x8*x4) - objvar + 3712*x9 + 5000*x10 =E= 0; * set non default bounds x2.lo = 0.2; x3.lo = 0.2; x4.lo = 0.2; * set non default levels x2.l = 0.5942; x3.l = 1.6167; x4.l = 1.31077; x5.l = 352; x6.l = 430; x7.l = 222; x8.l = 292; x9.l = 0.130670360422406; x10.l = 0.130670360422406; x20.l = 500.14934; x21.l = 638.25084; * 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;