s395.gms
* NLP written by GAMS Convert at 07/02/05 18:32:35
*
* Equation counts
* Total E G L N X C
* 2 2 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
* 101 1 100 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,e2;
e1.. sqr(x1) + sqr(x2) + sqr(x3) + sqr(x4) + sqr(x5) + sqr(x6) + sqr(x7) + sqr(
x8) + sqr(x9) + sqr(x10) + sqr(x11) + sqr(x12) + sqr(x13) + sqr(x14) +
sqr(x15) + sqr(x16) + sqr(x17) + sqr(x18) + sqr(x19) + sqr(x20) + sqr(x21)
+ sqr(x22) + sqr(x23) + sqr(x24) + sqr(x25) + sqr(x26) + sqr(x27) + sqr(
x28) + sqr(x29) + sqr(x30) + sqr(x31) + sqr(x32) + sqr(x33) + sqr(x34) +
sqr(x35) + sqr(x36) + sqr(x37) + sqr(x38) + sqr(x39) + sqr(x40) + sqr(x41)
+ sqr(x42) + sqr(x43) + sqr(x44) + sqr(x45) + sqr(x46) + sqr(x47) + sqr(
x48) + sqr(x49) + sqr(x50) =E= 1;
e2.. - (sqr(x1) + POWER(x1,4) + 2*(sqr(x2) + POWER(x2,4)) + 3*(sqr(x3) +
POWER(x3,4)) + 4*(sqr(x4) + POWER(x4,4)) + 5*(sqr(x5) + POWER(x5,4)) + 6*(
sqr(x6) + POWER(x6,4)) + 7*(sqr(x7) + POWER(x7,4)) + 8*(sqr(x8) + POWER(x8
,4)) + 9*(sqr(x9) + POWER(x9,4)) + 10*(sqr(x10) + POWER(x10,4)) + 11*(sqr(
x11) + POWER(x11,4)) + 12*(sqr(x12) + POWER(x12,4)) + 13*(sqr(x13) +
POWER(x13,4)) + 14*(sqr(x14) + POWER(x14,4)) + 15*(sqr(x15) + POWER(x15,4)
) + 16*(sqr(x16) + POWER(x16,4)) + 17*(sqr(x17) + POWER(x17,4)) + 18*(sqr(
x18) + POWER(x18,4)) + 19*(sqr(x19) + POWER(x19,4)) + 20*(sqr(x20) +
POWER(x20,4)) + 21*(sqr(x21) + POWER(x21,4)) + 22*(sqr(x22) + POWER(x22,4)
) + 23*(sqr(x23) + POWER(x23,4)) + 24*(sqr(x24) + POWER(x24,4)) + 25*(sqr(
x25) + POWER(x25,4)) + 26*(sqr(x26) + POWER(x26,4)) + 27*(sqr(x27) +
POWER(x27,4)) + 28*(sqr(x28) + POWER(x28,4)) + 29*(sqr(x29) + POWER(x29,4)
) + 30*(sqr(x30) + POWER(x30,4)) + 31*(sqr(x31) + POWER(x31,4)) + 32*(sqr(
x32) + POWER(x32,4)) + 33*(sqr(x33) + POWER(x33,4)) + 34*(sqr(x34) +
POWER(x34,4)) + 35*(sqr(x35) + POWER(x35,4)) + 36*(sqr(x36) + POWER(x36,4)
) + 37*(sqr(x37) + POWER(x37,4)) + 38*(sqr(x38) + POWER(x38,4)) + 39*(sqr(
x39) + POWER(x39,4)) + 40*(sqr(x40) + POWER(x40,4)) + 41*(sqr(x41) +
POWER(x41,4)) + 42*(sqr(x42) + POWER(x42,4)) + 43*(sqr(x43) + POWER(x43,4)
) + 44*(sqr(x44) + POWER(x44,4)) + 45*(sqr(x45) + POWER(x45,4)) + 46*(sqr(
x46) + POWER(x46,4)) + 47*(sqr(x47) + POWER(x47,4)) + 48*(sqr(x48) +
POWER(x48,4)) + 49*(sqr(x49) + POWER(x49,4)) + 50*(sqr(x50) + POWER(x50,4)
)) + objvar =E= 0;
* set non default bounds
* set non default levels
x1.l = 2;
x2.l = 2;
x3.l = 2;
x4.l = 2;
x5.l = 2;
x6.l = 2;
x7.l = 2;
x8.l = 2;
x9.l = 2;
x10.l = 2;
x11.l = 2;
x12.l = 2;
x13.l = 2;
x14.l = 2;
x15.l = 2;
x16.l = 2;
x17.l = 2;
x18.l = 2;
x19.l = 2;
x20.l = 2;
x21.l = 2;
x22.l = 2;
x23.l = 2;
x24.l = 2;
x25.l = 2;
x26.l = 2;
x27.l = 2;
x28.l = 2;
x29.l = 2;
x30.l = 2;
x31.l = 2;
x32.l = 2;
x33.l = 2;
x34.l = 2;
x35.l = 2;
x36.l = 2;
x37.l = 2;
x38.l = 2;
x39.l = 2;
x40.l = 2;
x41.l = 2;
x42.l = 2;
x43.l = 2;
x44.l = 2;
x45.l = 2;
x46.l = 2;
x47.l = 2;
x48.l = 2;
x49.l = 2;
x50.l = 2;
* 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;