chnrosnb.gms
* NLP written by GAMS Convert at 10/06/06 11:16:07
*
* Equation counts
* Total E G L N X C
* 1 1 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 51 51 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 51 1 50 0
*
* Solve m using NLP minimizing objvar;
Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19
,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36
,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48,x49,x50,objvar;
Equations e1;
e1.. - (31.36*sqr(x1 - sqr(x2)) + sqr((-1) + x2) + 92.16*sqr(x2 - sqr(x3)) +
sqr((-1) + x3) + 31.36*sqr(x3 - sqr(x4)) + sqr((-1) + x4) + 49*sqr(x4 -
sqr(x5)) + sqr((-1) + x5) + 23.04*sqr(x5 - sqr(x6)) + sqr((-1) + x6) + 81*
sqr(x6 - sqr(x7)) + sqr((-1) + x7) + 23.04*sqr(x7 - sqr(x8)) + sqr((-1) +
x8) + 16*sqr(x8 - sqr(x9)) + sqr((-1) + x9) + 19.36*sqr(x9 - sqr(x10)) +
sqr((-1) + x10) + 36*sqr(x10 - sqr(x11)) + sqr((-1) + x11) + 40.96*sqr(x11
- sqr(x12)) + sqr((-1) + x12) + 25*sqr(x12 - sqr(x13)) + sqr((-1) + x13)
+ 25*sqr(x13 - sqr(x14)) + sqr((-1) + x14) + 23.04*sqr(x14 - sqr(x15)) +
sqr((-1) + x15) + 23.04*sqr(x15 - sqr(x16)) + sqr((-1) + x16) + 31.36*sqr(
x16 - sqr(x17)) + sqr((-1) + x17) + 4*sqr(x17 - sqr(x18)) + sqr((-1) + x18
) + 4*sqr(x18 - sqr(x19)) + sqr((-1) + x19) + 25*sqr(x19 - sqr(x20)) +
sqr((-1) + x20) + 51.84*sqr(x20 - sqr(x21)) + sqr((-1) + x21) + 9*sqr(x21
- sqr(x22)) + sqr((-1) + x22) + 25*sqr(x22 - sqr(x23)) + sqr((-1) + x23)
+ 31.36*sqr(x23 - sqr(x24)) + sqr((-1) + x24) + 40.96*sqr(x24 - sqr(x25))
+ sqr((-1) + x25) + 64*sqr(x25 - sqr(x26)) + sqr((-1) + x26) + 16*sqr(x26
- sqr(x27)) + sqr((-1) + x27) + 40.96*sqr(x27 - sqr(x28)) + sqr((-1) +
x28) + 25*sqr(x28 - sqr(x29)) + sqr((-1) + x29) + 121*sqr(x29 - sqr(x30))
+ sqr((-1) + x30) + 25*sqr(x30 - sqr(x31)) + sqr((-1) + x31) + 25*sqr(x31
- sqr(x32)) + sqr((-1) + x32) + 25*sqr(x32 - sqr(x33)) + sqr((-1) + x33)
+ 144*sqr(x33 - sqr(x34)) + sqr((-1) + x34) + 36*sqr(x34 - sqr(x35)) +
sqr((-1) + x35) + 64*sqr(x35 - sqr(x36)) + sqr((-1) + x36) + 25*sqr(x36 -
sqr(x37)) + sqr((-1) + x37) + 31.36*sqr(x37 - sqr(x38)) + sqr((-1) + x38)
+ 51.84*sqr(x38 - sqr(x39)) + sqr((-1) + x39) + 36*sqr(x39 - sqr(x40)) +
sqr((-1) + x40) + 77.44*sqr(x40 - sqr(x41)) + sqr((-1) + x41) + 31.36*sqr(
x41 - sqr(x42)) + sqr((-1) + x42) + 36*sqr(x42 - sqr(x43)) + sqr((-1) +
x43) + 25*sqr(x43 - sqr(x44)) + sqr((-1) + x44) + 64*sqr(x44 - sqr(x45))
+ sqr((-1) + x45) + 36*sqr(x45 - sqr(x46)) + sqr((-1) + x46) + 25*sqr(x46
- sqr(x47)) + sqr((-1) + x47) + 31.36*sqr(x47 - sqr(x48)) + sqr((-1) +
x48) + 5.76*sqr(x48 - sqr(x49)) + sqr((-1) + x49) + 36*sqr(x49 - sqr(x50))
+ sqr((-1) + x50)) + objvar =E= 0;
* set non default bounds
* set non default levels
x1.l = -1;
x2.l = -1;
x3.l = -1;
x4.l = -1;
x5.l = -1;
x6.l = -1;
x7.l = -1;
x8.l = -1;
x9.l = -1;
x10.l = -1;
x11.l = -1;
x12.l = -1;
x13.l = -1;
x14.l = -1;
x15.l = -1;
x16.l = -1;
x17.l = -1;
x18.l = -1;
x19.l = -1;
x20.l = -1;
x21.l = -1;
x22.l = -1;
x23.l = -1;
x24.l = -1;
x25.l = -1;
x26.l = -1;
x27.l = -1;
x28.l = -1;
x29.l = -1;
x30.l = -1;
x31.l = -1;
x32.l = -1;
x33.l = -1;
x34.l = -1;
x35.l = -1;
x36.l = -1;
x37.l = -1;
x38.l = -1;
x39.l = -1;
x40.l = -1;
x41.l = -1;
x42.l = -1;
x43.l = -1;
x44.l = -1;
x45.l = -1;
x46.l = -1;
x47.l = -1;
x48.l = -1;
x49.l = -1;
x50.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;