lch.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:33:50
*
* 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
* 601 601 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 1001 1 1000 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,x111,x112,x113,x114,x115,x116
,x117,x118,x119,x120,x121,x122,x123,x124,x125,x126,x127,x128,x129
,x130,x131,x132,x133,x134,x135,x136,x137,x138,x139,x140,x141,x142
,x143,x144,x145,x146,x147,x148,x149,x150,x151,x152,x153,x154,x155
,x156,x157,x158,x159,x160,x161,x162,x163,x164,x165,x166,x167,x168
,x169,x170,x171,x172,x173,x174,x175,x176,x177,x178,x179,x180,x181
,x182,x183,x184,x185,x186,x187,x188,x189,x190,x191,x192,x193,x194
,x195,x196,x197,x198,x199,x200,x201,x202,x203,x204,x205,x206,x207
,x208,x209,x210,x211,x212,x213,x214,x215,x216,x217,x218,x219,x220
,x221,x222,x223,x224,x225,x226,x227,x228,x229,x230,x231,x232,x233
,x234,x235,x236,x237,x238,x239,x240,x241,x242,x243,x244,x245,x246
,x247,x248,x249,x250,x251,x252,x253,x254,x255,x256,x257,x258,x259
,x260,x261,x262,x263,x264,x265,x266,x267,x268,x269,x270,x271,x272
,x273,x274,x275,x276,x277,x278,x279,x280,x281,x282,x283,x284,x285
,x286,x287,x288,x289,x290,x291,x292,x293,x294,x295,x296,x297,x298
,x299,x300,x301,x302,x303,x304,x305,x306,x307,x308,x309,x310,x311
,x312,x313,x314,x315,x316,x317,x318,x319,x320,x321,x322,x323,x324
,x325,x326,x327,x328,x329,x330,x331,x332,x333,x334,x335,x336,x337
,x338,x339,x340,x341,x342,x343,x344,x345,x346,x347,x348,x349,x350
,x351,x352,x353,x354,x355,x356,x357,x358,x359,x360,x361,x362,x363
,x364,x365,x366,x367,x368,x369,x370,x371,x372,x373,x374,x375,x376
,x377,x378,x379,x380,x381,x382,x383,x384,x385,x386,x387,x388,x389
,x390,x391,x392,x393,x394,x395,x396,x397,x398,x399,x400,x401,x402
,x403,x404,x405,x406,x407,x408,x409,x410,x411,x412,x413,x414,x415
,x416,x417,x418,x419,x420,x421,x422,x423,x424,x425,x426,x427,x428
,x429,x430,x431,x432,x433,x434,x435,x436,x437,x438,x439,x440,x441
,x442,x443,x444,x445,x446,x447,x448,x449,x450,x451,x452,x453,x454
,x455,x456,x457,x458,x459,x460,x461,x462,x463,x464,x465,x466,x467
,x468,x469,x470,x471,x472,x473,x474,x475,x476,x477,x478,x479,x480
,x481,x482,x483,x484,x485,x486,x487,x488,x489,x490,x491,x492,x493
,x494,x495,x496,x497,x498,x499,x500,x501,x502,x503,x504,x505,x506
,x507,x508,x509,x510,x511,x512,x513,x514,x515,x516,x517,x518,x519
,x520,x521,x522,x523,x524,x525,x526,x527,x528,x529,x530,x531,x532
,x533,x534,x535,x536,x537,x538,x539,x540,x541,x542,x543,x544,x545
,x546,x547,x548,x549,x550,x551,x552,x553,x554,x555,x556,x557,x558
,x559,x560,x561,x562,x563,x564,x565,x566,x567,x568,x569,x570,x571
,x572,x573,x574,x575,x576,x577,x578,x579,x580,x581,x582,x583,x584
,x585,x586,x587,x588,x589,x590,x591,x592,x593,x594,x595,x596,x597
,x598,x599,x600,objvar;
Equations e1,e2;
e1.. x1*x1 + x2*x2 + x4*x4 + x5*x5 + x7*x7 + x8*x8 + x10*x10 + x11*x11 + x13*
x13 + x14*x14 + x16*x16 + x17*x17 + x19*x19 + x20*x20 + x22*x22 + x23*x23
+ x25*x25 + x26*x26 + x28*x28 + x29*x29 + x31*x31 + x32*x32 + x34*x34 +
x35*x35 + x37*x37 + x38*x38 + x40*x40 + x41*x41 + x43*x43 + x44*x44 + x46*
x46 + x47*x47 + x49*x49 + x50*x50 + x52*x52 + x53*x53 + x55*x55 + x56*x56
+ x58*x58 + x59*x59 + x61*x61 + x62*x62 + x64*x64 + x65*x65 + x67*x67 +
x68*x68 + x70*x70 + x71*x71 + x73*x73 + x74*x74 + x76*x76 + x77*x77 + x79*
x79 + x80*x80 + x82*x82 + x83*x83 + x85*x85 + x86*x86 + x88*x88 + x89*x89
+ x91*x91 + x92*x92 + x94*x94 + x95*x95 + x97*x97 + x98*x98 + x100*x100
+ x101*x101 + x103*x103 + x104*x104 + x106*x106 + x107*x107 + x109*x109
+ x110*x110 + x112*x112 + x113*x113 + x115*x115 + x116*x116 + x118*x118
+ x119*x119 + x121*x121 + x122*x122 + x124*x124 + x125*x125 + x127*x127
+ x128*x128 + x130*x130 + x131*x131 + x133*x133 + x134*x134 + x136*x136
+ x137*x137 + x139*x139 + x140*x140 + x142*x142 + x143*x143 + x145*x145
+ x146*x146 + x148*x148 + x149*x149 + x151*x151 + x152*x152 + x154*x154
+ x155*x155 + x157*x157 + x158*x158 + x160*x160 + x161*x161 + x163*x163
+ x164*x164 + x166*x166 + x167*x167 + x169*x169 + x170*x170 + x172*x172
+ x173*x173 + x175*x175 + x176*x176 + x178*x178 + x179*x179 + x181*x181
+ x182*x182 + x184*x184 + x185*x185 + x187*x187 + x188*x188 + x190*x190
+ x191*x191 + x193*x193 + x194*x194 + x196*x196 + x197*x197 + x199*x199
+ x200*x200 + x202*x202 + x203*x203 + x205*x205 + x206*x206 + x208*x208
+ x209*x209 + x211*x211 + x212*x212 + x214*x214 + x215*x215 + x217*x217
+ x218*x218 + x220*x220 + x221*x221 + x223*x223 + x224*x224 + x226*x226
+ x227*x227 + x229*x229 + x230*x230 + x232*x232 + x233*x233 + x235*x235
+ x236*x236 + x238*x238 + x239*x239 + x241*x241 + x242*x242 + x244*x244
+ x245*x245 + x247*x247 + x248*x248 + x250*x250 + x251*x251 + x253*x253
+ x254*x254 + x256*x256 + x257*x257 + x259*x259 + x260*x260 + x262*x262
+ x263*x263 + x265*x265 + x266*x266 + x268*x268 + x269*x269 + x271*x271
+ x272*x272 + x274*x274 + x275*x275 + x277*x277 + x278*x278 + x280*x280
+ x281*x281 + x283*x283 + x284*x284 + x286*x286 + x287*x287 + x289*x289
+ x290*x290 + x292*x292 + x293*x293 + x295*x295 + x296*x296 + x298*x298
+ x299*x299 + x301*x301 + x302*x302 + x304*x304 + x305*x305 + x307*x307
+ x308*x308 + x310*x310 + x311*x311 + x313*x313 + x314*x314 + x316*x316
+ x317*x317 + x319*x319 + x320*x320 + x322*x322 + x323*x323 + x325*x325
+ x326*x326 + x328*x328 + x329*x329 + x331*x331 + x332*x332 + x334*x334
+ x335*x335 + x337*x337 + x338*x338 + x340*x340 + x341*x341 + x343*x343
+ x344*x344 + x346*x346 + x347*x347 + x349*x349 + x350*x350 + x352*x352
+ x353*x353 + x355*x355 + x356*x356 + x358*x358 + x359*x359 + x361*x361
+ x362*x362 + x364*x364 + x365*x365 + x367*x367 + x368*x368 + x370*x370
+ x371*x371 + x373*x373 + x374*x374 + x376*x376 + x377*x377 + x379*x379
+ x380*x380 + x382*x382 + x383*x383 + x385*x385 + x386*x386 + x388*x388
+ x389*x389 + x391*x391 + x392*x392 + x394*x394 + x395*x395 + x397*x397
+ x398*x398 + x400*x400 + x401*x401 + x403*x403 + x404*x404 + x406*x406
+ x407*x407 + x409*x409 + x410*x410 + x412*x412 + x413*x413 + x415*x415
+ x416*x416 + x418*x418 + x419*x419 + x421*x421 + x422*x422 + x424*x424
+ x425*x425 + x427*x427 + x428*x428 + x430*x430 + x431*x431 + x433*x433
+ x434*x434 + x436*x436 + x437*x437 + x439*x439 + x440*x440 + x442*x442
+ x443*x443 + x445*x445 + x446*x446 + x448*x448 + x449*x449 + x451*x451
+ x452*x452 + x454*x454 + x455*x455 + x457*x457 + x458*x458 + x460*x460
+ x461*x461 + x463*x463 + x464*x464 + x466*x466 + x467*x467 + x469*x469
+ x470*x470 + x472*x472 + x473*x473 + x475*x475 + x476*x476 + x478*x478
+ x479*x479 + x481*x481 + x482*x482 + x484*x484 + x485*x485 + x487*x487
+ x488*x488 + x490*x490 + x491*x491 + x493*x493 + x494*x494 + x496*x496
+ x497*x497 + x499*x499 + x500*x500 + x502*x502 + x503*x503 + x505*x505
+ x506*x506 + x508*x508 + x509*x509 + x511*x511 + x512*x512 + x514*x514
+ x515*x515 + x517*x517 + x518*x518 + x520*x520 + x521*x521 + x523*x523
+ x524*x524 + x526*x526 + x527*x527 + x529*x529 + x530*x530 + x532*x532
+ x533*x533 + x535*x535 + x536*x536 + x538*x538 + x539*x539 + x541*x541
+ x542*x542 + x544*x544 + x545*x545 + x547*x547 + x548*x548 + x550*x550
+ x551*x551 + x553*x553 + x554*x554 + x556*x556 + x557*x557 + x559*x559
+ x560*x560 + x562*x562 + x563*x563 + x565*x565 + x566*x566 + x568*x568
+ x569*x569 + x571*x571 + x572*x572 + x574*x574 + x575*x575 + x577*x577
+ x578*x578 + x580*x580 + x581*x581 + x583*x583 + x584*x584 + x586*x586
+ x587*x587 + x589*x589 + x590*x590 + x592*x592 + x593*x593 + x595*x595
+ x596*x596 + x598*x598 + x599*x599 =E= 1;
e2.. - (-1.549429*x1*(x4 + x598) - 1.549429*x2*(x5 + x599) + (x1*x1 + x2*x2)*(
62*x6 - 62*x600) + (650*x600 - 650*x3)*(x600 - x3) - 1.549429*x4*(x1 + x7)
- 1.549429*x5*(x2 + x8) + (x4*x4 + x5*x5)*(62*x9 - 62*x3) + (650*x3 - 650
*x6)*(x3 - x6) - 1.549429*x7*(x4 + x10) - 1.549429*x8*(x5 + x11) + (x7*x7
+ x8*x8)*(62*x12 - 62*x6) + (650*x6 - 650*x9)*(x6 - x9) - 1.549429*x10*(
x7 + x13) - 1.549429*x11*(x8 + x14) + (x10*x10 + x11*x11)*(62*x15 - 62*x9)
+ (650*x9 - 650*x12)*(x9 - x12) - 1.549429*x13*(x10 + x16) - 1.549429*x14
*(x11 + x17) + (x13*x13 + x14*x14)*(62*x18 - 62*x12) + (650*x12 - 650*x15)
*(x12 - x15) - 1.549429*x16*(x13 + x19) - 1.549429*x17*(x14 + x20) + (x16*
x16 + x17*x17)*(62*x21 - 62*x15) + (650*x15 - 650*x18)*(x15 - x18) -
1.549429*x19*(x16 + x22) - 1.549429*x20*(x17 + x23) + (x19*x19 + x20*x20)*
(62*x24 - 62*x18) + (650*x18 - 650*x21)*(x18 - x21) - 1.549429*x22*(x19 +
x25) - 1.549429*x23*(x20 + x26) + (x22*x22 + x23*x23)*(62*x27 - 62*x21) +
(650*x21 - 650*x24)*(x21 - x24) - 1.549429*x25*(x22 + x28) - 1.549429*x26*
(x23 + x29) + (x25*x25 + x26*x26)*(62*x30 - 62*x24) + (650*x24 - 650*x27)*
(x24 - x27) - 1.549429*x28*(x25 + x31) - 1.549429*x29*(x26 + x32) + (x28*
x28 + x29*x29)*(62*x33 - 62*x27) + (650*x27 - 650*x30)*(x27 - x30) -
1.549429*x31*(x28 + x34) - 1.549429*x32*(x29 + x35) + (x31*x31 + x32*x32)*
(62*x36 - 62*x30) + (650*x30 - 650*x33)*(x30 - x33) - 1.549429*x34*(x31 +
x37) - 1.549429*x35*(x32 + x38) + (x34*x34 + x35*x35)*(62*x39 - 62*x33) +
(650*x33 - 650*x36)*(x33 - x36) - 1.549429*x37*(x34 + x40) - 1.549429*x38*
(x35 + x41) + (x37*x37 + x38*x38)*(62*x42 - 62*x36) + (650*x36 - 650*x39)*
(x36 - x39) - 1.549429*x40*(x37 + x43) - 1.549429*x41*(x38 + x44) + (x40*
x40 + x41*x41)*(62*x45 - 62*x39) + (650*x39 - 650*x42)*(x39 - x42) -
1.549429*x43*(x40 + x46) - 1.549429*x44*(x41 + x47) + (x43*x43 + x44*x44)*
(62*x48 - 62*x42) + (650*x42 - 650*x45)*(x42 - x45) - 1.549429*x46*(x43 +
x49) - 1.549429*x47*(x44 + x50) + (x46*x46 + x47*x47)*(62*x51 - 62*x45) +
(650*x45 - 650*x48)*(x45 - x48) - 1.549429*x49*(x46 + x52) - 1.549429*x50*
(x47 + x53) + (x49*x49 + x50*x50)*(62*x54 - 62*x48) + (650*x48 - 650*x51)*
(x48 - x51) - 1.549429*x52*(x49 + x55) - 1.549429*x53*(x50 + x56) + (x52*
x52 + x53*x53)*(62*x57 - 62*x51) + (650*x51 - 650*x54)*(x51 - x54) -
1.549429*x55*(x52 + x58) - 1.549429*x56*(x53 + x59) + (x55*x55 + x56*x56)*
(62*x60 - 62*x54) + (650*x54 - 650*x57)*(x54 - x57) - 1.549429*x58*(x55 +
x61) - 1.549429*x59*(x56 + x62) + (x58*x58 + x59*x59)*(62*x63 - 62*x57) +
(650*x57 - 650*x60)*(x57 - x60) - 1.549429*x61*(x58 + x64) - 1.549429*x62*
(x59 + x65) + (x61*x61 + x62*x62)*(62*x66 - 62*x60) + (650*x60 - 650*x63)*
(x60 - x63) - 1.549429*x64*(x61 + x67) - 1.549429*x65*(x62 + x68) + (x64*
x64 + x65*x65)*(62*x69 - 62*x63) + (650*x63 - 650*x66)*(x63 - x66) -
1.549429*x67*(x64 + x70) - 1.549429*x68*(x65 + x71) + (x67*x67 + x68*x68)*
(62*x72 - 62*x66) + (650*x66 - 650*x69)*(x66 - x69) - 1.549429*x70*(x67 +
x73) - 1.549429*x71*(x68 + x74) + (x70*x70 + x71*x71)*(62*x75 - 62*x69) +
(650*x69 - 650*x72)*(x69 - x72) - 1.549429*x73*(x70 + x76) - 1.549429*x74*
(x71 + x77) + (x73*x73 + x74*x74)*(62*x78 - 62*x72) + (650*x72 - 650*x75)*
(x72 - x75) - 1.549429*x76*(x73 + x79) - 1.549429*x77*(x74 + x80) + (x76*
x76 + x77*x77)*(62*x81 - 62*x75) + (650*x75 - 650*x78)*(x75 - x78) -
1.549429*x79*(x76 + x82) - 1.549429*x80*(x77 + x83) + (x79*x79 + x80*x80)*
(62*x84 - 62*x78) + (650*x78 - 650*x81)*(x78 - x81) - 1.549429*x82*(x79 +
x85) - 1.549429*x83*(x80 + x86) + (x82*x82 + x83*x83)*(62*x87 - 62*x81) +
(650*x81 - 650*x84)*(x81 - x84) - 1.549429*x85*(x82 + x88) - 1.549429*x86*
(x83 + x89) + (x85*x85 + x86*x86)*(62*x90 - 62*x84) + (650*x84 - 650*x87)*
(x84 - x87) - 1.549429*x88*(x85 + x91) - 1.549429*x89*(x86 + x92) + (x88*
x88 + x89*x89)*(62*x93 - 62*x87) + (650*x87 - 650*x90)*(x87 - x90) -
1.549429*x91*(x88 + x94) - 1.549429*x92*(x89 + x95) + (x91*x91 + x92*x92)*
(62*x96 - 62*x90) + (650*x90 - 650*x93)*(x90 - x93) - 1.549429*x94*(x91 +
x97) - 1.549429*x95*(x92 + x98) + (x94*x94 + x95*x95)*(62*x99 - 62*x93) +
(650*x93 - 650*x96)*(x93 - x96) - 1.549429*x97*(x94 + x100) - 1.549429*x98
*(x95 + x101) + (x97*x97 + x98*x98)*(62*x102 - 62*x96) + (650*x96 - 650*
x99)*(x96 - x99) - 1.549429*x100*(x97 + x103) - 1.549429*x101*(x98 + x104)
+ (x100*x100 + x101*x101)*(62*x105 - 62*x99) + (650*x99 - 650*x102)*(x99
- x102) - 1.549429*x103*(x100 + x106) - 1.549429*x104*(x101 + x107) + (
x103*x103 + x104*x104)*(62*x108 - 62*x102) + (650*x102 - 650*x105)*(x102
- x105) - 1.549429*x106*(x103 + x109) - 1.549429*x107*(x104 + x110) + (
x106*x106 + x107*x107)*(62*x111 - 62*x105) + (650*x105 - 650*x108)*(x105
- x108) - 1.549429*x109*(x106 + x112) - 1.549429*x110*(x107 + x113) + (
x109*x109 + x110*x110)*(62*x114 - 62*x108) + (650*x108 - 650*x111)*(x108
- x111) - 1.549429*x112*(x109 + x115) - 1.549429*x113*(x110 + x116) + (
x112*x112 + x113*x113)*(62*x117 - 62*x111) + (650*x111 - 650*x114)*(x111
- x114) - 1.549429*x115*(x112 + x118) - 1.549429*x116*(x113 + x119) + (
x115*x115 + x116*x116)*(62*x120 - 62*x114) + (650*x114 - 650*x117)*(x114
- x117) - 1.549429*x118*(x115 + x121) - 1.549429*x119*(x116 + x122) + (
x118*x118 + x119*x119)*(62*x123 - 62*x117) + (650*x117 - 650*x120)*(x117
- x120) - 1.549429*x121*(x118 + x124) - 1.549429*x122*(x119 + x125) + (
x121*x121 + x122*x122)*(62*x126 - 62*x120) + (650*x120 - 650*x123)*(x120
- x123) - 1.549429*x124*(x121 + x127) - 1.549429*x125*(x122 + x128) + (
x124*x124 + x125*x125)*(62*x129 - 62*x123) + (650*x123 - 650*x126)*(x123
- x126) - 1.549429*x127*(x124 + x130) - 1.549429*x128*(x125 + x131) + (
x127*x127 + x128*x128)*(62*x132 - 62*x126) + (650*x126 - 650*x129)*(x126
- x129) - 1.549429*x130*(x127 + x133) - 1.549429*x131*(x128 + x134) + (
x130*x130 + x131*x131)*(62*x135 - 62*x129) + (650*x129 - 650*x132)*(x129
- x132) - 1.549429*x133*(x130 + x136) - 1.549429*x134*(x131 + x137) + (
x133*x133 + x134*x134)*(62*x138 - 62*x132) + (650*x132 - 650*x135)*(x132
- x135) - 1.549429*x136*(x133 + x139) - 1.549429*x137*(x134 + x140) + (
x136*x136 + x137*x137)*(62*x141 - 62*x135) + (650*x135 - 650*x138)*(x135
- x138) - 1.549429*x139*(x136 + x142) - 1.549429*x140*(x137 + x143) + (
x139*x139 + x140*x140)*(62*x144 - 62*x138) + (650*x138 - 650*x141)*(x138
- x141) - 1.549429*x142*(x139 + x145) - 1.549429*x143*(x140 + x146) + (
x142*x142 + x143*x143)*(62*x147 - 62*x141) + (650*x141 - 650*x144)*(x141
- x144) - 1.549429*x145*(x142 + x148) - 1.549429*x146*(x143 + x149) + (
x145*x145 + x146*x146)*(62*x150 - 62*x144) + (650*x144 - 650*x147)*(x144
- x147) - 1.549429*x148*(x145 + x151) - 1.549429*x149*(x146 + x152) + (
x148*x148 + x149*x149)*(62*x153 - 62*x147) + (650*x147 - 650*x150)*(x147
- x150) - 1.549429*x151*(x148 + x154) - 1.549429*x152*(x149 + x155) + (
x151*x151 + x152*x152)*(62*x156 - 62*x150) + (650*x150 - 650*x153)*(x150
- x153) - 1.549429*x154*(x151 + x157) - 1.549429*x155*(x152 + x158) + (
x154*x154 + x155*x155)*(62*x159 - 62*x153) + (650*x153 - 650*x156)*(x153
- x156) - 1.549429*x157*(x154 + x160) - 1.549429*x158*(x155 + x161) + (
x157*x157 + x158*x158)*(62*x162 - 62*x156) + (650*x156 - 650*x159)*(x156
- x159) - 1.549429*x160*(x157 + x163) - 1.549429*x161*(x158 + x164) + (
x160*x160 + x161*x161)*(62*x165 - 62*x159) + (650*x159 - 650*x162)*(x159
- x162) - 1.549429*x163*(x160 + x166) - 1.549429*x164*(x161 + x167) + (
x163*x163 + x164*x164)*(62*x168 - 62*x162) + (650*x162 - 650*x165)*(x162
- x165) - 1.549429*x166*(x163 + x169) - 1.549429*x167*(x164 + x170) + (
x166*x166 + x167*x167)*(62*x171 - 62*x165) + (650*x165 - 650*x168)*(x165
- x168) - 1.549429*x169*(x166 + x172) - 1.549429*x170*(x167 + x173) + (
x169*x169 + x170*x170)*(62*x174 - 62*x168) + (650*x168 - 650*x171)*(x168
- x171) - 1.549429*x172*(x169 + x175) - 1.549429*x173*(x170 + x176) + (
x172*x172 + x173*x173)*(62*x177 - 62*x171) + (650*x171 - 650*x174)*(x171
- x174) - 1.549429*x175*(x172 + x178) - 1.549429*x176*(x173 + x179) + (
x175*x175 + x176*x176)*(62*x180 - 62*x174) + (650*x174 - 650*x177)*(x174
- x177) - 1.549429*x178*(x175 + x181) - 1.549429*x179*(x176 + x182) + (
x178*x178 + x179*x179)*(62*x183 - 62*x177) + (650*x177 - 650*x180)*(x177
- x180) - 1.549429*x181*(x178 + x184) - 1.549429*x182*(x179 + x185) + (
x181*x181 + x182*x182)*(62*x186 - 62*x180) + (650*x180 - 650*x183)*(x180
- x183) - 1.549429*x184*(x181 + x187) - 1.549429*x185*(x182 + x188) + (
x184*x184 + x185*x185)*(62*x189 - 62*x183) + (650*x183 - 650*x186)*(x183
- x186) - 1.549429*x187*(x184 + x190) - 1.549429*x188*(x185 + x191) + (
x187*x187 + x188*x188)*(62*x192 - 62*x186) + (650*x186 - 650*x189)*(x186
- x189) - 1.549429*x190*(x187 + x193) - 1.549429*x191*(x188 + x194) + (
x190*x190 + x191*x191)*(62*x195 - 62*x189) + (650*x189 - 650*x192)*(x189
- x192) - 1.549429*x193*(x190 + x196) - 1.549429*x194*(x191 + x197) + (
x193*x193 + x194*x194)*(62*x198 - 62*x192) + (650*x192 - 650*x195)*(x192
- x195) - 1.549429*x196*(x193 + x199) - 1.549429*x197*(x194 + x200) + (
x196*x196 + x197*x197)*(62*x201 - 62*x195) + (650*x195 - 650*x198)*(x195
- x198) - 1.549429*x199*(x196 + x202) - 1.549429*x200*(x197 + x203) + (
x199*x199 + x200*x200)*(62*x204 - 62*x198) + (650*x198 - 650*x201)*(x198
- x201) - 1.549429*x202*(x199 + x205) - 1.549429*x203*(x200 + x206) + (
x202*x202 + x203*x203)*(62*x207 - 62*x201) + (650*x201 - 650*x204)*(x201
- x204) - 1.549429*x205*(x202 + x208) - 1.549429*x206*(x203 + x209) + (
x205*x205 + x206*x206)*(62*x210 - 62*x204) + (650*x204 - 650*x207)*(x204
- x207) - 1.549429*x208*(x205 + x211) - 1.549429*x209*(x206 + x212) + (
x208*x208 + x209*x209)*(62*x213 - 62*x207) + (650*x207 - 650*x210)*(x207
- x210) - 1.549429*x211*(x208 + x214) - 1.549429*x212*(x209 + x215) + (
x211*x211 + x212*x212)*(62*x216 - 62*x210) + (650*x210 - 650*x213)*(x210
- x213) - 1.549429*x214*(x211 + x217) - 1.549429*x215*(x212 + x218) + (
x214*x214 + x215*x215)*(62*x219 - 62*x213) + (650*x213 - 650*x216)*(x213
- x216) - 1.549429*x217*(x214 + x220) - 1.549429*x218*(x215 + x221) + (
x217*x217 + x218*x218)*(62*x222 - 62*x216) + (650*x216 - 650*x219)*(x216
- x219) - 1.549429*x220*(x217 + x223) - 1.549429*x221*(x218 + x224) + (
x220*x220 + x221*x221)*(62*x225 - 62*x219) + (650*x219 - 650*x222)*(x219
- x222) - 1.549429*x223*(x220 + x226) - 1.549429*x224*(x221 + x227) + (
x223*x223 + x224*x224)*(62*x228 - 62*x222) + (650*x222 - 650*x225)*(x222
- x225) - 1.549429*x226*(x223 + x229) - 1.549429*x227*(x224 + x230) + (
x226*x226 + x227*x227)*(62*x231 - 62*x225) + (650*x225 - 650*x228)*(x225
- x228) - 1.549429*x229*(x226 + x232) - 1.549429*x230*(x227 + x233) + (
x229*x229 + x230*x230)*(62*x234 - 62*x228) + (650*x228 - 650*x231)*(x228
- x231) - 1.549429*x232*(x229 + x235) - 1.549429*x233*(x230 + x236) + (
x232*x232 + x233*x233)*(62*x237 - 62*x231) + (650*x231 - 650*x234)*(x231
- x234) - 1.549429*x235*(x232 + x238) - 1.549429*x236*(x233 + x239) + (
x235*x235 + x236*x236)*(62*x240 - 62*x234) + (650*x234 - 650*x237)*(x234
- x237) - 1.549429*x238*(x235 + x241) - 1.549429*x239*(x236 + x242) + (
x238*x238 + x239*x239)*(62*x243 - 62*x237) + (650*x237 - 650*x240)*(x237
- x240) - 1.549429*x241*(x238 + x244) - 1.549429*x242*(x239 + x245) + (
x241*x241 + x242*x242)*(62*x246 - 62*x240) + (650*x240 - 650*x243)*(x240
- x243) - 1.549429*x244*(x241 + x247) - 1.549429*x245*(x242 + x248) + (
x244*x244 + x245*x245)*(62*x249 - 62*x243) + (650*x243 - 650*x246)*(x243
- x246) - 1.549429*x247*(x244 + x250) - 1.549429*x248*(x245 + x251) + (
x247*x247 + x248*x248)*(62*x252 - 62*x246) + (650*x246 - 650*x249)*(x246
- x249) - 1.549429*x250*(x247 + x253) - 1.549429*x251*(x248 + x254) + (
x250*x250 + x251*x251)*(62*x255 - 62*x249) + (650*x249 - 650*x252)*(x249
- x252) - 1.549429*x253*(x250 + x256) - 1.549429*x254*(x251 + x257) + (
x253*x253 + x254*x254)*(62*x258 - 62*x252) + (650*x252 - 650*x255)*(x252
- x255) - 1.549429*x256*(x253 + x259) - 1.549429*x257*(x254 + x260) + (
x256*x256 + x257*x257)*(62*x261 - 62*x255) + (650*x255 - 650*x258)*(x255
- x258) - 1.549429*x259*(x256 + x262) - 1.549429*x260*(x257 + x263) + (
x259*x259 + x260*x260)*(62*x264 - 62*x258) + (650*x258 - 650*x261)*(x258
- x261) - 1.549429*x262*(x259 + x265) - 1.549429*x263*(x260 + x266) + (
x262*x262 + x263*x263)*(62*x267 - 62*x261) + (650*x261 - 650*x264)*(x261
- x264) - 1.549429*x265*(x262 + x268) - 1.549429*x266*(x263 + x269) + (
x265*x265 + x266*x266)*(62*x270 - 62*x264) + (650*x264 - 650*x267)*(x264
- x267) - 1.549429*x268*(x265 + x271) - 1.549429*x269*(x266 + x272) + (
x268*x268 + x269*x269)*(62*x273 - 62*x267) + (650*x267 - 650*x270)*(x267
- x270) - 1.549429*x271*(x268 + x274) - 1.549429*x272*(x269 + x275) + (
x271*x271 + x272*x272)*(62*x276 - 62*x270) + (650*x270 - 650*x273)*(x270
- x273) - 1.549429*x274*(x271 + x277) - 1.549429*x275*(x272 + x278) + (
x274*x274 + x275*x275)*(62*x279 - 62*x273) + (650*x273 - 650*x276)*(x273
- x276) - 1.549429*x277*(x274 + x280) - 1.549429*x278*(x275 + x281) + (
x277*x277 + x278*x278)*(62*x282 - 62*x276) + (650*x276 - 650*x279)*(x276
- x279) - 1.549429*x280*(x277 + x283) - 1.549429*x281*(x278 + x284) + (
x280*x280 + x281*x281)*(62*x285 - 62*x279) + (650*x279 - 650*x282)*(x279
- x282) - 1.549429*x283*(x280 + x286) - 1.549429*x284*(x281 + x287) + (
x283*x283 + x284*x284)*(62*x288 - 62*x282) + (650*x282 - 650*x285)*(x282
- x285) - 1.549429*x286*(x283 + x289) - 1.549429*x287*(x284 + x290) + (
x286*x286 + x287*x287)*(62*x291 - 62*x285) + (650*x285 - 650*x288)*(x285
- x288) - 1.549429*x289*(x286 + x292) - 1.549429*x290*(x287 + x293) + (
x289*x289 + x290*x290)*(62*x294 - 62*x288) + (650*x288 - 650*x291)*(x288
- x291) - 1.549429*x292*(x289 + x295) - 1.549429*x293*(x290 + x296) + (
x292*x292 + x293*x293)*(62*x297 - 62*x291) + (650*x291 - 650*x294)*(x291
- x294) - 1.549429*x295*(x292 + x298) - 1.549429*x296*(x293 + x299) + (
x295*x295 + x296*x296)*(62*x300 - 62*x294) + (650*x294 - 650*x297)*(x294
- x297) - 1.549429*x298*(x295 + x301) - 1.549429*x299*(x296 + x302) + (
x298*x298 + x299*x299)*(62*x303 - 62*x297) + (650*x297 - 650*x300)*(x297
- x300) - 1.549429*x301*(x298 + x304) - 1.549429*x302*(x299 + x305) + (
x301*x301 + x302*x302)*(62*x306 - 62*x300) + (650*x300 - 650*x303)*(x300
- x303) - 1.549429*x304*(x301 + x307) - 1.549429*x305*(x302 + x308) + (
x304*x304 + x305*x305)*(62*x309 - 62*x303) + (650*x303 - 650*x306)*(x303
- x306) - 1.549429*x307*(x304 + x310) - 1.549429*x308*(x305 + x311) + (
x307*x307 + x308*x308)*(62*x312 - 62*x306) + (650*x306 - 650*x309)*(x306
- x309) - 1.549429*x310*(x307 + x313) - 1.549429*x311*(x308 + x314) + (
x310*x310 + x311*x311)*(62*x315 - 62*x309) + (650*x309 - 650*x312)*(x309
- x312) - 1.549429*x313*(x310 + x316) - 1.549429*x314*(x311 + x317) + (
x313*x313 + x314*x314)*(62*x318 - 62*x312) + (650*x312 - 650*x315)*(x312
- x315) - 1.549429*x316*(x313 + x319) - 1.549429*x317*(x314 + x320) + (
x316*x316 + x317*x317)*(62*x321 - 62*x315) + (650*x315 - 650*x318)*(x315
- x318) - 1.549429*x319*(x316 + x322) - 1.549429*x320*(x317 + x323) + (
x319*x319 + x320*x320)*(62*x324 - 62*x318) + (650*x318 - 650*x321)*(x318
- x321) - 1.549429*x322*(x319 + x325) - 1.549429*x323*(x320 + x326) + (
x322*x322 + x323*x323)*(62*x327 - 62*x321) + (650*x321 - 650*x324)*(x321
- x324) - 1.549429*x325*(x322 + x328) - 1.549429*x326*(x323 + x329) + (
x325*x325 + x326*x326)*(62*x330 - 62*x324) + (650*x324 - 650*x327)*(x324
- x327) - 1.549429*x328*(x325 + x331) - 1.549429*x329*(x326 + x332) + (
x328*x328 + x329*x329)*(62*x333 - 62*x327) + (650*x327 - 650*x330)*(x327
- x330) - 1.549429*x331*(x328 + x334) - 1.549429*x332*(x329 + x335) + (
x331*x331 + x332*x332)*(62*x336 - 62*x330) + (650*x330 - 650*x333)*(x330
- x333) - 1.549429*x334*(x331 + x337) - 1.549429*x335*(x332 + x338) + (
x334*x334 + x335*x335)*(62*x339 - 62*x333) + (650*x333 - 650*x336)*(x333
- x336) - 1.549429*x337*(x334 + x340) - 1.549429*x338*(x335 + x341) + (
x337*x337 + x338*x338)*(62*x342 - 62*x336) + (650*x336 - 650*x339)*(x336
- x339) - 1.549429*x340*(x337 + x343) - 1.549429*x341*(x338 + x344) + (
x340*x340 + x341*x341)*(62*x345 - 62*x339) + (650*x339 - 650*x342)*(x339
- x342) - 1.549429*x343*(x340 + x346) - 1.549429*x344*(x341 + x347) + (
x343*x343 + x344*x344)*(62*x348 - 62*x342) + (650*x342 - 650*x345)*(x342
- x345) - 1.549429*x346*(x343 + x349) - 1.549429*x347*(x344 + x350) + (
x346*x346 + x347*x347)*(62*x351 - 62*x345) + (650*x345 - 650*x348)*(x345
- x348) - 1.549429*x349*(x346 + x352) - 1.549429*x350*(x347 + x353) + (
x349*x349 + x350*x350)*(62*x354 - 62*x348) + (650*x348 - 650*x351)*(x348
- x351) - 1.549429*x352*(x349 + x355) - 1.549429*x353*(x350 + x356) + (
x352*x352 + x353*x353)*(62*x357 - 62*x351) + (650*x351 - 650*x354)*(x351
- x354) - 1.549429*x355*(x352 + x358) - 1.549429*x356*(x353 + x359) + (
x355*x355 + x356*x356)*(62*x360 - 62*x354) + (650*x354 - 650*x357)*(x354
- x357) - 1.549429*x358*(x355 + x361) - 1.549429*x359*(x356 + x362) + (
x358*x358 + x359*x359)*(62*x363 - 62*x357) + (650*x357 - 650*x360)*(x357
- x360) - 1.549429*x361*(x358 + x364) - 1.549429*x362*(x359 + x365) + (
x361*x361 + x362*x362)*(62*x366 - 62*x360) + (650*x360 - 650*x363)*(x360
- x363) - 1.549429*x364*(x361 + x367) - 1.549429*x365*(x362 + x368) + (
x364*x364 + x365*x365)*(62*x369 - 62*x363) + (650*x363 - 650*x366)*(x363
- x366) - 1.549429*x367*(x364 + x370) - 1.549429*x368*(x365 + x371) + (
x367*x367 + x368*x368)*(62*x372 - 62*x366) + (650*x366 - 650*x369)*(x366
- x369) - 1.549429*x370*(x367 + x373) - 1.549429*x371*(x368 + x374) + (
x370*x370 + x371*x371)*(62*x375 - 62*x369) + (650*x369 - 650*x372)*(x369
- x372) - 1.549429*x373*(x370 + x376) - 1.549429*x374*(x371 + x377) + (
x373*x373 + x374*x374)*(62*x378 - 62*x372) + (650*x372 - 650*x375)*(x372
- x375) - 1.549429*x376*(x373 + x379) - 1.549429*x377*(x374 + x380) + (
x376*x376 + x377*x377)*(62*x381 - 62*x375) + (650*x375 - 650*x378)*(x375
- x378) - 1.549429*x379*(x376 + x382) - 1.549429*x380*(x377 + x383) + (
x379*x379 + x380*x380)*(62*x384 - 62*x378) + (650*x378 - 650*x381)*(x378
- x381) - 1.549429*x382*(x379 + x385) - 1.549429*x383*(x380 + x386) + (
x382*x382 + x383*x383)*(62*x387 - 62*x381) + (650*x381 - 650*x384)*(x381
- x384) - 1.549429*x385*(x382 + x388) - 1.549429*x386*(x383 + x389) + (
x385*x385 + x386*x386)*(62*x390 - 62*x384) + (650*x384 - 650*x387)*(x384
- x387) - 1.549429*x388*(x385 + x391) - 1.549429*x389*(x386 + x392) + (
x388*x388 + x389*x389)*(62*x393 - 62*x387) + (650*x387 - 650*x390)*(x387
- x390) - 1.549429*x391*(x388 + x394) - 1.549429*x392*(x389 + x395) + (
x391*x391 + x392*x392)*(62*x396 - 62*x390) + (650*x390 - 650*x393)*(x390
- x393) - 1.549429*x394*(x391 + x397) - 1.549429*x395*(x392 + x398) + (
x394*x394 + x395*x395)*(62*x399 - 62*x393) + (650*x393 - 650*x396)*(x393
- x396) - 1.549429*x397*(x394 + x400) - 1.549429*x398*(x395 + x401) + (
x397*x397 + x398*x398)*(62*x402 - 62*x396) + (650*x396 - 650*x399)*(x396
- x399) - 1.549429*x400*(x397 + x403) - 1.549429*x401*(x398 + x404) + (
x400*x400 + x401*x401)*(62*x405 - 62*x399) + (650*x399 - 650*x402)*(x399
- x402) - 1.549429*x403*(x400 + x406) - 1.549429*x404*(x401 + x407) + (
x403*x403 + x404*x404)*(62*x408 - 62*x402) + (650*x402 - 650*x405)*(x402
- x405) - 1.549429*x406*(x403 + x409) - 1.549429*x407*(x404 + x410) + (
x406*x406 + x407*x407)*(62*x411 - 62*x405) + (650*x405 - 650*x408)*(x405
- x408) - 1.549429*x409*(x406 + x412) - 1.549429*x410*(x407 + x413) + (
x409*x409 + x410*x410)*(62*x414 - 62*x408) + (650*x408 - 650*x411)*(x408
- x411) - 1.549429*x412*(x409 + x415) - 1.549429*x413*(x410 + x416) + (
x412*x412 + x413*x413)*(62*x417 - 62*x411) + (650*x411 - 650*x414)*(x411
- x414) - 1.549429*x415*(x412 + x418) - 1.549429*x416*(x413 + x419) + (
x415*x415 + x416*x416)*(62*x420 - 62*x414) + (650*x414 - 650*x417)*(x414
- x417) - 1.549429*x418*(x415 + x421) - 1.549429*x419*(x416 + x422) + (
x418*x418 + x419*x419)*(62*x423 - 62*x417) + (650*x417 - 650*x420)*(x417
- x420) - 1.549429*x421*(x418 + x424) - 1.549429*x422*(x419 + x425) + (
x421*x421 + x422*x422)*(62*x426 - 62*x420) + (650*x420 - 650*x423)*(x420
- x423) - 1.549429*x424*(x421 + x427) - 1.549429*x425*(x422 + x428) + (
x424*x424 + x425*x425)*(62*x429 - 62*x423) + (650*x423 - 650*x426)*(x423
- x426) - 1.549429*x427*(x424 + x430) - 1.549429*x428*(x425 + x431) + (
x427*x427 + x428*x428)*(62*x432 - 62*x426) + (650*x426 - 650*x429)*(x426
- x429) - 1.549429*x430*(x427 + x433) - 1.549429*x431*(x428 + x434) + (
x430*x430 + x431*x431)*(62*x435 - 62*x429) + (650*x429 - 650*x432)*(x429
- x432) - 1.549429*x433*(x430 + x436) - 1.549429*x434*(x431 + x437) + (
x433*x433 + x434*x434)*(62*x438 - 62*x432) + (650*x432 - 650*x435)*(x432
- x435) - 1.549429*x436*(x433 + x439) - 1.549429*x437*(x434 + x440) + (
x436*x436 + x437*x437)*(62*x441 - 62*x435) + (650*x435 - 650*x438)*(x435
- x438) - 1.549429*x439*(x436 + x442) - 1.549429*x440*(x437 + x443) + (
x439*x439 + x440*x440)*(62*x444 - 62*x438) + (650*x438 - 650*x441)*(x438
- x441) - 1.549429*x442*(x439 + x445) - 1.549429*x443*(x440 + x446) + (
x442*x442 + x443*x443)*(62*x447 - 62*x441) + (650*x441 - 650*x444)*(x441
- x444) - 1.549429*x445*(x442 + x448) - 1.549429*x446*(x443 + x449) + (
x445*x445 + x446*x446)*(62*x450 - 62*x444) + (650*x444 - 650*x447)*(x444
- x447) - 1.549429*x448*(x445 + x451) - 1.549429*x449*(x446 + x452) + (
x448*x448 + x449*x449)*(62*x453 - 62*x447) + (650*x447 - 650*x450)*(x447
- x450) - 1.549429*x451*(x448 + x454) - 1.549429*x452*(x449 + x455) + (
x451*x451 + x452*x452)*(62*x456 - 62*x450) + (650*x450 - 650*x453)*(x450
- x453) - 1.549429*x454*(x451 + x457) - 1.549429*x455*(x452 + x458) + (
x454*x454 + x455*x455)*(62*x459 - 62*x453) + (650*x453 - 650*x456)*(x453
- x456) - 1.549429*x457*(x454 + x460) - 1.549429*x458*(x455 + x461) + (
x457*x457 + x458*x458)*(62*x462 - 62*x456) + (650*x456 - 650*x459)*(x456
- x459) - 1.549429*x460*(x457 + x463) - 1.549429*x461*(x458 + x464) + (
x460*x460 + x461*x461)*(62*x465 - 62*x459) + (650*x459 - 650*x462)*(x459
- x462) - 1.549429*x463*(x460 + x466) - 1.549429*x464*(x461 + x467) + (
x463*x463 + x464*x464)*(62*x468 - 62*x462) + (650*x462 - 650*x465)*(x462
- x465) - 1.549429*x466*(x463 + x469) - 1.549429*x467*(x464 + x470) + (
x466*x466 + x467*x467)*(62*x471 - 62*x465) + (650*x465 - 650*x468)*(x465
- x468) - 1.549429*x469*(x466 + x472) - 1.549429*x470*(x467 + x473) + (
x469*x469 + x470*x470)*(62*x474 - 62*x468) + (650*x468 - 650*x471)*(x468
- x471) - 1.549429*x472*(x469 + x475) - 1.549429*x473*(x470 + x476) + (
x472*x472 + x473*x473)*(62*x477 - 62*x471) + (650*x471 - 650*x474)*(x471
- x474) - 1.549429*x475*(x472 + x478) - 1.549429*x476*(x473 + x479) + (
x475*x475 + x476*x476)*(62*x480 - 62*x474) + (650*x474 - 650*x477)*(x474
- x477) - 1.549429*x478*(x475 + x481) - 1.549429*x479*(x476 + x482) + (
x478*x478 + x479*x479)*(62*x483 - 62*x477) + (650*x477 - 650*x480)*(x477
- x480) - 1.549429*x481*(x478 + x484) - 1.549429*x482*(x479 + x485) + (
x481*x481 + x482*x482)*(62*x486 - 62*x480) + (650*x480 - 650*x483)*(x480
- x483) - 1.549429*x484*(x481 + x487) - 1.549429*x485*(x482 + x488) + (
x484*x484 + x485*x485)*(62*x489 - 62*x483) + (650*x483 - 650*x486)*(x483
- x486) - 1.549429*x487*(x484 + x490) - 1.549429*x488*(x485 + x491) + (
x487*x487 + x488*x488)*(62*x492 - 62*x486) + (650*x486 - 650*x489)*(x486
- x489) - 1.549429*x490*(x487 + x493) - 1.549429*x491*(x488 + x494) + (
x490*x490 + x491*x491)*(62*x495 - 62*x489) + (650*x489 - 650*x492)*(x489
- x492) - 1.549429*x493*(x490 + x496) - 1.549429*x494*(x491 + x497) + (
x493*x493 + x494*x494)*(62*x498 - 62*x492) + (650*x492 - 650*x495)*(x492
- x495) - 1.549429*x496*(x493 + x499) - 1.549429*x497*(x494 + x500) + (
x496*x496 + x497*x497)*(62*x501 - 62*x495) + (650*x495 - 650*x498)*(x495
- x498) - 1.549429*x499*(x496 + x502) - 1.549429*x500*(x497 + x503) + (
x499*x499 + x500*x500)*(62*x504 - 62*x498) + (650*x498 - 650*x501)*(x498
- x501) - 1.549429*x502*(x499 + x505) - 1.549429*x503*(x500 + x506) + (
x502*x502 + x503*x503)*(62*x507 - 62*x501) + (650*x501 - 650*x504)*(x501
- x504) - 1.549429*x505*(x502 + x508) - 1.549429*x506*(x503 + x509) + (
x505*x505 + x506*x506)*(62*x510 - 62*x504) + (650*x504 - 650*x507)*(x504
- x507) - 1.549429*x508*(x505 + x511) - 1.549429*x509*(x506 + x512) + (
x508*x508 + x509*x509)*(62*x513 - 62*x507) + (650*x507 - 650*x510)*(x507
- x510) - 1.549429*x511*(x508 + x514) - 1.549429*x512*(x509 + x515) + (
x511*x511 + x512*x512)*(62*x516 - 62*x510) + (650*x510 - 650*x513)*(x510
- x513) - 1.549429*x514*(x511 + x517) - 1.549429*x515*(x512 + x518) + (
x514*x514 + x515*x515)*(62*x519 - 62*x513) + (650*x513 - 650*x516)*(x513
- x516) - 1.549429*x517*(x514 + x520) - 1.549429*x518*(x515 + x521) + (
x517*x517 + x518*x518)*(62*x522 - 62*x516) + (650*x516 - 650*x519)*(x516
- x519) - 1.549429*x520*(x517 + x523) - 1.549429*x521*(x518 + x524) + (
x520*x520 + x521*x521)*(62*x525 - 62*x519) + (650*x519 - 650*x522)*(x519
- x522) - 1.549429*x523*(x520 + x526) - 1.549429*x524*(x521 + x527) + (
x523*x523 + x524*x524)*(62*x528 - 62*x522) + (650*x522 - 650*x525)*(x522
- x525) - 1.549429*x526*(x523 + x529) - 1.549429*x527*(x524 + x530) + (
x526*x526 + x527*x527)*(62*x531 - 62*x525) + (650*x525 - 650*x528)*(x525
- x528) - 1.549429*x529*(x526 + x532) - 1.549429*x530*(x527 + x533) + (
x529*x529 + x530*x530)*(62*x534 - 62*x528) + (650*x528 - 650*x531)*(x528
- x531) - 1.549429*x532*(x529 + x535) - 1.549429*x533*(x530 + x536) + (
x532*x532 + x533*x533)*(62*x537 - 62*x531) + (650*x531 - 650*x534)*(x531
- x534) - 1.549429*x535*(x532 + x538) - 1.549429*x536*(x533 + x539) + (
x535*x535 + x536*x536)*(62*x540 - 62*x534) + (650*x534 - 650*x537)*(x534
- x537) - 1.549429*x538*(x535 + x541) - 1.549429*x539*(x536 + x542) + (
x538*x538 + x539*x539)*(62*x543 - 62*x537) + (650*x537 - 650*x540)*(x537
- x540) - 1.549429*x541*(x538 + x544) - 1.549429*x542*(x539 + x545) + (
x541*x541 + x542*x542)*(62*x546 - 62*x540) + (650*x540 - 650*x543)*(x540
- x543) - 1.549429*x544*(x541 + x547) - 1.549429*x545*(x542 + x548) + (
x544*x544 + x545*x545)*(62*x549 - 62*x543) + (650*x543 - 650*x546)*(x543
- x546) - 1.549429*x547*(x544 + x550) - 1.549429*x548*(x545 + x551) + (
x547*x547 + x548*x548)*(62*x552 - 62*x546) + (650*x546 - 650*x549)*(x546
- x549) - 1.549429*x550*(x547 + x553) - 1.549429*x551*(x548 + x554) + (
x550*x550 + x551*x551)*(62*x555 - 62*x549) + (650*x549 - 650*x552)*(x549
- x552) - 1.549429*x553*(x550 + x556) - 1.549429*x554*(x551 + x557) + (
x553*x553 + x554*x554)*(62*x558 - 62*x552) + (650*x552 - 650*x555)*(x552
- x555) - 1.549429*x556*(x553 + x559) - 1.549429*x557*(x554 + x560) + (
x556*x556 + x557*x557)*(62*x561 - 62*x555) + (650*x555 - 650*x558)*(x555
- x558) - 1.549429*x559*(x556 + x562) - 1.549429*x560*(x557 + x563) + (
x559*x559 + x560*x560)*(62*x564 - 62*x558) + (650*x558 - 650*x561)*(x558
- x561) - 1.549429*x562*(x559 + x565) - 1.549429*x563*(x560 + x566) + (
x562*x562 + x563*x563)*(62*x567 - 62*x561) + (650*x561 - 650*x564)*(x561
- x564) - 1.549429*x565*(x562 + x568) - 1.549429*x566*(x563 + x569) + (
x565*x565 + x566*x566)*(62*x570 - 62*x564) + (650*x564 - 650*x567)*(x564
- x567) - 1.549429*x568*(x565 + x571) - 1.549429*x569*(x566 + x572) + (
x568*x568 + x569*x569)*(62*x573 - 62*x567) + (650*x567 - 650*x570)*(x567
- x570) - 1.549429*x571*(x568 + x574) - 1.549429*x572*(x569 + x575) + (
x571*x571 + x572*x572)*(62*x576 - 62*x570) + (650*x570 - 650*x573)*(x570
- x573) - 1.549429*x574*(x571 + x577) - 1.549429*x575*(x572 + x578) + (
x574*x574 + x575*x575)*(62*x579 - 62*x573) + (650*x573 - 650*x576)*(x573
- x576) - 1.549429*x577*(x574 + x580) - 1.549429*x578*(x575 + x581) + (
x577*x577 + x578*x578)*(62*x582 - 62*x576) + (650*x576 - 650*x579)*(x576
- x579) - 1.549429*x580*(x577 + x583) - 1.549429*x581*(x578 + x584) + (
x580*x580 + x581*x581)*(62*x585 - 62*x579) + (650*x579 - 650*x582)*(x579
- x582) - 1.549429*x583*(x580 + x586) - 1.549429*x584*(x581 + x587) + (
x583*x583 + x584*x584)*(62*x588 - 62*x582) + (650*x582 - 650*x585)*(x582
- x585) - 1.549429*x586*(x583 + x589) - 1.549429*x587*(x584 + x590) + (
x586*x586 + x587*x587)*(62*x591 - 62*x585) + (650*x585 - 650*x588)*(x585
- x588) - 1.549429*x589*(x586 + x592) - 1.549429*x590*(x587 + x593) + (
x589*x589 + x590*x590)*(62*x594 - 62*x588) + (650*x588 - 650*x591)*(x588
- x591) - 1.549429*x592*(x589 + x595) - 1.549429*x593*(x590 + x596) + (
x592*x592 + x593*x593)*(62*x597 - 62*x591) + (650*x591 - 650*x594)*(x591
- x594) - 1.549429*x595*(x592 + x598) - 1.549429*x596*(x593 + x599) + (
x595*x595 + x596*x596)*(62*x600 - 62*x594) + (650*x594 - 650*x597)*(x594
- x597) - 1.549429*x598*(x1 + x595) - 1.549429*x599*(x2 + x596) + (x598*
x598 + x599*x599)*(62*x3 - 62*x597) + (650*x597 - 650*x600)*(x597 - x600))
+ objvar =E= 0;
* set non default bounds
* set non default levels
x1.l = 0.05;
x2.l = 0.05;
x3.l = 0.1;
x4.l = 0.05;
x5.l = 0.05;
x6.l = 0.2;
x7.l = 0.05;
x8.l = 0.05;
x9.l = 0.3;
x10.l = 0.05;
x11.l = 0.05;
x12.l = 0.4;
x13.l = 0.05;
x14.l = 0.05;
x15.l = 0.5;
x16.l = 0.05;
x17.l = 0.05;
x18.l = 0.6;
x19.l = 0.05;
x20.l = 0.05;
x21.l = 0.7;
x22.l = 0.05;
x23.l = 0.05;
x24.l = 0.8;
x25.l = 0.05;
x26.l = 0.05;
x27.l = 0.9;
x28.l = 0.05;
x29.l = 0.05;
x30.l = 1;
x31.l = 0.05;
x32.l = 0.05;
x33.l = 1.1;
x34.l = 0.05;
x35.l = 0.05;
x36.l = 1.2;
x37.l = 0.05;
x38.l = 0.05;
x39.l = 1.3;
x40.l = 0.05;
x41.l = 0.05;
x42.l = 1.4;
x43.l = 0.05;
x44.l = 0.05;
x45.l = 1.5;
x46.l = 0.05;
x47.l = 0.05;
x48.l = 1.6;
x49.l = 0.05;
x50.l = 0.05;
x51.l = 1.7;
x52.l = 0.05;
x53.l = 0.05;
x54.l = 1.8;
x55.l = 0.05;
x56.l = 0.05;
x57.l = 1.9;
x58.l = 0.05;
x59.l = 0.05;
x60.l = 2;
x61.l = 0.05;
x62.l = 0.05;
x63.l = 2.1;
x64.l = 0.05;
x65.l = 0.05;
x66.l = 2.2;
x67.l = 0.05;
x68.l = 0.05;
x69.l = 2.3;
x70.l = 0.05;
x71.l = 0.05;
x72.l = 2.4;
x73.l = 0.05;
x74.l = 0.05;
x75.l = 2.5;
x76.l = 0.05;
x77.l = 0.05;
x78.l = 2.6;
x79.l = 0.05;
x80.l = 0.05;
x81.l = 2.7;
x82.l = 0.05;
x83.l = 0.05;
x84.l = 2.8;
x85.l = 0.05;
x86.l = 0.05;
x87.l = 2.9;
x88.l = 0.05;
x89.l = 0.05;
x90.l = 3;
x91.l = 0.05;
x92.l = 0.05;
x93.l = 3.1;
x94.l = 0.05;
x95.l = 0.05;
x96.l = 3.2;
x97.l = 0.05;
x98.l = 0.05;
x99.l = 3.3;
x100.l = 0.05;
x101.l = 0.05;
x102.l = 3.4;
x103.l = 0.05;
x104.l = 0.05;
x105.l = 3.5;
x106.l = 0.05;
x107.l = 0.05;
x108.l = 3.6;
x109.l = 0.05;
x110.l = 0.05;
x111.l = 3.7;
x112.l = 0.05;
x113.l = 0.05;
x114.l = 3.8;
x115.l = 0.05;
x116.l = 0.05;
x117.l = 3.9;
x118.l = 0.05;
x119.l = 0.05;
x120.l = 4;
x121.l = 0.05;
x122.l = 0.05;
x123.l = 4.1;
x124.l = 0.05;
x125.l = 0.05;
x126.l = 4.2;
x127.l = 0.05;
x128.l = 0.05;
x129.l = 4.3;
x130.l = 0.05;
x131.l = 0.05;
x132.l = 4.4;
x133.l = 0.05;
x134.l = 0.05;
x135.l = 4.5;
x136.l = 0.05;
x137.l = 0.05;
x138.l = 4.6;
x139.l = 0.05;
x140.l = 0.05;
x141.l = 4.7;
x142.l = 0.05;
x143.l = 0.05;
x144.l = 4.8;
x145.l = 0.05;
x146.l = 0.05;
x147.l = 4.9;
x148.l = 0.05;
x149.l = 0.05;
x150.l = 5;
x151.l = 0.05;
x152.l = 0.05;
x153.l = 5.1;
x154.l = 0.05;
x155.l = 0.05;
x156.l = 5.2;
x157.l = 0.05;
x158.l = 0.05;
x159.l = 5.3;
x160.l = 0.05;
x161.l = 0.05;
x162.l = 5.4;
x163.l = 0.05;
x164.l = 0.05;
x165.l = 5.5;
x166.l = 0.05;
x167.l = 0.05;
x168.l = 5.6;
x169.l = 0.05;
x170.l = 0.05;
x171.l = 5.7;
x172.l = 0.05;
x173.l = 0.05;
x174.l = 5.8;
x175.l = 0.05;
x176.l = 0.05;
x177.l = 5.9;
x178.l = 0.05;
x179.l = 0.05;
x180.l = 5.99999999999999;
x181.l = 0.05;
x182.l = 0.05;
x183.l = 6.09999999999999;
x184.l = 0.05;
x185.l = 0.05;
x186.l = 6.19999999999999;
x187.l = 0.05;
x188.l = 0.05;
x189.l = 6.29999999999999;
x190.l = 0.05;
x191.l = 0.05;
x192.l = 6.39999999999999;
x193.l = 0.05;
x194.l = 0.05;
x195.l = 6.49999999999999;
x196.l = 0.05;
x197.l = 0.05;
x198.l = 6.59999999999999;
x199.l = 0.05;
x200.l = 0.05;
x201.l = 6.69999999999999;
x202.l = 0.05;
x203.l = 0.05;
x204.l = 6.79999999999999;
x205.l = 0.05;
x206.l = 0.05;
x207.l = 6.89999999999999;
x208.l = 0.05;
x209.l = 0.05;
x210.l = 6.99999999999999;
x211.l = 0.05;
x212.l = 0.05;
x213.l = 7.09999999999999;
x214.l = 0.05;
x215.l = 0.05;
x216.l = 7.19999999999999;
x217.l = 0.05;
x218.l = 0.05;
x219.l = 7.29999999999999;
x220.l = 0.05;
x221.l = 0.05;
x222.l = 7.39999999999999;
x223.l = 0.05;
x224.l = 0.05;
x225.l = 7.49999999999999;
x226.l = 0.05;
x227.l = 0.05;
x228.l = 7.59999999999999;
x229.l = 0.05;
x230.l = 0.05;
x231.l = 7.69999999999999;
x232.l = 0.05;
x233.l = 0.05;
x234.l = 7.79999999999999;
x235.l = 0.05;
x236.l = 0.05;
x237.l = 7.89999999999999;
x238.l = 0.05;
x239.l = 0.05;
x240.l = 7.99999999999999;
x241.l = 0.05;
x242.l = 0.05;
x243.l = 8.09999999999999;
x244.l = 0.05;
x245.l = 0.05;
x246.l = 8.19999999999999;
x247.l = 0.05;
x248.l = 0.05;
x249.l = 8.29999999999999;
x250.l = 0.05;
x251.l = 0.05;
x252.l = 8.39999999999999;
x253.l = 0.05;
x254.l = 0.05;
x255.l = 8.49999999999999;
x256.l = 0.05;
x257.l = 0.05;
x258.l = 8.59999999999999;
x259.l = 0.05;
x260.l = 0.05;
x261.l = 8.69999999999999;
x262.l = 0.05;
x263.l = 0.05;
x264.l = 8.79999999999998;
x265.l = 0.05;
x266.l = 0.05;
x267.l = 8.89999999999998;
x268.l = 0.05;
x269.l = 0.05;
x270.l = 8.99999999999998;
x271.l = 0.05;
x272.l = 0.05;
x273.l = 9.09999999999998;
x274.l = 0.05;
x275.l = 0.05;
x276.l = 9.19999999999998;
x277.l = 0.05;
x278.l = 0.05;
x279.l = 9.29999999999998;
x280.l = 0.05;
x281.l = 0.05;
x282.l = 9.39999999999998;
x283.l = 0.05;
x284.l = 0.05;
x285.l = 9.49999999999998;
x286.l = 0.05;
x287.l = 0.05;
x288.l = 9.59999999999998;
x289.l = 0.05;
x290.l = 0.05;
x291.l = 9.69999999999998;
x292.l = 0.05;
x293.l = 0.05;
x294.l = 9.79999999999998;
x295.l = 0.05;
x296.l = 0.05;
x297.l = 9.89999999999998;
x298.l = 0.05;
x299.l = 0.05;
x300.l = 9.99999999999998;
x301.l = 0.05;
x302.l = 0.05;
x303.l = 10.1;
x304.l = 0.05;
x305.l = 0.05;
x306.l = 10.2;
x307.l = 0.05;
x308.l = 0.05;
x309.l = 10.3;
x310.l = 0.05;
x311.l = 0.05;
x312.l = 10.4;
x313.l = 0.05;
x314.l = 0.05;
x315.l = 10.5;
x316.l = 0.05;
x317.l = 0.05;
x318.l = 10.6;
x319.l = 0.05;
x320.l = 0.05;
x321.l = 10.7;
x322.l = 0.05;
x323.l = 0.05;
x324.l = 10.8;
x325.l = 0.05;
x326.l = 0.05;
x327.l = 10.9;
x328.l = 0.05;
x329.l = 0.05;
x330.l = 11;
x331.l = 0.05;
x332.l = 0.05;
x333.l = 11.1;
x334.l = 0.05;
x335.l = 0.05;
x336.l = 11.2;
x337.l = 0.05;
x338.l = 0.05;
x339.l = 11.3;
x340.l = 0.05;
x341.l = 0.05;
x342.l = 11.4;
x343.l = 0.05;
x344.l = 0.05;
x345.l = 11.5;
x346.l = 0.05;
x347.l = 0.05;
x348.l = 11.6;
x349.l = 0.05;
x350.l = 0.05;
x351.l = 11.7;
x352.l = 0.05;
x353.l = 0.05;
x354.l = 11.8;
x355.l = 0.05;
x356.l = 0.05;
x357.l = 11.9;
x358.l = 0.05;
x359.l = 0.05;
x360.l = 12;
x361.l = 0.05;
x362.l = 0.05;
x363.l = 12.1;
x364.l = 0.05;
x365.l = 0.05;
x366.l = 12.2;
x367.l = 0.05;
x368.l = 0.05;
x369.l = 12.3;
x370.l = 0.05;
x371.l = 0.05;
x372.l = 12.4;
x373.l = 0.05;
x374.l = 0.05;
x375.l = 12.5;
x376.l = 0.05;
x377.l = 0.05;
x378.l = 12.6;
x379.l = 0.05;
x380.l = 0.05;
x381.l = 12.7;
x382.l = 0.05;
x383.l = 0.05;
x384.l = 12.8;
x385.l = 0.05;
x386.l = 0.05;
x387.l = 12.9;
x388.l = 0.05;
x389.l = 0.05;
x390.l = 13;
x391.l = 0.05;
x392.l = 0.05;
x393.l = 13.1;
x394.l = 0.05;
x395.l = 0.05;
x396.l = 13.2;
x397.l = 0.05;
x398.l = 0.05;
x399.l = 13.3;
x400.l = 0.05;
x401.l = 0.05;
x402.l = 13.4;
x403.l = 0.05;
x404.l = 0.05;
x405.l = 13.5;
x406.l = 0.05;
x407.l = 0.05;
x408.l = 13.6;
x409.l = 0.05;
x410.l = 0.05;
x411.l = 13.7;
x412.l = 0.05;
x413.l = 0.05;
x414.l = 13.8;
x415.l = 0.05;
x416.l = 0.05;
x417.l = 13.9;
x418.l = 0.05;
x419.l = 0.05;
x420.l = 14;
x421.l = 0.05;
x422.l = 0.05;
x423.l = 14.1;
x424.l = 0.05;
x425.l = 0.05;
x426.l = 14.2;
x427.l = 0.05;
x428.l = 0.05;
x429.l = 14.3;
x430.l = 0.05;
x431.l = 0.05;
x432.l = 14.4;
x433.l = 0.05;
x434.l = 0.05;
x435.l = 14.5;
x436.l = 0.05;
x437.l = 0.05;
x438.l = 14.6;
x439.l = 0.05;
x440.l = 0.05;
x441.l = 14.7;
x442.l = 0.05;
x443.l = 0.05;
x444.l = 14.8;
x445.l = 0.05;
x446.l = 0.05;
x447.l = 14.9;
x448.l = 0.05;
x449.l = 0.05;
x450.l = 15;
x451.l = 0.05;
x452.l = 0.05;
x453.l = 15.1;
x454.l = 0.05;
x455.l = 0.05;
x456.l = 15.2;
x457.l = 0.05;
x458.l = 0.05;
x459.l = 15.3;
x460.l = 0.05;
x461.l = 0.05;
x462.l = 15.4;
x463.l = 0.05;
x464.l = 0.05;
x465.l = 15.5;
x466.l = 0.05;
x467.l = 0.05;
x468.l = 15.6;
x469.l = 0.05;
x470.l = 0.05;
x471.l = 15.7;
x472.l = 0.05;
x473.l = 0.05;
x474.l = 15.8;
x475.l = 0.05;
x476.l = 0.05;
x477.l = 15.9;
x478.l = 0.05;
x479.l = 0.05;
x480.l = 16;
x481.l = 0.05;
x482.l = 0.05;
x483.l = 16.1;
x484.l = 0.05;
x485.l = 0.05;
x486.l = 16.2;
x487.l = 0.05;
x488.l = 0.05;
x489.l = 16.3;
x490.l = 0.05;
x491.l = 0.05;
x492.l = 16.4;
x493.l = 0.05;
x494.l = 0.05;
x495.l = 16.5;
x496.l = 0.05;
x497.l = 0.05;
x498.l = 16.6;
x499.l = 0.05;
x500.l = 0.05;
x501.l = 16.7;
x502.l = 0.05;
x503.l = 0.05;
x504.l = 16.8;
x505.l = 0.05;
x506.l = 0.05;
x507.l = 16.9;
x508.l = 0.05;
x509.l = 0.05;
x510.l = 17;
x511.l = 0.05;
x512.l = 0.05;
x513.l = 17.1;
x514.l = 0.05;
x515.l = 0.05;
x516.l = 17.2;
x517.l = 0.05;
x518.l = 0.05;
x519.l = 17.3;
x520.l = 0.05;
x521.l = 0.05;
x522.l = 17.4;
x523.l = 0.05;
x524.l = 0.05;
x525.l = 17.5;
x526.l = 0.05;
x527.l = 0.05;
x528.l = 17.6;
x529.l = 0.05;
x530.l = 0.05;
x531.l = 17.7;
x532.l = 0.05;
x533.l = 0.05;
x534.l = 17.8;
x535.l = 0.05;
x536.l = 0.05;
x537.l = 17.9;
x538.l = 0.05;
x539.l = 0.05;
x540.l = 18;
x541.l = 0.05;
x542.l = 0.05;
x543.l = 18.1;
x544.l = 0.05;
x545.l = 0.05;
x546.l = 18.2;
x547.l = 0.05;
x548.l = 0.05;
x549.l = 18.3;
x550.l = 0.05;
x551.l = 0.05;
x552.l = 18.4;
x553.l = 0.05;
x554.l = 0.05;
x555.l = 18.5;
x556.l = 0.05;
x557.l = 0.05;
x558.l = 18.6;
x559.l = 0.05;
x560.l = 0.05;
x561.l = 18.7;
x562.l = 0.05;
x563.l = 0.05;
x564.l = 18.8;
x565.l = 0.05;
x566.l = 0.05;
x567.l = 18.9;
x568.l = 0.05;
x569.l = 0.05;
x570.l = 19;
x571.l = 0.05;
x572.l = 0.05;
x573.l = 19.1;
x574.l = 0.05;
x575.l = 0.05;
x576.l = 19.2;
x577.l = 0.05;
x578.l = 0.05;
x579.l = 19.3;
x580.l = 0.05;
x581.l = 0.05;
x582.l = 19.4;
x583.l = 0.05;
x584.l = 0.05;
x585.l = 19.5;
x586.l = 0.05;
x587.l = 0.05;
x588.l = 19.6;
x589.l = 0.05;
x590.l = 0.05;
x591.l = 19.7;
x592.l = 0.05;
x593.l = 0.05;
x594.l = 19.8;
x595.l = 0.05;
x596.l = 0.05;
x597.l = 19.9;
x598.l = 0.05;
x599.l = 0.05;
x600.l = 20;
* 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;
]]></plaintext></body>