springs_100_nlp.gms:
References:
- R. Vanderbei, Nonlinear Optimization Models (AMPL), See http://www.princeton.edu/~rvdb/ampl/nlmodels/.
- Yu-Ju Kuo, and H. D. Mittelmann, Interior Point Methods for Second Order Cone Programming and OR Applications. Computational Optimization and Applications (to appear).
- Original source: GAMS Model spring.gms from http://www.princeton.edu/~rvdb/ampl/nlmodels/springs/springs.mod
Point:
* LP written by GAMS Convert at 10/14/03 16:55:42
*
* Equation counts
* Total E G L N X C
* 406 305 0 0 0 0 101
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 609 609 0 0 0 0 0 0
* FX 4 4 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 1310 1310 0 0
*
* Solve m using LP minimizing objvar;
Variables objvar,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,x601,x602,x603,x604,x605,x606,x607,x608,x609;
Positive Variables x2,x103,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,x601,x602,x603,x604,x605,x606,x607,x608
,x609;
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,e57,e58,e59,e60,e61,e62,e63,e64,e65,e66,e67,e68,e69,e70
,e71,e72,e73,e74,e75,e76,e77,e78,e79,e80,e81,e82,e83,e84,e85,e86,e87
,e88,e89,e90,e91,e92,e93,e94,e95,e96,e97,e98,e99,e100,e101,e102,e103
,e104,e105,e106,e107,e108,e109,e110,e111,e112,e113,e114,e115,e116
,e117,e118,e119,e120,e121,e122,e123,e124,e125,e126,e127,e128,e129
,e130,e131,e132,e133,e134,e135,e136,e137,e138,e139,e140,e141,e142
,e143,e144,e145,e146,e147,e148,e149,e150,e151,e152,e153,e154,e155
,e156,e157,e158,e159,e160,e161,e162,e163,e164,e165,e166,e167,e168
,e169,e170,e171,e172,e173,e174,e175,e176,e177,e178,e179,e180,e181
,e182,e183,e184,e185,e186,e187,e188,e189,e190,e191,e192,e193,e194
,e195,e196,e197,e198,e199,e200,e201,e202,e203,e204,e205,e206,e207
,e208,e209,e210,e211,e212,e213,e214,e215,e216,e217,e218,e219,e220
,e221,e222,e223,e224,e225,e226,e227,e228,e229,e230,e231,e232,e233
,e234,e235,e236,e237,e238,e239,e240,e241,e242,e243,e244,e245,e246
,e247,e248,e249,e250,e251,e252,e253,e254,e255,e256,e257,e258,e259
,e260,e261,e262,e263,e264,e265,e266,e267,e268,e269,e270,e271,e272
,e273,e274,e275,e276,e277,e278,e279,e280,e281,e282,e283,e284,e285
,e286,e287,e288,e289,e290,e291,e292,e293,e294,e295,e296,e297,e298
,e299,e300,e301,e302,e303,e304,e305,e306,e307,e308,e309,e310,e311
,e312,e313,e314,e315,e316,e317,e318,e319,e320,e321,e322,e323,e324
,e325,e326,e327,e328,e329,e330,e331,e332,e333,e334,e335,e336,e337
,e338,e339,e340,e341,e342,e343,e344,e345,e346,e347,e348,e349,e350
,e351,e352,e353,e354,e355,e356,e357,e358,e359,e360,e361,e362,e363
,e364,e365,e366,e367,e368,e369,e370,e371,e372,e373,e374,e375,e376
,e377,e378,e379,e380,e381,e382,e383,e384,e385,e386,e387,e388,e389
,e390,e391,e392,e393,e394,e395,e396,e397,e398,e399,e400,e401,e402
,e403,e404,e405,e406;
e1.. objvar - 9.8*x104 - 9.8*x105 - 9.8*x106 - 9.8*x107 - 9.8*x108
- 9.8*x109 - 9.8*x110 - 9.8*x111 - 9.8*x112 - 9.8*x113 - 9.8*x114
- 9.8*x115 - 9.8*x116 - 9.8*x117 - 9.8*x118 - 9.8*x119 - 9.8*x120
- 9.8*x121 - 9.8*x122 - 9.8*x123 - 9.8*x124 - 9.8*x125 - 9.8*x126
- 9.8*x127 - 9.8*x128 - 9.8*x129 - 9.8*x130 - 9.8*x131 - 9.8*x132
- 9.8*x133 - 9.8*x134 - 9.8*x135 - 9.8*x136 - 9.8*x137 - 9.8*x138
- 9.8*x139 - 9.8*x140 - 9.8*x141 - 9.8*x142 - 9.8*x143 - 9.8*x144
- 9.8*x145 - 9.8*x146 - 9.8*x147 - 9.8*x148 - 9.8*x149 - 9.8*x150
- 9.8*x151 - 9.8*x152 - 9.8*x153 - 9.8*x154 - 9.8*x155 - 9.8*x156
- 9.8*x157 - 9.8*x158 - 9.8*x159 - 9.8*x160 - 9.8*x161 - 9.8*x162
- 9.8*x163 - 9.8*x164 - 9.8*x165 - 9.8*x166 - 9.8*x167 - 9.8*x168
- 9.8*x169 - 9.8*x170 - 9.8*x171 - 9.8*x172 - 9.8*x173 - 9.8*x174
- 9.8*x175 - 9.8*x176 - 9.8*x177 - 9.8*x178 - 9.8*x179 - 9.8*x180
- 9.8*x181 - 9.8*x182 - 9.8*x183 - 9.8*x184 - 9.8*x185 - 9.8*x186
- 9.8*x187 - 9.8*x188 - 9.8*x189 - 9.8*x190 - 9.8*x191 - 9.8*x192
- 9.8*x193 - 9.8*x194 - 9.8*x195 - 9.8*x196 - 9.8*x197 - 9.8*x198
- 9.8*x199 - 9.8*x200 - 9.8*x201 - 9.8*x202 - 100*x609 =E= 0;
e2.. - x2 + x306 =E= 0;
e3.. x2 - x3 + x307 =E= 0;
e4.. x3 - x4 + x308 =E= 0;
e5.. x4 - x5 + x309 =E= 0;
e6.. x5 - x6 + x310 =E= 0;
e7.. x6 - x7 + x311 =E= 0;
e8.. x7 - x8 + x312 =E= 0;
e9.. x8 - x9 + x313 =E= 0;
e10.. x9 - x10 + x314 =E= 0;
e11.. x10 - x11 + x315 =E= 0;
e12.. x11 - x12 + x316 =E= 0;
e13.. x12 - x13 + x317 =E= 0;
e14.. x13 - x14 + x318 =E= 0;
e15.. x14 - x15 + x319 =E= 0;
e16.. x15 - x16 + x320 =E= 0;
e17.. x16 - x17 + x321 =E= 0;
e18.. x17 - x18 + x322 =E= 0;
e19.. x18 - x19 + x323 =E= 0;
e20.. x19 - x20 + x324 =E= 0;
e21.. x20 - x21 + x325 =E= 0;
e22.. x21 - x22 + x326 =E= 0;
e23.. x22 - x23 + x327 =E= 0;
e24.. x23 - x24 + x328 =E= 0;
e25.. x24 - x25 + x329 =E= 0;
e26.. x25 - x26 + x330 =E= 0;
e27.. x26 - x27 + x331 =E= 0;
e28.. x27 - x28 + x332 =E= 0;
e29.. x28 - x29 + x333 =E= 0;
e30.. x29 - x30 + x334 =E= 0;
e31.. x30 - x31 + x335 =E= 0;
e32.. x31 - x32 + x336 =E= 0;
e33.. x32 - x33 + x337 =E= 0;
e34.. x33 - x34 + x338 =E= 0;
e35.. x34 - x35 + x339 =E= 0;
e36.. x35 - x36 + x340 =E= 0;
e37.. x36 - x37 + x341 =E= 0;
e38.. x37 - x38 + x342 =E= 0;
e39.. x38 - x39 + x343 =E= 0;
e40.. x39 - x40 + x344 =E= 0;
e41.. x40 - x41 + x345 =E= 0;
e42.. x41 - x42 + x346 =E= 0;
e43.. x42 - x43 + x347 =E= 0;
e44.. x43 - x44 + x348 =E= 0;
e45.. x44 - x45 + x349 =E= 0;
e46.. x45 - x46 + x350 =E= 0;
e47.. x46 - x47 + x351 =E= 0;
e48.. x47 - x48 + x352 =E= 0;
e49.. x48 - x49 + x353 =E= 0;
e50.. x49 - x50 + x354 =E= 0;
e51.. x50 - x51 + x355 =E= 0;
e52.. x51 - x52 + x356 =E= 0;
e53.. x52 - x53 + x357 =E= 0;
e54.. x53 - x54 + x358 =E= 0;
e55.. x54 - x55 + x359 =E= 0;
e56.. x55 - x56 + x360 =E= 0;
e57.. x56 - x57 + x361 =E= 0;
e58.. x57 - x58 + x362 =E= 0;
e59.. x58 - x59 + x363 =E= 0;
e60.. x59 - x60 + x364 =E= 0;
e61.. x60 - x61 + x365 =E= 0;
e62.. x61 - x62 + x366 =E= 0;
e63.. x62 - x63 + x367 =E= 0;
e64.. x63 - x64 + x368 =E= 0;
e65.. x64 - x65 + x369 =E= 0;
e66.. x65 - x66 + x370 =E= 0;
e67.. x66 - x67 + x371 =E= 0;
e68.. x67 - x68 + x372 =E= 0;
e69.. x68 - x69 + x373 =E= 0;
e70.. x69 - x70 + x374 =E= 0;
e71.. x70 - x71 + x375 =E= 0;
e72.. x71 - x72 + x376 =E= 0;
e73.. x72 - x73 + x377 =E= 0;
e74.. x73 - x74 + x378 =E= 0;
e75.. x74 - x75 + x379 =E= 0;
e76.. x75 - x76 + x380 =E= 0;
e77.. x76 - x77 + x381 =E= 0;
e78.. x77 - x78 + x382 =E= 0;
e79.. x78 - x79 + x383 =E= 0;
e80.. x79 - x80 + x384 =E= 0;
e81.. x80 - x81 + x385 =E= 0;
e82.. x81 - x82 + x386 =E= 0;
e83.. x82 - x83 + x387 =E= 0;
e84.. x83 - x84 + x388 =E= 0;
e85.. x84 - x85 + x389 =E= 0;
e86.. x85 - x86 + x390 =E= 0;
e87.. x86 - x87 + x391 =E= 0;
e88.. x87 - x88 + x392 =E= 0;
e89.. x88 - x89 + x393 =E= 0;
e90.. x89 - x90 + x394 =E= 0;
e91.. x90 - x91 + x395 =E= 0;
e92.. x91 - x92 + x396 =E= 0;
e93.. x92 - x93 + x397 =E= 0;
e94.. x93 - x94 + x398 =E= 0;
e95.. x94 - x95 + x399 =E= 0;
e96.. x95 - x96 + x400 =E= 0;
e97.. x96 - x97 + x401 =E= 0;
e98.. x97 - x98 + x402 =E= 0;
e99.. x98 - x99 + x403 =E= 0;
e100.. x99 - x100 + x404 =E= 0;
e101.. x100 - x101 + x405 =E= 0;
e102.. x101 - x102 + x406 =E= 0;
e103.. - x103 + x407 =E= 0;
e104.. x103 - x104 + x408 =E= 0;
e105.. x104 - x105 + x409 =E= 0;
e106.. x105 - x106 + x410 =E= 0;
e107.. x106 - x107 + x411 =E= 0;
e108.. x107 - x108 + x412 =E= 0;
e109.. x108 - x109 + x413 =E= 0;
e110.. x109 - x110 + x414 =E= 0;
e111.. x110 - x111 + x415 =E= 0;
e112.. x111 - x112 + x416 =E= 0;
e113.. x112 - x113 + x417 =E= 0;
e114.. x113 - x114 + x418 =E= 0;
e115.. x114 - x115 + x419 =E= 0;
e116.. x115 - x116 + x420 =E= 0;
e117.. x116 - x117 + x421 =E= 0;
e118.. x117 - x118 + x422 =E= 0;
e119.. x118 - x119 + x423 =E= 0;
e120.. x119 - x120 + x424 =E= 0;
e121.. x120 - x121 + x425 =E= 0;
e122.. x121 - x122 + x426 =E= 0;
e123.. x122 - x123 + x427 =E= 0;
e124.. x123 - x124 + x428 =E= 0;
e125.. x124 - x125 + x429 =E= 0;
e126.. x125 - x126 + x430 =E= 0;
e127.. x126 - x127 + x431 =E= 0;
e128.. x127 - x128 + x432 =E= 0;
e129.. x128 - x129 + x433 =E= 0;
e130.. x129 - x130 + x434 =E= 0;
e131.. x130 - x131 + x435 =E= 0;
e132.. x131 - x132 + x436 =E= 0;
e133.. x132 - x133 + x437 =E= 0;
e134.. x133 - x134 + x438 =E= 0;
e135.. x134 - x135 + x439 =E= 0;
e136.. x135 - x136 + x440 =E= 0;
e137.. x136 - x137 + x441 =E= 0;
e138.. x137 - x138 + x442 =E= 0;
e139.. x138 - x139 + x443 =E= 0;
e140.. x139 - x140 + x444 =E= 0;
e141.. x140 - x141 + x445 =E= 0;
e142.. x141 - x142 + x446 =E= 0;
e143.. x142 - x143 + x447 =E= 0;
e144.. x143 - x144 + x448 =E= 0;
e145.. x144 - x145 + x449 =E= 0;
e146.. x145 - x146 + x450 =E= 0;
e147.. x146 - x147 + x451 =E= 0;
e148.. x147 - x148 + x452 =E= 0;
e149.. x148 - x149 + x453 =E= 0;
e150.. x149 - x150 + x454 =E= 0;
e151.. x150 - x151 + x455 =E= 0;
e152.. x151 - x152 + x456 =E= 0;
e153.. x152 - x153 + x457 =E= 0;
e154.. x153 - x154 + x458 =E= 0;
e155.. x154 - x155 + x459 =E= 0;
e156.. x155 - x156 + x460 =E= 0;
e157.. x156 - x157 + x461 =E= 0;
e158.. x157 - x158 + x462 =E= 0;
e159.. x158 - x159 + x463 =E= 0;
e160.. x159 - x160 + x464 =E= 0;
e161.. x160 - x161 + x465 =E= 0;
e162.. x161 - x162 + x466 =E= 0;
e163.. x162 - x163 + x467 =E= 0;
e164.. x163 - x164 + x468 =E= 0;
e165.. x164 - x165 + x469 =E= 0;
e166.. x165 - x166 + x470 =E= 0;
e167.. x166 - x167 + x471 =E= 0;
e168.. x167 - x168 + x472 =E= 0;
e169.. x168 - x169 + x473 =E= 0;
e170.. x169 - x170 + x474 =E= 0;
e171.. x170 - x171 + x475 =E= 0;
e172.. x171 - x172 + x476 =E= 0;
e173.. x172 - x173 + x477 =E= 0;
e174.. x173 - x174 + x478 =E= 0;
e175.. x174 - x175 + x479 =E= 0;
e176.. x175 - x176 + x480 =E= 0;
e177.. x176 - x177 + x481 =E= 0;
e178.. x177 - x178 + x482 =E= 0;
e179.. x178 - x179 + x483 =E= 0;
e180.. x179 - x180 + x484 =E= 0;
e181.. x180 - x181 + x485 =E= 0;
e182.. x181 - x182 + x486 =E= 0;
e183.. x182 - x183 + x487 =E= 0;
e184.. x183 - x184 + x488 =E= 0;
e185.. x184 - x185 + x489 =E= 0;
e186.. x185 - x186 + x490 =E= 0;
e187.. x186 - x187 + x491 =E= 0;
e188.. x187 - x188 + x492 =E= 0;
e189.. x188 - x189 + x493 =E= 0;
e190.. x189 - x190 + x494 =E= 0;
e191.. x190 - x191 + x495 =E= 0;
e192.. x191 - x192 + x496 =E= 0;
e193.. x192 - x193 + x497 =E= 0;
e194.. x193 - x194 + x498 =E= 0;
e195.. x194 - x195 + x499 =E= 0;
e196.. x195 - x196 + x500 =E= 0;
e197.. x196 - x197 + x501 =E= 0;
e198.. x197 - x198 + x502 =E= 0;
e199.. x198 - x199 + x503 =E= 0;
e200.. x199 - x200 + x504 =E= 0;
e201.. x200 - x201 + x505 =E= 0;
e202.. x201 - x202 + x506 =E= 0;
e203.. x202 - x203 + x507 =E= 0;
e204.. x305 =E= 1;
e205.. x204 - x508 =E= 0.0447213595499958;
e206.. x205 - x509 =E= 0.0447213595499958;
e207.. x206 - x510 =E= 0.0447213595499958;
e208.. x207 - x511 =E= 0.0447213595499958;
e209.. x208 - x512 =E= 0.0447213595499958;
e210.. x209 - x513 =E= 0.0447213595499958;
e211.. x210 - x514 =E= 0.0447213595499958;
e212.. x211 - x515 =E= 0.0447213595499958;
e213.. x212 - x516 =E= 0.0447213595499958;
e214.. x213 - x517 =E= 0.0447213595499958;
e215.. x214 - x518 =E= 0.0447213595499958;
e216.. x215 - x519 =E= 0.0447213595499958;
e217.. x216 - x520 =E= 0.0447213595499958;
e218.. x217 - x521 =E= 0.0447213595499958;
e219.. x218 - x522 =E= 0.0447213595499958;
e220.. x219 - x523 =E= 0.0447213595499958;
e221.. x220 - x524 =E= 0.0447213595499958;
e222.. x221 - x525 =E= 0.0447213595499958;
e223.. x222 - x526 =E= 0.0447213595499958;
e224.. x223 - x527 =E= 0.0447213595499958;
e225.. x224 - x528 =E= 0.0447213595499958;
e226.. x225 - x529 =E= 0.0447213595499958;
e227.. x226 - x530 =E= 0.0447213595499958;
e228.. x227 - x531 =E= 0.0447213595499958;
e229.. x228 - x532 =E= 0.0447213595499958;
e230.. x229 - x533 =E= 0.0447213595499958;
e231.. x230 - x534 =E= 0.0447213595499958;
e232.. x231 - x535 =E= 0.0447213595499958;
e233.. x232 - x536 =E= 0.0447213595499958;
e234.. x233 - x537 =E= 0.0447213595499958;
e235.. x234 - x538 =E= 0.0447213595499958;
e236.. x235 - x539 =E= 0.0447213595499958;
e237.. x236 - x540 =E= 0.0447213595499958;
e238.. x237 - x541 =E= 0.0447213595499958;
e239.. x238 - x542 =E= 0.0447213595499958;
e240.. x239 - x543 =E= 0.0447213595499958;
e241.. x240 - x544 =E= 0.0447213595499958;
e242.. x241 - x545 =E= 0.0447213595499958;
e243.. x242 - x546 =E= 0.0447213595499958;
e244.. x243 - x547 =E= 0.0447213595499958;
e245.. x244 - x548 =E= 0.0447213595499958;
e246.. x245 - x549 =E= 0.0447213595499958;
e247.. x246 - x550 =E= 0.0447213595499958;
e248.. x247 - x551 =E= 0.0447213595499958;
e249.. x248 - x552 =E= 0.0447213595499958;
e250.. x249 - x553 =E= 0.0447213595499958;
e251.. x250 - x554 =E= 0.0447213595499958;
e252.. x251 - x555 =E= 0.0447213595499958;
e253.. x252 - x556 =E= 0.0447213595499958;
e254.. x253 - x557 =E= 0.0447213595499958;
e255.. x254 - x558 =E= 0.0447213595499958;
e256.. x255 - x559 =E= 0.0447213595499958;
e257.. x256 - x560 =E= 0.0447213595499958;
e258.. x257 - x561 =E= 0.0447213595499958;
e259.. x258 - x562 =E= 0.0447213595499958;
e260.. x259 - x563 =E= 0.0447213595499958;
e261.. x260 - x564 =E= 0.0447213595499958;
e262.. x261 - x565 =E= 0.0447213595499958;
e263.. x262 - x566 =E= 0.0447213595499958;
e264.. x263 - x567 =E= 0.0447213595499958;
e265.. x264 - x568 =E= 0.0447213595499958;
e266.. x265 - x569 =E= 0.0447213595499958;
e267.. x266 - x570 =E= 0.0447213595499958;
e268.. x267 - x571 =E= 0.0447213595499958;
e269.. x268 - x572 =E= 0.0447213595499958;
e270.. x269 - x573 =E= 0.0447213595499958;
e271.. x270 - x574 =E= 0.0447213595499958;
e272.. x271 - x575 =E= 0.0447213595499958;
e273.. x272 - x576 =E= 0.0447213595499958;
e274.. x273 - x577 =E= 0.0447213595499958;
e275.. x274 - x578 =E= 0.0447213595499958;
e276.. x275 - x579 =E= 0.0447213595499958;
e277.. x276 - x580 =E= 0.0447213595499958;
e278.. x277 - x581 =E= 0.0447213595499958;
e279.. x278 - x582 =E= 0.0447213595499958;
e280.. x279 - x583 =E= 0.0447213595499958;
e281.. x280 - x584 =E= 0.0447213595499958;
e282.. x281 - x585 =E= 0.0447213595499958;
e283.. x282 - x586 =E= 0.0447213595499958;
e284.. x283 - x587 =E= 0.0447213595499958;
e285.. x284 - x588 =E= 0.0447213595499958;
e286.. x285 - x589 =E= 0.0447213595499958;
e287.. x286 - x590 =E= 0.0447213595499958;
e288.. x287 - x591 =E= 0.0447213595499958;
e289.. x288 - x592 =E= 0.0447213595499958;
e290.. x289 - x593 =E= 0.0447213595499958;
e291.. x290 - x594 =E= 0.0447213595499958;
e292.. x291 - x595 =E= 0.0447213595499958;
e293.. x292 - x596 =E= 0.0447213595499958;
e294.. x293 - x597 =E= 0.0447213595499958;
e295.. x294 - x598 =E= 0.0447213595499958;
e296.. x295 - x599 =E= 0.0447213595499958;
e297.. x296 - x600 =E= 0.0447213595499958;
e298.. x297 - x601 =E= 0.0447213595499958;
e299.. x298 - x602 =E= 0.0447213595499958;
e300.. x299 - x603 =E= 0.0447213595499958;
e301.. x300 - x604 =E= 0.0447213595499958;
e302.. x301 - x605 =E= 0.0447213595499958;
e303.. x302 - x606 =E= 0.0447213595499958;
e304.. x303 - x607 =E= 0.0447213595499958;
e305.. x304 - x608 =E= 0.0447213595499958;
e306.. x205 =G= SQRT(SQR(x307)+SQR(x408));
e307.. x206 =G= SQRT(SQR(x308)+SQR(x409));
e308.. x207 =G= SQRT(SQR(x309)+SQR(x410));
e309.. x208 =G= SQRT(SQR(x310)+SQR(x411));
e310.. x209 =G= SQRT(SQR(x311)+SQR(x412));
e311.. x210 =G= SQRT(SQR(x312)+SQR(x413));
e312.. x211 =G= SQRT(SQR(x313)+SQR(x414));
e313.. x212 =G= SQRT(SQR(x314)+SQR(x415));
e314.. x213 =G= SQRT(SQR(x315)+SQR(x416));
e315.. x214 =G= SQRT(SQR(x316)+SQR(x417));
e316.. x215 =G= SQRT(SQR(x317)+SQR(x418));
e317.. x216 =G= SQRT(SQR(x318)+SQR(x419));
e318.. x217 =G= SQRT(SQR(x319)+SQR(x420));
e319.. x218 =G= SQRT(SQR(x320)+SQR(x421));
e320.. x219 =G= SQRT(SQR(x321)+SQR(x422));
e321.. x220 =G= SQRT(SQR(x322)+SQR(x423));
e322.. x221 =G= SQRT(SQR(x323)+SQR(x424));
e323.. x222 =G= SQRT(SQR(x324)+SQR(x425));
e324.. x223 =G= SQRT(SQR(x325)+SQR(x426));
e325.. x224 =G= SQRT(SQR(x326)+SQR(x427));
e326.. x225 =G= SQRT(SQR(x327)+SQR(x428));
e327.. x226 =G= SQRT(SQR(x328)+SQR(x429));
e328.. x227 =G= SQRT(SQR(x329)+SQR(x430));
e329.. x228 =G= SQRT(SQR(x330)+SQR(x431));
e330.. x229 =G= SQRT(SQR(x331)+SQR(x432));
e331.. x230 =G= SQRT(SQR(x332)+SQR(x433));
e332.. x231 =G= SQRT(SQR(x333)+SQR(x434));
e333.. x232 =G= SQRT(SQR(x334)+SQR(x435));
e334.. x233 =G= SQRT(SQR(x335)+SQR(x436));
e335.. x234 =G= SQRT(SQR(x336)+SQR(x437));
e336.. x235 =G= SQRT(SQR(x337)+SQR(x438));
e337.. x236 =G= SQRT(SQR(x338)+SQR(x439));
e338.. x237 =G= SQRT(SQR(x339)+SQR(x440));
e339.. x238 =G= SQRT(SQR(x340)+SQR(x441));
e340.. x239 =G= SQRT(SQR(x341)+SQR(x442));
e341.. x240 =G= SQRT(SQR(x342)+SQR(x443));
e342.. x241 =G= SQRT(SQR(x343)+SQR(x444));
e343.. x242 =G= SQRT(SQR(x344)+SQR(x445));
e344.. x243 =G= SQRT(SQR(x345)+SQR(x446));
e345.. x244 =G= SQRT(SQR(x346)+SQR(x447));
e346.. x245 =G= SQRT(SQR(x347)+SQR(x448));
e347.. x246 =G= SQRT(SQR(x348)+SQR(x449));
e348.. x247 =G= SQRT(SQR(x349)+SQR(x450));
e349.. x248 =G= SQRT(SQR(x350)+SQR(x451));
e350.. x249 =G= SQRT(SQR(x351)+SQR(x452));
e351.. x250 =G= SQRT(SQR(x352)+SQR(x453));
e352.. x251 =G= SQRT(SQR(x353)+SQR(x454));
e353.. x252 =G= SQRT(SQR(x354)+SQR(x455));
e354.. x253 =G= SQRT(SQR(x355)+SQR(x456));
e355.. x254 =G= SQRT(SQR(x356)+SQR(x457));
e356.. x255 =G= SQRT(SQR(x357)+SQR(x458));
e357.. x256 =G= SQRT(SQR(x358)+SQR(x459));
e358.. x257 =G= SQRT(SQR(x359)+SQR(x460));
e359.. x258 =G= SQRT(SQR(x360)+SQR(x461));
e360.. x259 =G= SQRT(SQR(x361)+SQR(x462));
e361.. x260 =G= SQRT(SQR(x362)+SQR(x463));
e362.. x261 =G= SQRT(SQR(x363)+SQR(x464));
e363.. x262 =G= SQRT(SQR(x364)+SQR(x465));
e364.. x263 =G= SQRT(SQR(x365)+SQR(x466));
e365.. x264 =G= SQRT(SQR(x366)+SQR(x467));
e366.. x265 =G= SQRT(SQR(x367)+SQR(x468));
e367.. x266 =G= SQRT(SQR(x368)+SQR(x469));
e368.. x267 =G= SQRT(SQR(x369)+SQR(x470));
e369.. x268 =G= SQRT(SQR(x370)+SQR(x471));
e370.. x269 =G= SQRT(SQR(x371)+SQR(x472));
e371.. x270 =G= SQRT(SQR(x372)+SQR(x473));
e372.. x271 =G= SQRT(SQR(x373)+SQR(x474));
e373.. x272 =G= SQRT(SQR(x374)+SQR(x475));
e374.. x273 =G= SQRT(SQR(x375)+SQR(x476));
e375.. x274 =G= SQRT(SQR(x376)+SQR(x477));
e376.. x275 =G= SQRT(SQR(x377)+SQR(x478));
e377.. x276 =G= SQRT(SQR(x378)+SQR(x479));
e378.. x277 =G= SQRT(SQR(x379)+SQR(x480));
e379.. x278 =G= SQRT(SQR(x380)+SQR(x481));
e380.. x279 =G= SQRT(SQR(x381)+SQR(x482));
e381.. x280 =G= SQRT(SQR(x382)+SQR(x483));
e382.. x281 =G= SQRT(SQR(x383)+SQR(x484));
e383.. x282 =G= SQRT(SQR(x384)+SQR(x485));
e384.. x283 =G= SQRT(SQR(x385)+SQR(x486));
e385.. x284 =G= SQRT(SQR(x386)+SQR(x487));
e386.. x285 =G= SQRT(SQR(x387)+SQR(x488));
e387.. x286 =G= SQRT(SQR(x388)+SQR(x489));
e388.. x287 =G= SQRT(SQR(x389)+SQR(x490));
e389.. x288 =G= SQRT(SQR(x390)+SQR(x491));
e390.. x289 =G= SQRT(SQR(x391)+SQR(x492));
e391.. x290 =G= SQRT(SQR(x392)+SQR(x493));
e392.. x291 =G= SQRT(SQR(x393)+SQR(x494));
e393.. x292 =G= SQRT(SQR(x394)+SQR(x495));
e394.. x293 =G= SQRT(SQR(x395)+SQR(x496));
e395.. x294 =G= SQRT(SQR(x396)+SQR(x497));
e396.. x295 =G= SQRT(SQR(x397)+SQR(x498));
e397.. x296 =G= SQRT(SQR(x398)+SQR(x499));
e398.. x297 =G= SQRT(SQR(x399)+SQR(x500));
e399.. x298 =G= SQRT(SQR(x400)+SQR(x501));
e400.. x299 =G= SQRT(SQR(x401)+SQR(x502));
e401.. x300 =G= SQRT(SQR(x402)+SQR(x503));
e402.. x301 =G= SQRT(SQR(x403)+SQR(x504));
e403.. x302 =G= SQRT(SQR(x404)+SQR(x505));
e404.. x303 =G= SQRT(SQR(x405)+SQR(x506));
e405.. x304 =G= SQRT(SQR(x406)+SQR(x507));
e406.. x305*x609*2 =G= SQR(x509)+SQR(x510)+SQR(x511)+SQR(x512)+SQR(x513)
+SQR(x514)+SQR(x515)+SQR(x516)+SQR(x517)+SQR(x518)+SQR(x519)+SQR(x520)
+SQR(x521)+SQR(x522)+SQR(x523)+SQR(x524)+SQR(x525)+SQR(x526)+SQR(x527)
+SQR(x528)+SQR(x529)+SQR(x530)+SQR(x531)+SQR(x532)+SQR(x533)+SQR(x534)
+SQR(x535)+SQR(x536)+SQR(x537)+SQR(x538)+SQR(x539)+SQR(x540)+SQR(x541)
+SQR(x542)+SQR(x543)+SQR(x544)+SQR(x545)+SQR(x546)+SQR(x547)+SQR(x548)
+SQR(x549)+SQR(x550)+SQR(x551)+SQR(x552)+SQR(x553)+SQR(x554)+SQR(x555)
+SQR(x556)+SQR(x557)+SQR(x558)+SQR(x559)+SQR(x560)+SQR(x561)+SQR(x562)
+SQR(x563)+SQR(x564)+SQR(x565)+SQR(x566)+SQR(x567)+SQR(x568)+SQR(x569)
+SQR(x570)+SQR(x571)+SQR(x572)+SQR(x573)+SQR(x574)+SQR(x575)+SQR(x576)
+SQR(x577)+SQR(x578)+SQR(x579)+SQR(x580)+SQR(x581)+SQR(x582)+SQR(x583)
+SQR(x584)+SQR(x585)+SQR(x586)+SQR(x587)+SQR(x588)+SQR(x589)+SQR(x590)
+SQR(x591)+SQR(x592)+SQR(x593)+SQR(x594)+SQR(x595)+SQR(x596)+SQR(x597)
+SQR(x598)+SQR(x599)+SQR(x600)+SQR(x601)+SQR(x602)+SQR(x603)+SQR(x604)
+SQR(x605)+SQR(x606)+SQR(x607)+SQR(x608);
* set non default bounds
x2.fx = 0;
x102.fx = 2;
x103.fx = 0;
x203.fx = -1;
* set non default levels
x3.l = 0.02;
x4.l = 0.04;
x5.l = 0.06;
x6.l = 0.08;
x7.l = 0.1;
x8.l = 0.12;
x9.l = 0.14;
x10.l = 0.16;
x11.l = 0.18;
x12.l = 0.2;
x13.l = 0.22;
x14.l = 0.24;
x15.l = 0.26;
x16.l = 0.28;
x17.l = 0.3;
x18.l = 0.32;
x19.l = 0.34;
x20.l = 0.36;
x21.l = 0.38;
x22.l = 0.4;
x23.l = 0.42;
x24.l = 0.44;
x25.l = 0.46;
x26.l = 0.48;
x27.l = 0.5;
x28.l = 0.52;
x29.l = 0.54;
x30.l = 0.56;
x31.l = 0.58;
x32.l = 0.6;
x33.l = 0.62;
x34.l = 0.64;
x35.l = 0.66;
x36.l = 0.68;
x37.l = 0.7;
x38.l = 0.72;
x39.l = 0.74;
x40.l = 0.76;
x41.l = 0.78;
x42.l = 0.8;
x43.l = 0.82;
x44.l = 0.84;
x45.l = 0.86;
x46.l = 0.88;
x47.l = 0.9;
x48.l = 0.92;
x49.l = 0.94;
x50.l = 0.96;
x51.l = 0.98;
x52.l = 1;
x53.l = 1.02;
x54.l = 1.04;
x55.l = 1.06;
x56.l = 1.08;
x57.l = 1.1;
x58.l = 1.12;
x59.l = 1.14;
x60.l = 1.16;
x61.l = 1.18;
x62.l = 1.2;
x63.l = 1.22;
x64.l = 1.24;
x65.l = 1.26;
x66.l = 1.28;
x67.l = 1.3;
x68.l = 1.32;
x69.l = 1.34;
x70.l = 1.36;
x71.l = 1.38;
x72.l = 1.4;
x73.l = 1.42;
x74.l = 1.44;
x75.l = 1.46;
x76.l = 1.48;
x77.l = 1.5;
x78.l = 1.52;
x79.l = 1.54;
x80.l = 1.56;
x81.l = 1.58;
x82.l = 1.6;
x83.l = 1.62;
x84.l = 1.64;
x85.l = 1.66;
x86.l = 1.68;
x87.l = 1.7;
x88.l = 1.72;
x89.l = 1.74;
x90.l = 1.76;
x91.l = 1.78;
x92.l = 1.8;
x93.l = 1.82;
x94.l = 1.84;
x95.l = 1.86;
x96.l = 1.88;
x97.l = 1.9;
x98.l = 1.92;
x99.l = 1.94;
x100.l = 1.96;
x101.l = 1.98;
x104.l = -0.01;
x105.l = -0.02;
x106.l = -0.03;
x107.l = -0.04;
x108.l = -0.05;
x109.l = -0.06;
x110.l = -0.07;
x111.l = -0.08;
x112.l = -0.09;
x113.l = -0.1;
x114.l = -0.11;
x115.l = -0.12;
x116.l = -0.13;
x117.l = -0.14;
x118.l = -0.15;
x119.l = -0.16;
x120.l = -0.17;
x121.l = -0.18;
x122.l = -0.19;
x123.l = -0.2;
x124.l = -0.21;
x125.l = -0.22;
x126.l = -0.23;
x127.l = -0.24;
x128.l = -0.25;
x129.l = -0.26;
x130.l = -0.27;
x131.l = -0.28;
x132.l = -0.29;
x133.l = -0.3;
x134.l = -0.31;
x135.l = -0.32;
x136.l = -0.33;
x137.l = -0.34;
x138.l = -0.35;
x139.l = -0.36;
x140.l = -0.37;
x141.l = -0.38;
x142.l = -0.39;
x143.l = -0.4;
x144.l = -0.41;
x145.l = -0.42;
x146.l = -0.43;
x147.l = -0.44;
x148.l = -0.45;
x149.l = -0.46;
x150.l = -0.47;
x151.l = -0.48;
x152.l = -0.49;
x153.l = -0.5;
x154.l = -0.51;
x155.l = -0.52;
x156.l = -0.53;
x157.l = -0.54;
x158.l = -0.55;
x159.l = -0.56;
x160.l = -0.57;
x161.l = -0.58;
x162.l = -0.59;
x163.l = -0.6;
x164.l = -0.61;
x165.l = -0.62;
x166.l = -0.63;
x167.l = -0.64;
x168.l = -0.65;
x169.l = -0.66;
x170.l = -0.67;
x171.l = -0.68;
x172.l = -0.69;
x173.l = -0.7;
x174.l = -0.71;
x175.l = -0.72;
x176.l = -0.73;
x177.l = -0.74;
x178.l = -0.75;
x179.l = -0.76;
x180.l = -0.77;
x181.l = -0.78;
x182.l = -0.79;
x183.l = -0.8;
x184.l = -0.81;
x185.l = -0.82;
x186.l = -0.83;
x187.l = -0.84;
x188.l = -0.85;
x189.l = -0.86;
x190.l = -0.87;
x191.l = -0.88;
x192.l = -0.89;
x193.l = -0.9;
x194.l = -0.91;
x195.l = -0.92;
x196.l = -0.93;
x197.l = -0.94;
x198.l = -0.95;
x199.l = -0.96;
x200.l = -0.97;
x201.l = -0.98;
x202.l = -0.99;
* 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;