eigena2.gms
* NLP written by GAMS Convert at 10/06/06 11:32:00
*
* Equation counts
* Total E G L N X C
* 56 56 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 111 111 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 1111 1 1110 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,x51,x52,x53
,x54,x55,x56,x57,x58,x59,x60,x61,x62,x63,x64,x65,x66,x67,x68,x69,x70
,x71,x72,x73,x74,x75,x76,x77,x78,x79,x80,x81,x82,x83,x84,x85,x86,x87
,x88,x89,x90,x91,x92,x93,x94,x95,x96,x97,x98,x99,x100,x101,x102,x103
,x104,x105,x106,x107,x108,x109,x110,objvar;
Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19
,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36
,e37,e38,e39,e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,e50,e51,e52,e53
,e54,e55,e56;
e1.. x12*x11 + x22*x21 + x32*x31 + x42*x41 + x52*x51 + x62*x61 + x72*x71 + x82*
x81 + x92*x91 + x102*x101 =E= 0;
e2.. x13*x11 + x23*x21 + x33*x31 + x43*x41 + x53*x51 + x63*x61 + x73*x71 + x83*
x81 + x93*x91 + x103*x101 =E= 0;
e3.. x13*x12 + x23*x22 + x33*x32 + x43*x42 + x53*x52 + x63*x62 + x73*x72 + x83*
x82 + x93*x92 + x103*x102 =E= 0;
e4.. x14*x11 + x24*x21 + x34*x31 + x44*x41 + x54*x51 + x64*x61 + x74*x71 + x84*
x81 + x94*x91 + x104*x101 =E= 0;
e5.. x14*x12 + x24*x22 + x34*x32 + x44*x42 + x54*x52 + x64*x62 + x74*x72 + x84*
x82 + x94*x92 + x104*x102 =E= 0;
e6.. x14*x13 + x24*x23 + x34*x33 + x44*x43 + x54*x53 + x64*x63 + x74*x73 + x84*
x83 + x94*x93 + x104*x103 =E= 0;
e7.. x15*x11 + x25*x21 + x35*x31 + x45*x41 + x55*x51 + x65*x61 + x75*x71 + x85*
x81 + x95*x91 + x105*x101 =E= 0;
e8.. x15*x12 + x25*x22 + x35*x32 + x45*x42 + x55*x52 + x65*x62 + x75*x72 + x85*
x82 + x95*x92 + x105*x102 =E= 0;
e9.. x15*x13 + x25*x23 + x35*x33 + x45*x43 + x55*x53 + x65*x63 + x75*x73 + x85*
x83 + x95*x93 + x105*x103 =E= 0;
e10.. x15*x14 + x25*x24 + x35*x34 + x45*x44 + x55*x54 + x65*x64 + x75*x74 + x85
*x84 + x95*x94 + x105*x104 =E= 0;
e11.. x16*x11 + x26*x21 + x36*x31 + x46*x41 + x56*x51 + x66*x61 + x76*x71 + x86
*x81 + x96*x91 + x106*x101 =E= 0;
e12.. x16*x12 + x26*x22 + x36*x32 + x46*x42 + x56*x52 + x66*x62 + x76*x72 + x86
*x82 + x96*x92 + x106*x102 =E= 0;
e13.. x16*x13 + x26*x23 + x36*x33 + x46*x43 + x56*x53 + x66*x63 + x76*x73 + x86
*x83 + x96*x93 + x106*x103 =E= 0;
e14.. x16*x14 + x26*x24 + x36*x34 + x46*x44 + x56*x54 + x66*x64 + x76*x74 + x86
*x84 + x96*x94 + x106*x104 =E= 0;
e15.. x16*x15 + x26*x25 + x36*x35 + x46*x45 + x56*x55 + x66*x65 + x76*x75 + x86
*x85 + x96*x95 + x106*x105 =E= 0;
e16.. x17*x11 + x27*x21 + x37*x31 + x47*x41 + x57*x51 + x67*x61 + x77*x71 + x87
*x81 + x97*x91 + x107*x101 =E= 0;
e17.. x17*x12 + x27*x22 + x37*x32 + x47*x42 + x57*x52 + x67*x62 + x77*x72 + x87
*x82 + x97*x92 + x107*x102 =E= 0;
e18.. x17*x13 + x27*x23 + x37*x33 + x47*x43 + x57*x53 + x67*x63 + x77*x73 + x87
*x83 + x97*x93 + x107*x103 =E= 0;
e19.. x17*x14 + x27*x24 + x37*x34 + x47*x44 + x57*x54 + x67*x64 + x77*x74 + x87
*x84 + x97*x94 + x107*x104 =E= 0;
e20.. x17*x15 + x27*x25 + x37*x35 + x47*x45 + x57*x55 + x67*x65 + x77*x75 + x87
*x85 + x97*x95 + x107*x105 =E= 0;
e21.. x17*x16 + x27*x26 + x37*x36 + x47*x46 + x57*x56 + x67*x66 + x77*x76 + x87
*x86 + x97*x96 + x107*x106 =E= 0;
e22.. x18*x11 + x28*x21 + x38*x31 + x48*x41 + x58*x51 + x68*x61 + x78*x71 + x88
*x81 + x98*x91 + x108*x101 =E= 0;
e23.. x18*x12 + x28*x22 + x38*x32 + x48*x42 + x58*x52 + x68*x62 + x78*x72 + x88
*x82 + x98*x92 + x108*x102 =E= 0;
e24.. x18*x13 + x28*x23 + x38*x33 + x48*x43 + x58*x53 + x68*x63 + x78*x73 + x88
*x83 + x98*x93 + x108*x103 =E= 0;
e25.. x18*x14 + x28*x24 + x38*x34 + x48*x44 + x58*x54 + x68*x64 + x78*x74 + x88
*x84 + x98*x94 + x108*x104 =E= 0;
e26.. x18*x15 + x28*x25 + x38*x35 + x48*x45 + x58*x55 + x68*x65 + x78*x75 + x88
*x85 + x98*x95 + x108*x105 =E= 0;
e27.. x18*x16 + x28*x26 + x38*x36 + x48*x46 + x58*x56 + x68*x66 + x78*x76 + x88
*x86 + x98*x96 + x108*x106 =E= 0;
e28.. x18*x17 + x28*x27 + x38*x37 + x48*x47 + x58*x57 + x68*x67 + x78*x77 + x88
*x87 + x98*x97 + x108*x107 =E= 0;
e29.. x19*x11 + x29*x21 + x39*x31 + x49*x41 + x59*x51 + x69*x61 + x79*x71 + x89
*x81 + x99*x91 + x109*x101 =E= 0;
e30.. x19*x12 + x29*x22 + x39*x32 + x49*x42 + x59*x52 + x69*x62 + x79*x72 + x89
*x82 + x99*x92 + x109*x102 =E= 0;
e31.. x19*x13 + x29*x23 + x39*x33 + x49*x43 + x59*x53 + x69*x63 + x79*x73 + x89
*x83 + x99*x93 + x109*x103 =E= 0;
e32.. x19*x14 + x29*x24 + x39*x34 + x49*x44 + x59*x54 + x69*x64 + x79*x74 + x89
*x84 + x99*x94 + x109*x104 =E= 0;
e33.. x19*x15 + x29*x25 + x39*x35 + x49*x45 + x59*x55 + x69*x65 + x79*x75 + x89
*x85 + x99*x95 + x109*x105 =E= 0;
e34.. x19*x16 + x29*x26 + x39*x36 + x49*x46 + x59*x56 + x69*x66 + x79*x76 + x89
*x86 + x99*x96 + x109*x106 =E= 0;
e35.. x19*x17 + x29*x27 + x39*x37 + x49*x47 + x59*x57 + x69*x67 + x79*x77 + x89
*x87 + x99*x97 + x109*x107 =E= 0;
e36.. x19*x18 + x29*x28 + x39*x38 + x49*x48 + x59*x58 + x69*x68 + x79*x78 + x89
*x88 + x99*x98 + x109*x108 =E= 0;
e37.. x20*x11 + x30*x21 + x40*x31 + x50*x41 + x60*x51 + x70*x61 + x80*x71 + x90
*x81 + x100*x91 + x110*x101 =E= 0;
e38.. x20*x12 + x30*x22 + x40*x32 + x50*x42 + x60*x52 + x70*x62 + x80*x72 + x90
*x82 + x100*x92 + x110*x102 =E= 0;
e39.. x20*x13 + x30*x23 + x40*x33 + x50*x43 + x60*x53 + x70*x63 + x80*x73 + x90
*x83 + x100*x93 + x110*x103 =E= 0;
e40.. x20*x14 + x30*x24 + x40*x34 + x50*x44 + x60*x54 + x70*x64 + x80*x74 + x90
*x84 + x100*x94 + x110*x104 =E= 0;
e41.. x20*x15 + x30*x25 + x40*x35 + x50*x45 + x60*x55 + x70*x65 + x80*x75 + x90
*x85 + x100*x95 + x110*x105 =E= 0;
e42.. x20*x16 + x30*x26 + x40*x36 + x50*x46 + x60*x56 + x70*x66 + x80*x76 + x90
*x86 + x100*x96 + x110*x106 =E= 0;
e43.. x20*x17 + x30*x27 + x40*x37 + x50*x47 + x60*x57 + x70*x67 + x80*x77 + x90
*x87 + x100*x97 + x110*x107 =E= 0;
e44.. x20*x18 + x30*x28 + x40*x38 + x50*x48 + x60*x58 + x70*x68 + x80*x78 + x90
*x88 + x100*x98 + x110*x108 =E= 0;
e45.. x20*x19 + x30*x29 + x40*x39 + x50*x49 + x60*x59 + x70*x69 + x80*x79 + x90
*x89 + x100*x99 + x110*x109 =E= 0;
e46.. x11*x11 + x21*x21 + x31*x31 + x41*x41 + x51*x51 + x61*x61 + x71*x71 + x81
*x81 + x91*x91 + x101*x101 =E= 1;
e47.. x12*x12 + x22*x22 + x32*x32 + x42*x42 + x52*x52 + x62*x62 + x72*x72 + x82
*x82 + x92*x92 + x102*x102 =E= 1;
e48.. x13*x13 + x23*x23 + x33*x33 + x43*x43 + x53*x53 + x63*x63 + x73*x73 + x83
*x83 + x93*x93 + x103*x103 =E= 1;
e49.. x14*x14 + x24*x24 + x34*x34 + x44*x44 + x54*x54 + x64*x64 + x74*x74 + x84
*x84 + x94*x94 + x104*x104 =E= 1;
e50.. x15*x15 + x25*x25 + x35*x35 + x45*x45 + x55*x55 + x65*x65 + x75*x75 + x85
*x85 + x95*x95 + x105*x105 =E= 1;
e51.. x16*x16 + x26*x26 + x36*x36 + x46*x46 + x56*x56 + x66*x66 + x76*x76 + x86
*x86 + x96*x96 + x106*x106 =E= 1;
e52.. x17*x17 + x27*x27 + x37*x37 + x47*x47 + x57*x57 + x67*x67 + x77*x77 + x87
*x87 + x97*x97 + x107*x107 =E= 1;
e53.. x18*x18 + x28*x28 + x38*x38 + x48*x48 + x58*x58 + x68*x68 + x78*x78 + x88
*x88 + x98*x98 + x108*x108 =E= 1;
e54.. x19*x19 + x29*x29 + x39*x39 + x49*x49 + x59*x59 + x69*x69 + x79*x79 + x89
*x89 + x99*x99 + x109*x109 =E= 1;
e55.. x20*x20 + x30*x30 + x40*x40 + x50*x50 + x60*x60 + x70*x70 + x80*x80 + x90
*x90 + x100*x100 + x110*x110 =E= 1;
e56.. - (sqr(x11*x1 - x11 + x11*x1 + x11*x1 + x11*x1 + x11*x1 + x11*x1 + x11*
x1 + x11*x1 + x11*x1 + x11*x1) + sqr(x21*x2 - x21 + x21*x2 + x21*x2 + x21
*x2 + x21*x2 + x21*x2 + x21*x2 + x21*x2 + x21*x2 + x21*x2) + sqr(x31*x3
- x31 + x31*x3 + x31*x3 + x31*x3 + x31*x3 + x31*x3 + x31*x3 + x31*x3 +
x31*x3 + x31*x3) + sqr(x41*x4 - x41 + x41*x4 + x41*x4 + x41*x4 + x41*x4
+ x41*x4 + x41*x4 + x41*x4 + x41*x4 + x41*x4) + sqr(x51*x5 - x51 + x51*
x5 + x51*x5 + x51*x5 + x51*x5 + x51*x5 + x51*x5 + x51*x5 + x51*x5 + x51*
x5) + sqr(x61*x6 - x61 + x61*x6 + x61*x6 + x61*x6 + x61*x6 + x61*x6 + x61
*x6 + x61*x6 + x61*x6 + x61*x6) + sqr(x71*x7 - x71 + x71*x7 + x71*x7 +
x71*x7 + x71*x7 + x71*x7 + x71*x7 + x71*x7 + x71*x7 + x71*x7) + sqr(x81*
x8 - x81 + x81*x8 + x81*x8 + x81*x8 + x81*x8 + x81*x8 + x81*x8 + x81*x8
+ x81*x8 + x81*x8) + sqr(x91*x9 - x91 + x91*x9 + x91*x9 + x91*x9 + x91*
x9 + x91*x9 + x91*x9 + x91*x9 + x91*x9 + x91*x9) + sqr(x101*x10 - x101 +
x101*x10 + x101*x10 + x101*x10 + x101*x10 + x101*x10 + x101*x10 + x101*
x10 + x101*x10 + x101*x10) + sqr(x12*x1 + x12*x1 - 2*x12 + x12*x1 + x12*
x1 + x12*x1 + x12*x1 + x12*x1 + x12*x1 + x12*x1 + x12*x1) + sqr(x22*x2 +
x22*x2 - 2*x22 + x22*x2 + x22*x2 + x22*x2 + x22*x2 + x22*x2 + x22*x2 +
x22*x2 + x22*x2) + sqr(x32*x3 + x32*x3 - 2*x32 + x32*x3 + x32*x3 + x32*x3
+ x32*x3 + x32*x3 + x32*x3 + x32*x3 + x32*x3) + sqr(x42*x4 + x42*x4 - 2*
x42 + x42*x4 + x42*x4 + x42*x4 + x42*x4 + x42*x4 + x42*x4 + x42*x4 + x42*
x4) + sqr(x52*x5 + x52*x5 - 2*x52 + x52*x5 + x52*x5 + x52*x5 + x52*x5 +
x52*x5 + x52*x5 + x52*x5 + x52*x5) + sqr(x62*x6 + x62*x6 - 2*x62 + x62*x6
+ x62*x6 + x62*x6 + x62*x6 + x62*x6 + x62*x6 + x62*x6 + x62*x6) + sqr(
x72*x7 + x72*x7 - 2*x72 + x72*x7 + x72*x7 + x72*x7 + x72*x7 + x72*x7 +
x72*x7 + x72*x7 + x72*x7) + sqr(x82*x8 + x82*x8 - 2*x82 + x82*x8 + x82*x8
+ x82*x8 + x82*x8 + x82*x8 + x82*x8 + x82*x8 + x82*x8) + sqr(x92*x9 +
x92*x9 - 2*x92 + x92*x9 + x92*x9 + x92*x9 + x92*x9 + x92*x9 + x92*x9 +
x92*x9 + x92*x9) + sqr(x102*x10 + x102*x10 - 2*x102 + x102*x10 + x102*x10
+ x102*x10 + x102*x10 + x102*x10 + x102*x10 + x102*x10 + x102*x10) +
sqr(x13*x1 + x13*x1 + x13*x1 - 3*x13 + x13*x1 + x13*x1 + x13*x1 + x13*x1
+ x13*x1 + x13*x1 + x13*x1) + sqr(x23*x2 + x23*x2 + x23*x2 - 3*x23 + x23
*x2 + x23*x2 + x23*x2 + x23*x2 + x23*x2 + x23*x2 + x23*x2) + sqr(x33*x3
+ x33*x3 + x33*x3 - 3*x33 + x33*x3 + x33*x3 + x33*x3 + x33*x3 + x33*x3
+ x33*x3 + x33*x3) + sqr(x43*x4 + x43*x4 + x43*x4 - 3*x43 + x43*x4 + x43
*x4 + x43*x4 + x43*x4 + x43*x4 + x43*x4 + x43*x4) + sqr(x53*x5 + x53*x5
+ x53*x5 - 3*x53 + x53*x5 + x53*x5 + x53*x5 + x53*x5 + x53*x5 + x53*x5
+ x53*x5) + sqr(x63*x6 + x63*x6 + x63*x6 - 3*x63 + x63*x6 + x63*x6 + x63
*x6 + x63*x6 + x63*x6 + x63*x6 + x63*x6) + sqr(x73*x7 + x73*x7 + x73*x7
- 3*x73 + x73*x7 + x73*x7 + x73*x7 + x73*x7 + x73*x7 + x73*x7 + x73*x7)
+ sqr(x83*x8 + x83*x8 + x83*x8 - 3*x83 + x83*x8 + x83*x8 + x83*x8 + x83*
x8 + x83*x8 + x83*x8 + x83*x8) + sqr(x93*x9 + x93*x9 + x93*x9 - 3*x93 +
x93*x9 + x93*x9 + x93*x9 + x93*x9 + x93*x9 + x93*x9 + x93*x9) + sqr(x103*
x10 + x103*x10 + x103*x10 - 3*x103 + x103*x10 + x103*x10 + x103*x10 +
x103*x10 + x103*x10 + x103*x10 + x103*x10) + sqr(x14*x1 + x14*x1 + x14*x1
+ x14*x1 - 4*x14 + x14*x1 + x14*x1 + x14*x1 + x14*x1 + x14*x1 + x14*x1)
+ sqr(x24*x2 + x24*x2 + x24*x2 + x24*x2 - 4*x24 + x24*x2 + x24*x2 + x24*
x2 + x24*x2 + x24*x2 + x24*x2) + sqr(x34*x3 + x34*x3 + x34*x3 + x34*x3 -
4*x34 + x34*x3 + x34*x3 + x34*x3 + x34*x3 + x34*x3 + x34*x3) + sqr(x44*x4
+ x44*x4 + x44*x4 + x44*x4 - 4*x44 + x44*x4 + x44*x4 + x44*x4 + x44*x4
+ x44*x4 + x44*x4) + sqr(x54*x5 + x54*x5 + x54*x5 + x54*x5 - 4*x54 + x54
*x5 + x54*x5 + x54*x5 + x54*x5 + x54*x5 + x54*x5) + sqr(x64*x6 + x64*x6
+ x64*x6 + x64*x6 - 4*x64 + x64*x6 + x64*x6 + x64*x6 + x64*x6 + x64*x6
+ x64*x6) + sqr(x74*x7 + x74*x7 + x74*x7 + x74*x7 - 4*x74 + x74*x7 + x74
*x7 + x74*x7 + x74*x7 + x74*x7 + x74*x7) + sqr(x84*x8 + x84*x8 + x84*x8
+ x84*x8 - 4*x84 + x84*x8 + x84*x8 + x84*x8 + x84*x8 + x84*x8 + x84*x8)
+ sqr(x94*x9 + x94*x9 + x94*x9 + x94*x9 - 4*x94 + x94*x9 + x94*x9 + x94*
x9 + x94*x9 + x94*x9 + x94*x9) + sqr(x104*x10 + x104*x10 + x104*x10 +
x104*x10 - 4*x104 + x104*x10 + x104*x10 + x104*x10 + x104*x10 + x104*x10
+ x104*x10) + sqr(x15*x1 + x15*x1 + x15*x1 + x15*x1 + x15*x1 - 5*x15 +
x15*x1 + x15*x1 + x15*x1 + x15*x1 + x15*x1) + sqr(x25*x2 + x25*x2 + x25*
x2 + x25*x2 + x25*x2 - 5*x25 + x25*x2 + x25*x2 + x25*x2 + x25*x2 + x25*x2
) + sqr(x35*x3 + x35*x3 + x35*x3 + x35*x3 + x35*x3 - 5*x35 + x35*x3 + x35
*x3 + x35*x3 + x35*x3 + x35*x3) + sqr(x45*x4 + x45*x4 + x45*x4 + x45*x4
+ x45*x4 - 5*x45 + x45*x4 + x45*x4 + x45*x4 + x45*x4 + x45*x4) + sqr(x55
*x5 + x55*x5 + x55*x5 + x55*x5 + x55*x5 - 5*x55 + x55*x5 + x55*x5 + x55*
x5 + x55*x5 + x55*x5) + sqr(x65*x6 + x65*x6 + x65*x6 + x65*x6 + x65*x6 -
5*x65 + x65*x6 + x65*x6 + x65*x6 + x65*x6 + x65*x6) + sqr(x75*x7 + x75*x7
+ x75*x7 + x75*x7 + x75*x7 - 5*x75 + x75*x7 + x75*x7 + x75*x7 + x75*x7
+ x75*x7) + sqr(x85*x8 + x85*x8 + x85*x8 + x85*x8 + x85*x8 - 5*x85 + x85
*x8 + x85*x8 + x85*x8 + x85*x8 + x85*x8) + sqr(x95*x9 + x95*x9 + x95*x9
+ x95*x9 + x95*x9 - 5*x95 + x95*x9 + x95*x9 + x95*x9 + x95*x9 + x95*x9)
+ sqr(x105*x10 + x105*x10 + x105*x10 + x105*x10 + x105*x10 - 5*x105 +
x105*x10 + x105*x10 + x105*x10 + x105*x10 + x105*x10) + sqr(x16*x1 + x16*
x1 + x16*x1 + x16*x1 + x16*x1 + x16*x1 - 6*x16 + x16*x1 + x16*x1 + x16*x1
+ x16*x1) + sqr(x26*x2 + x26*x2 + x26*x2 + x26*x2 + x26*x2 + x26*x2 - 6*
x26 + x26*x2 + x26*x2 + x26*x2 + x26*x2) + sqr(x36*x3 + x36*x3 + x36*x3
+ x36*x3 + x36*x3 + x36*x3 - 6*x36 + x36*x3 + x36*x3 + x36*x3 + x36*x3)
+ sqr(x46*x4 + x46*x4 + x46*x4 + x46*x4 + x46*x4 + x46*x4 - 6*x46 + x46*
x4 + x46*x4 + x46*x4 + x46*x4) + sqr(x56*x5 + x56*x5 + x56*x5 + x56*x5 +
x56*x5 + x56*x5 - 6*x56 + x56*x5 + x56*x5 + x56*x5 + x56*x5) + sqr(x66*x6
+ x66*x6 + x66*x6 + x66*x6 + x66*x6 + x66*x6 - 6*x66 + x66*x6 + x66*x6
+ x66*x6 + x66*x6) + sqr(x76*x7 + x76*x7 + x76*x7 + x76*x7 + x76*x7 +
x76*x7 - 6*x76 + x76*x7 + x76*x7 + x76*x7 + x76*x7) + sqr(x86*x8 + x86*x8
+ x86*x8 + x86*x8 + x86*x8 + x86*x8 - 6*x86 + x86*x8 + x86*x8 + x86*x8
+ x86*x8) + sqr(x96*x9 + x96*x9 + x96*x9 + x96*x9 + x96*x9 + x96*x9 - 6*
x96 + x96*x9 + x96*x9 + x96*x9 + x96*x9) + sqr(x106*x10 + x106*x10 + x106
*x10 + x106*x10 + x106*x10 + x106*x10 - 6*x106 + x106*x10 + x106*x10 +
x106*x10 + x106*x10) + sqr(x17*x1 + x17*x1 + x17*x1 + x17*x1 + x17*x1 +
x17*x1 + x17*x1 - 7*x17 + x17*x1 + x17*x1 + x17*x1) + sqr(x27*x2 + x27*x2
+ x27*x2 + x27*x2 + x27*x2 + x27*x2 + x27*x2 - 7*x27 + x27*x2 + x27*x2
+ x27*x2) + sqr(x37*x3 + x37*x3 + x37*x3 + x37*x3 + x37*x3 + x37*x3 +
x37*x3 - 7*x37 + x37*x3 + x37*x3 + x37*x3) + sqr(x47*x4 + x47*x4 + x47*x4
+ x47*x4 + x47*x4 + x47*x4 + x47*x4 - 7*x47 + x47*x4 + x47*x4 + x47*x4)
+ sqr(x57*x5 + x57*x5 + x57*x5 + x57*x5 + x57*x5 + x57*x5 + x57*x5 - 7*
x57 + x57*x5 + x57*x5 + x57*x5) + sqr(x67*x6 + x67*x6 + x67*x6 + x67*x6
+ x67*x6 + x67*x6 + x67*x6 - 7*x67 + x67*x6 + x67*x6 + x67*x6) + sqr(x77
*x7 + x77*x7 + x77*x7 + x77*x7 + x77*x7 + x77*x7 + x77*x7 - 7*x77 + x77*
x7 + x77*x7 + x77*x7) + sqr(x87*x8 + x87*x8 + x87*x8 + x87*x8 + x87*x8 +
x87*x8 + x87*x8 - 7*x87 + x87*x8 + x87*x8 + x87*x8) + sqr(x97*x9 + x97*x9
+ x97*x9 + x97*x9 + x97*x9 + x97*x9 + x97*x9 - 7*x97 + x97*x9 + x97*x9
+ x97*x9) + sqr(x107*x10 + x107*x10 + x107*x10 + x107*x10 + x107*x10 +
x107*x10 + x107*x10 - 7*x107 + x107*x10 + x107*x10 + x107*x10) + sqr(x18*
x1 + x18*x1 + x18*x1 + x18*x1 + x18*x1 + x18*x1 + x18*x1 + x18*x1 - 8*x18
+ x18*x1 + x18*x1) + sqr(x28*x2 + x28*x2 + x28*x2 + x28*x2 + x28*x2 +
x28*x2 + x28*x2 + x28*x2 - 8*x28 + x28*x2 + x28*x2) + sqr(x38*x3 + x38*x3
+ x38*x3 + x38*x3 + x38*x3 + x38*x3 + x38*x3 + x38*x3 - 8*x38 + x38*x3
+ x38*x3) + sqr(x48*x4 + x48*x4 + x48*x4 + x48*x4 + x48*x4 + x48*x4 +
x48*x4 + x48*x4 - 8*x48 + x48*x4 + x48*x4) + sqr(x58*x5 + x58*x5 + x58*x5
+ x58*x5 + x58*x5 + x58*x5 + x58*x5 + x58*x5 - 8*x58 + x58*x5 + x58*x5)
+ sqr(x68*x6 + x68*x6 + x68*x6 + x68*x6 + x68*x6 + x68*x6 + x68*x6 + x68
*x6 - 8*x68 + x68*x6 + x68*x6) + sqr(x78*x7 + x78*x7 + x78*x7 + x78*x7 +
x78*x7 + x78*x7 + x78*x7 + x78*x7 - 8*x78 + x78*x7 + x78*x7) + sqr(x88*x8
+ x88*x8 + x88*x8 + x88*x8 + x88*x8 + x88*x8 + x88*x8 + x88*x8 - 8*x88
+ x88*x8 + x88*x8) + sqr(x98*x9 + x98*x9 + x98*x9 + x98*x9 + x98*x9 +
x98*x9 + x98*x9 + x98*x9 - 8*x98 + x98*x9 + x98*x9) + sqr(x108*x10 + x108
*x10 + x108*x10 + x108*x10 + x108*x10 + x108*x10 + x108*x10 + x108*x10 -
8*x108 + x108*x10 + x108*x10) + sqr(x19*x1 + x19*x1 + x19*x1 + x19*x1 +
x19*x1 + x19*x1 + x19*x1 + x19*x1 + x19*x1 - 9*x19 + x19*x1) + sqr(x29*x2
+ x29*x2 + x29*x2 + x29*x2 + x29*x2 + x29*x2 + x29*x2 + x29*x2 + x29*x2
- 9*x29 + x29*x2) + sqr(x39*x3 + x39*x3 + x39*x3 + x39*x3 + x39*x3 + x39
*x3 + x39*x3 + x39*x3 + x39*x3 - 9*x39 + x39*x3) + sqr(x49*x4 + x49*x4 +
x49*x4 + x49*x4 + x49*x4 + x49*x4 + x49*x4 + x49*x4 + x49*x4 - 9*x49 +
x49*x4) + sqr(x59*x5 + x59*x5 + x59*x5 + x59*x5 + x59*x5 + x59*x5 + x59*
x5 + x59*x5 + x59*x5 - 9*x59 + x59*x5) + sqr(x69*x6 + x69*x6 + x69*x6 +
x69*x6 + x69*x6 + x69*x6 + x69*x6 + x69*x6 + x69*x6 - 9*x69 + x69*x6) +
sqr(x79*x7 + x79*x7 + x79*x7 + x79*x7 + x79*x7 + x79*x7 + x79*x7 + x79*x7
+ x79*x7 - 9*x79 + x79*x7) + sqr(x89*x8 + x89*x8 + x89*x8 + x89*x8 + x89
*x8 + x89*x8 + x89*x8 + x89*x8 + x89*x8 - 9*x89 + x89*x8) + sqr(x99*x9 +
x99*x9 + x99*x9 + x99*x9 + x99*x9 + x99*x9 + x99*x9 + x99*x9 + x99*x9 - 9
*x99 + x99*x9) + sqr(x109*x10 + x109*x10 + x109*x10 + x109*x10 + x109*x10
+ x109*x10 + x109*x10 + x109*x10 + x109*x10 - 9*x109 + x109*x10) + sqr(
x20*x1 + x20*x1 + x20*x1 + x20*x1 + x20*x1 + x20*x1 + x20*x1 + x20*x1 +
x20*x1 + x20*x1 - 10*x20) + sqr(x30*x2 + x30*x2 + x30*x2 + x30*x2 + x30*
x2 + x30*x2 + x30*x2 + x30*x2 + x30*x2 + x30*x2 - 10*x30) + sqr(x40*x3 +
x40*x3 + x40*x3 + x40*x3 + x40*x3 + x40*x3 + x40*x3 + x40*x3 + x40*x3 +
x40*x3 - 10*x40) + sqr(x50*x4 + x50*x4 + x50*x4 + x50*x4 + x50*x4 + x50*
x4 + x50*x4 + x50*x4 + x50*x4 + x50*x4 - 10*x50) + sqr(x60*x5 + x60*x5 +
x60*x5 + x60*x5 + x60*x5 + x60*x5 + x60*x5 + x60*x5 + x60*x5 + x60*x5 -
10*x60) + sqr(x70*x6 + x70*x6 + x70*x6 + x70*x6 + x70*x6 + x70*x6 + x70*
x6 + x70*x6 + x70*x6 + x70*x6 - 10*x70) + sqr(x80*x7 + x80*x7 + x80*x7 +
x80*x7 + x80*x7 + x80*x7 + x80*x7 + x80*x7 + x80*x7 + x80*x7 - 10*x80) +
sqr(x90*x8 + x90*x8 + x90*x8 + x90*x8 + x90*x8 + x90*x8 + x90*x8 + x90*x8
+ x90*x8 + x90*x8 - 10*x90) + sqr(x100*x9 + x100*x9 + x100*x9 + x100*x9
+ x100*x9 + x100*x9 + x100*x9 + x100*x9 + x100*x9 + x100*x9 - 10*x100)
+ sqr(x110*x10 + x110*x10 + x110*x10 + x110*x10 + x110*x10 + x110*x10 +
x110*x10 + x110*x10 + x110*x10 + x110*x10 - 10*x110)) + 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;
x22.l = 1;
x33.l = 1;
x44.l = 1;
x55.l = 1;
x66.l = 1;
x77.l = 1;
x88.l = 1;
x99.l = 1;
x110.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;