probpenl.gms
* NLP written by GAMS Convert at 10/06/06 11:47:15
*
* Equation counts
* Total E G L N X C
* 1 1 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 501 501 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 501 1 500 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,objvar;
Positive 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;
Equations e1;
e1.. - (0.002*exp(-x1*x2)*(0.0001*x1 + 0.0001*x2) + 0.002*exp(-x2*x3)*(0.0001*
x2 + 0.0001*x3) + 0.002*exp(-x3*x4)*(0.0001*x3 + 0.0001*x4) + 0.002*exp(-
x4*x5)*(0.0001*x4 + 0.0001*x5) + 0.002*exp(-x5*x6)*(0.0001*x5 + 0.0001*x6)
+ 0.002*exp(-x6*x7)*(0.0001*x6 + 0.0001*x7) + 0.002*exp(-x7*x8)*(0.0001*
x7 + 0.0001*x8) + 0.002*exp(-x8*x9)*(0.0001*x8 + 0.0001*x9) + 0.002*exp(-
x9*x10)*(0.0001*x9 + 0.0001*x10) + 0.002*exp(-x10*x11)*(0.0001*x10 +
0.0001*x11) + 0.002*exp(-x11*x12)*(0.0001*x11 + 0.0001*x12) + 0.002*exp(-
x12*x13)*(0.0001*x12 + 0.0001*x13) + 0.002*exp(-x13*x14)*(0.0001*x13 +
0.0001*x14) + 0.002*exp(-x14*x15)*(0.0001*x14 + 0.0001*x15) + 0.002*exp(-
x15*x16)*(0.0001*x15 + 0.0001*x16) + 0.002*exp(-x16*x17)*(0.0001*x16 +
0.0001*x17) + 0.002*exp(-x17*x18)*(0.0001*x17 + 0.0001*x18) + 0.002*exp(-
x18*x19)*(0.0001*x18 + 0.0001*x19) + 0.002*exp(-x19*x20)*(0.0001*x19 +
0.0001*x20) + 0.002*exp(-x20*x21)*(0.0001*x20 + 0.0001*x21) + 0.002*exp(-
x21*x22)*(0.0001*x21 + 0.0001*x22) + 0.002*exp(-x22*x23)*(0.0001*x22 +
0.0001*x23) + 0.002*exp(-x23*x24)*(0.0001*x23 + 0.0001*x24) + 0.002*exp(-
x24*x25)*(0.0001*x24 + 0.0001*x25) + 0.002*exp(-x25*x26)*(0.0001*x25 +
0.0001*x26) + 0.002*exp(-x26*x27)*(0.0001*x26 + 0.0001*x27) + 0.002*exp(-
x27*x28)*(0.0001*x27 + 0.0001*x28) + 0.002*exp(-x28*x29)*(0.0001*x28 +
0.0001*x29) + 0.002*exp(-x29*x30)*(0.0001*x29 + 0.0001*x30) + 0.002*exp(-
x30*x31)*(0.0001*x30 + 0.0001*x31) + 0.002*exp(-x31*x32)*(0.0001*x31 +
0.0001*x32) + 0.002*exp(-x32*x33)*(0.0001*x32 + 0.0001*x33) + 0.002*exp(-
x33*x34)*(0.0001*x33 + 0.0001*x34) + 0.002*exp(-x34*x35)*(0.0001*x34 +
0.0001*x35) + 0.002*exp(-x35*x36)*(0.0001*x35 + 0.0001*x36) + 0.002*exp(-
x36*x37)*(0.0001*x36 + 0.0001*x37) + 0.002*exp(-x37*x38)*(0.0001*x37 +
0.0001*x38) + 0.002*exp(-x38*x39)*(0.0001*x38 + 0.0001*x39) + 0.002*exp(-
x39*x40)*(0.0001*x39 + 0.0001*x40) + 0.002*exp(-x40*x41)*(0.0001*x40 +
0.0001*x41) + 0.002*exp(-x41*x42)*(0.0001*x41 + 0.0001*x42) + 0.002*exp(-
x42*x43)*(0.0001*x42 + 0.0001*x43) + 0.002*exp(-x43*x44)*(0.0001*x43 +
0.0001*x44) + 0.002*exp(-x44*x45)*(0.0001*x44 + 0.0001*x45) + 0.002*exp(-
x45*x46)*(0.0001*x45 + 0.0001*x46) + 0.002*exp(-x46*x47)*(0.0001*x46 +
0.0001*x47) + 0.002*exp(-x47*x48)*(0.0001*x47 + 0.0001*x48) + 0.002*exp(-
x48*x49)*(0.0001*x48 + 0.0001*x49) + 0.002*exp(-x49*x50)*(0.0001*x49 +
0.0001*x50) + 0.002*exp(-x50*x51)*(0.0001*x50 + 0.0001*x51) + 0.002*exp(-
x51*x52)*(0.0001*x51 + 0.0001*x52) + 0.002*exp(-x52*x53)*(0.0001*x52 +
0.0001*x53) + 0.002*exp(-x53*x54)*(0.0001*x53 + 0.0001*x54) + 0.002*exp(-
x54*x55)*(0.0001*x54 + 0.0001*x55) + 0.002*exp(-x55*x56)*(0.0001*x55 +
0.0001*x56) + 0.002*exp(-x56*x57)*(0.0001*x56 + 0.0001*x57) + 0.002*exp(-
x57*x58)*(0.0001*x57 + 0.0001*x58) + 0.002*exp(-x58*x59)*(0.0001*x58 +
0.0001*x59) + 0.002*exp(-x59*x60)*(0.0001*x59 + 0.0001*x60) + 0.002*exp(-
x60*x61)*(0.0001*x60 + 0.0001*x61) + 0.002*exp(-x61*x62)*(0.0001*x61 +
0.0001*x62) + 0.002*exp(-x62*x63)*(0.0001*x62 + 0.0001*x63) + 0.002*exp(-
x63*x64)*(0.0001*x63 + 0.0001*x64) + 0.002*exp(-x64*x65)*(0.0001*x64 +
0.0001*x65) + 0.002*exp(-x65*x66)*(0.0001*x65 + 0.0001*x66) + 0.002*exp(-
x66*x67)*(0.0001*x66 + 0.0001*x67) + 0.002*exp(-x67*x68)*(0.0001*x67 +
0.0001*x68) + 0.002*exp(-x68*x69)*(0.0001*x68 + 0.0001*x69) + 0.002*exp(-
x69*x70)*(0.0001*x69 + 0.0001*x70) + 0.002*exp(-x70*x71)*(0.0001*x70 +
0.0001*x71) + 0.002*exp(-x71*x72)*(0.0001*x71 + 0.0001*x72) + 0.002*exp(-
x72*x73)*(0.0001*x72 + 0.0001*x73) + 0.002*exp(-x73*x74)*(0.0001*x73 +
0.0001*x74) + 0.002*exp(-x74*x75)*(0.0001*x74 + 0.0001*x75) + 0.002*exp(-
x75*x76)*(0.0001*x75 + 0.0001*x76) + 0.002*exp(-x76*x77)*(0.0001*x76 +
0.0001*x77) + 0.002*exp(-x77*x78)*(0.0001*x77 + 0.0001*x78) + 0.002*exp(-
x78*x79)*(0.0001*x78 + 0.0001*x79) + 0.002*exp(-x79*x80)*(0.0001*x79 +
0.0001*x80) + 0.002*exp(-x80*x81)*(0.0001*x80 + 0.0001*x81) + 0.002*exp(-
x81*x82)*(0.0001*x81 + 0.0001*x82) + 0.002*exp(-x82*x83)*(0.0001*x82 +
0.0001*x83) + 0.002*exp(-x83*x84)*(0.0001*x83 + 0.0001*x84) + 0.002*exp(-
x84*x85)*(0.0001*x84 + 0.0001*x85) + 0.002*exp(-x85*x86)*(0.0001*x85 +
0.0001*x86) + 0.002*exp(-x86*x87)*(0.0001*x86 + 0.0001*x87) + 0.002*exp(-
x87*x88)*(0.0001*x87 + 0.0001*x88) + 0.002*exp(-x88*x89)*(0.0001*x88 +
0.0001*x89) + 0.002*exp(-x89*x90)*(0.0001*x89 + 0.0001*x90) + 0.002*exp(-
x90*x91)*(0.0001*x90 + 0.0001*x91) + 0.002*exp(-x91*x92)*(0.0001*x91 +
0.0001*x92) + 0.002*exp(-x92*x93)*(0.0001*x92 + 0.0001*x93) + 0.002*exp(-
x93*x94)*(0.0001*x93 + 0.0001*x94) + 0.002*exp(-x94*x95)*(0.0001*x94 +
0.0001*x95) + 0.002*exp(-x95*x96)*(0.0001*x95 + 0.0001*x96) + 0.002*exp(-
x96*x97)*(0.0001*x96 + 0.0001*x97) + 0.002*exp(-x97*x98)*(0.0001*x97 +
0.0001*x98) + 0.002*exp(-x98*x99)*(0.0001*x98 + 0.0001*x99) + 0.002*exp(-
x99*x100)*(0.0001*x99 + 0.0001*x100) + 0.002*exp(-x100*x101)*(0.0001*x100
+ 0.0001*x101) + 0.002*exp(-x101*x102)*(0.0001*x101 + 0.0001*x102) +
0.002*exp(-x102*x103)*(0.0001*x102 + 0.0001*x103) + 0.002*exp(-x103*x104)*
(0.0001*x103 + 0.0001*x104) + 0.002*exp(-x104*x105)*(0.0001*x104 + 0.0001*
x105) + 0.002*exp(-x105*x106)*(0.0001*x105 + 0.0001*x106) + 0.002*exp(-
x106*x107)*(0.0001*x106 + 0.0001*x107) + 0.002*exp(-x107*x108)*(0.0001*
x107 + 0.0001*x108) + 0.002*exp(-x108*x109)*(0.0001*x108 + 0.0001*x109) +
0.002*exp(-x109*x110)*(0.0001*x109 + 0.0001*x110) + 0.002*exp(-x110*x111)*
(0.0001*x110 + 0.0001*x111) + 0.002*exp(-x111*x112)*(0.0001*x111 + 0.0001*
x112) + 0.002*exp(-x112*x113)*(0.0001*x112 + 0.0001*x113) + 0.002*exp(-
x113*x114)*(0.0001*x113 + 0.0001*x114) + 0.002*exp(-x114*x115)*(0.0001*
x114 + 0.0001*x115) + 0.002*exp(-x115*x116)*(0.0001*x115 + 0.0001*x116) +
0.002*exp(-x116*x117)*(0.0001*x116 + 0.0001*x117) + 0.002*exp(-x117*x118)*
(0.0001*x117 + 0.0001*x118) + 0.002*exp(-x118*x119)*(0.0001*x118 + 0.0001*
x119) + 0.002*exp(-x119*x120)*(0.0001*x119 + 0.0001*x120) + 0.002*exp(-
x120*x121)*(0.0001*x120 + 0.0001*x121) + 0.002*exp(-x121*x122)*(0.0001*
x121 + 0.0001*x122) + 0.002*exp(-x122*x123)*(0.0001*x122 + 0.0001*x123) +
0.002*exp(-x123*x124)*(0.0001*x123 + 0.0001*x124) + 0.002*exp(-x124*x125)*
(0.0001*x124 + 0.0001*x125) + 0.002*exp(-x125*x126)*(0.0001*x125 + 0.0001*
x126) + 0.002*exp(-x126*x127)*(0.0001*x126 + 0.0001*x127) + 0.002*exp(-
x127*x128)*(0.0001*x127 + 0.0001*x128) + 0.002*exp(-x128*x129)*(0.0001*
x128 + 0.0001*x129) + 0.002*exp(-x129*x130)*(0.0001*x129 + 0.0001*x130) +
0.002*exp(-x130*x131)*(0.0001*x130 + 0.0001*x131) + 0.002*exp(-x131*x132)*
(0.0001*x131 + 0.0001*x132) + 0.002*exp(-x132*x133)*(0.0001*x132 + 0.0001*
x133) + 0.002*exp(-x133*x134)*(0.0001*x133 + 0.0001*x134) + 0.002*exp(-
x134*x135)*(0.0001*x134 + 0.0001*x135) + 0.002*exp(-x135*x136)*(0.0001*
x135 + 0.0001*x136) + 0.002*exp(-x136*x137)*(0.0001*x136 + 0.0001*x137) +
0.002*exp(-x137*x138)*(0.0001*x137 + 0.0001*x138) + 0.002*exp(-x138*x139)*
(0.0001*x138 + 0.0001*x139) + 0.002*exp(-x139*x140)*(0.0001*x139 + 0.0001*
x140) + 0.002*exp(-x140*x141)*(0.0001*x140 + 0.0001*x141) + 0.002*exp(-
x141*x142)*(0.0001*x141 + 0.0001*x142) + 0.002*exp(-x142*x143)*(0.0001*
x142 + 0.0001*x143) + 0.002*exp(-x143*x144)*(0.0001*x143 + 0.0001*x144) +
0.002*exp(-x144*x145)*(0.0001*x144 + 0.0001*x145) + 0.002*exp(-x145*x146)*
(0.0001*x145 + 0.0001*x146) + 0.002*exp(-x146*x147)*(0.0001*x146 + 0.0001*
x147) + 0.002*exp(-x147*x148)*(0.0001*x147 + 0.0001*x148) + 0.002*exp(-
x148*x149)*(0.0001*x148 + 0.0001*x149) + 0.002*exp(-x149*x150)*(0.0001*
x149 + 0.0001*x150) + 0.002*exp(-x150*x151)*(0.0001*x150 + 0.0001*x151) +
0.002*exp(-x151*x152)*(0.0001*x151 + 0.0001*x152) + 0.002*exp(-x152*x153)*
(0.0001*x152 + 0.0001*x153) + 0.002*exp(-x153*x154)*(0.0001*x153 + 0.0001*
x154) + 0.002*exp(-x154*x155)*(0.0001*x154 + 0.0001*x155) + 0.002*exp(-
x155*x156)*(0.0001*x155 + 0.0001*x156) + 0.002*exp(-x156*x157)*(0.0001*
x156 + 0.0001*x157) + 0.002*exp(-x157*x158)*(0.0001*x157 + 0.0001*x158) +
0.002*exp(-x158*x159)*(0.0001*x158 + 0.0001*x159) + 0.002*exp(-x159*x160)*
(0.0001*x159 + 0.0001*x160) + 0.002*exp(-x160*x161)*(0.0001*x160 + 0.0001*
x161) + 0.002*exp(-x161*x162)*(0.0001*x161 + 0.0001*x162) + 0.002*exp(-
x162*x163)*(0.0001*x162 + 0.0001*x163) + 0.002*exp(-x163*x164)*(0.0001*
x163 + 0.0001*x164) + 0.002*exp(-x164*x165)*(0.0001*x164 + 0.0001*x165) +
0.002*exp(-x165*x166)*(0.0001*x165 + 0.0001*x166) + 0.002*exp(-x166*x167)*
(0.0001*x166 + 0.0001*x167) + 0.002*exp(-x167*x168)*(0.0001*x167 + 0.0001*
x168) + 0.002*exp(-x168*x169)*(0.0001*x168 + 0.0001*x169) + 0.002*exp(-
x169*x170)*(0.0001*x169 + 0.0001*x170) + 0.002*exp(-x170*x171)*(0.0001*
x170 + 0.0001*x171) + 0.002*exp(-x171*x172)*(0.0001*x171 + 0.0001*x172) +
0.002*exp(-x172*x173)*(0.0001*x172 + 0.0001*x173) + 0.002*exp(-x173*x174)*
(0.0001*x173 + 0.0001*x174) + 0.002*exp(-x174*x175)*(0.0001*x174 + 0.0001*
x175) + 0.002*exp(-x175*x176)*(0.0001*x175 + 0.0001*x176) + 0.002*exp(-
x176*x177)*(0.0001*x176 + 0.0001*x177) + 0.002*exp(-x177*x178)*(0.0001*
x177 + 0.0001*x178) + 0.002*exp(-x178*x179)*(0.0001*x178 + 0.0001*x179) +
0.002*exp(-x179*x180)*(0.0001*x179 + 0.0001*x180) + 0.002*exp(-x180*x181)*
(0.0001*x180 + 0.0001*x181) + 0.002*exp(-x181*x182)*(0.0001*x181 + 0.0001*
x182) + 0.002*exp(-x182*x183)*(0.0001*x182 + 0.0001*x183) + 0.002*exp(-
x183*x184)*(0.0001*x183 + 0.0001*x184) + 0.002*exp(-x184*x185)*(0.0001*
x184 + 0.0001*x185) + 0.002*exp(-x185*x186)*(0.0001*x185 + 0.0001*x186) +
0.002*exp(-x186*x187)*(0.0001*x186 + 0.0001*x187) + 0.002*exp(-x187*x188)*
(0.0001*x187 + 0.0001*x188) + 0.002*exp(-x188*x189)*(0.0001*x188 + 0.0001*
x189) + 0.002*exp(-x189*x190)*(0.0001*x189 + 0.0001*x190) + 0.002*exp(-
x190*x191)*(0.0001*x190 + 0.0001*x191) + 0.002*exp(-x191*x192)*(0.0001*
x191 + 0.0001*x192) + 0.002*exp(-x192*x193)*(0.0001*x192 + 0.0001*x193) +
0.002*exp(-x193*x194)*(0.0001*x193 + 0.0001*x194) + 0.002*exp(-x194*x195)*
(0.0001*x194 + 0.0001*x195) + 0.002*exp(-x195*x196)*(0.0001*x195 + 0.0001*
x196) + 0.002*exp(-x196*x197)*(0.0001*x196 + 0.0001*x197) + 0.002*exp(-
x197*x198)*(0.0001*x197 + 0.0001*x198) + 0.002*exp(-x198*x199)*(0.0001*
x198 + 0.0001*x199) + 0.002*exp(-x199*x200)*(0.0001*x199 + 0.0001*x200) +
0.002*exp(-x200*x201)*(0.0001*x200 + 0.0001*x201) + 0.002*exp(-x201*x202)*
(0.0001*x201 + 0.0001*x202) + 0.002*exp(-x202*x203)*(0.0001*x202 + 0.0001*
x203) + 0.002*exp(-x203*x204)*(0.0001*x203 + 0.0001*x204) + 0.002*exp(-
x204*x205)*(0.0001*x204 + 0.0001*x205) + 0.002*exp(-x205*x206)*(0.0001*
x205 + 0.0001*x206) + 0.002*exp(-x206*x207)*(0.0001*x206 + 0.0001*x207) +
0.002*exp(-x207*x208)*(0.0001*x207 + 0.0001*x208) + 0.002*exp(-x208*x209)*
(0.0001*x208 + 0.0001*x209) + 0.002*exp(-x209*x210)*(0.0001*x209 + 0.0001*
x210) + 0.002*exp(-x210*x211)*(0.0001*x210 + 0.0001*x211) + 0.002*exp(-
x211*x212)*(0.0001*x211 + 0.0001*x212) + 0.002*exp(-x212*x213)*(0.0001*
x212 + 0.0001*x213) + 0.002*exp(-x213*x214)*(0.0001*x213 + 0.0001*x214) +
0.002*exp(-x214*x215)*(0.0001*x214 + 0.0001*x215) + 0.002*exp(-x215*x216)*
(0.0001*x215 + 0.0001*x216) + 0.002*exp(-x216*x217)*(0.0001*x216 + 0.0001*
x217) + 0.002*exp(-x217*x218)*(0.0001*x217 + 0.0001*x218) + 0.002*exp(-
x218*x219)*(0.0001*x218 + 0.0001*x219) + 0.002*exp(-x219*x220)*(0.0001*
x219 + 0.0001*x220) + 0.002*exp(-x220*x221)*(0.0001*x220 + 0.0001*x221) +
0.002*exp(-x221*x222)*(0.0001*x221 + 0.0001*x222) + 0.002*exp(-x222*x223)*
(0.0001*x222 + 0.0001*x223) + 0.002*exp(-x223*x224)*(0.0001*x223 + 0.0001*
x224) + 0.002*exp(-x224*x225)*(0.0001*x224 + 0.0001*x225) + 0.002*exp(-
x225*x226)*(0.0001*x225 + 0.0001*x226) + 0.002*exp(-x226*x227)*(0.0001*
x226 + 0.0001*x227) + 0.002*exp(-x227*x228)*(0.0001*x227 + 0.0001*x228) +
0.002*exp(-x228*x229)*(0.0001*x228 + 0.0001*x229) + 0.002*exp(-x229*x230)*
(0.0001*x229 + 0.0001*x230) + 0.002*exp(-x230*x231)*(0.0001*x230 + 0.0001*
x231) + 0.002*exp(-x231*x232)*(0.0001*x231 + 0.0001*x232) + 0.002*exp(-
x232*x233)*(0.0001*x232 + 0.0001*x233) + 0.002*exp(-x233*x234)*(0.0001*
x233 + 0.0001*x234) + 0.002*exp(-x234*x235)*(0.0001*x234 + 0.0001*x235) +
0.002*exp(-x235*x236)*(0.0001*x235 + 0.0001*x236) + 0.002*exp(-x236*x237)*
(0.0001*x236 + 0.0001*x237) + 0.002*exp(-x237*x238)*(0.0001*x237 + 0.0001*
x238) + 0.002*exp(-x238*x239)*(0.0001*x238 + 0.0001*x239) + 0.002*exp(-
x239*x240)*(0.0001*x239 + 0.0001*x240) + 0.002*exp(-x240*x241)*(0.0001*
x240 + 0.0001*x241) + 0.002*exp(-x241*x242)*(0.0001*x241 + 0.0001*x242) +
0.002*exp(-x242*x243)*(0.0001*x242 + 0.0001*x243) + 0.002*exp(-x243*x244)*
(0.0001*x243 + 0.0001*x244) + 0.002*exp(-x244*x245)*(0.0001*x244 + 0.0001*
x245) + 0.002*exp(-x245*x246)*(0.0001*x245 + 0.0001*x246) + 0.002*exp(-
x246*x247)*(0.0001*x246 + 0.0001*x247) + 0.002*exp(-x247*x248)*(0.0001*
x247 + 0.0001*x248) + 0.002*exp(-x248*x249)*(0.0001*x248 + 0.0001*x249) +
0.002*exp(-x249*x250)*(0.0001*x249 + 0.0001*x250) + 0.002*exp(-x250*x251)*
(0.0001*x250 + 0.0001*x251) + 0.002*exp(-x251*x252)*(0.0001*x251 + 0.0001*
x252) + 0.002*exp(-x252*x253)*(0.0001*x252 + 0.0001*x253) + 0.002*exp(-
x253*x254)*(0.0001*x253 + 0.0001*x254) + 0.002*exp(-x254*x255)*(0.0001*
x254 + 0.0001*x255) + 0.002*exp(-x255*x256)*(0.0001*x255 + 0.0001*x256) +
0.002*exp(-x256*x257)*(0.0001*x256 + 0.0001*x257) + 0.002*exp(-x257*x258)*
(0.0001*x257 + 0.0001*x258) + 0.002*exp(-x258*x259)*(0.0001*x258 + 0.0001*
x259) + 0.002*exp(-x259*x260)*(0.0001*x259 + 0.0001*x260) + 0.002*exp(-
x260*x261)*(0.0001*x260 + 0.0001*x261) + 0.002*exp(-x261*x262)*(0.0001*
x261 + 0.0001*x262) + 0.002*exp(-x262*x263)*(0.0001*x262 + 0.0001*x263) +
0.002*exp(-x263*x264)*(0.0001*x263 + 0.0001*x264) + 0.002*exp(-x264*x265)*
(0.0001*x264 + 0.0001*x265) + 0.002*exp(-x265*x266)*(0.0001*x265 + 0.0001*
x266) + 0.002*exp(-x266*x267)*(0.0001*x266 + 0.0001*x267) + 0.002*exp(-
x267*x268)*(0.0001*x267 + 0.0001*x268) + 0.002*exp(-x268*x269)*(0.0001*
x268 + 0.0001*x269) + 0.002*exp(-x269*x270)*(0.0001*x269 + 0.0001*x270) +
0.002*exp(-x270*x271)*(0.0001*x270 + 0.0001*x271) + 0.002*exp(-x271*x272)*
(0.0001*x271 + 0.0001*x272) + 0.002*exp(-x272*x273)*(0.0001*x272 + 0.0001*
x273) + 0.002*exp(-x273*x274)*(0.0001*x273 + 0.0001*x274) + 0.002*exp(-
x274*x275)*(0.0001*x274 + 0.0001*x275) + 0.002*exp(-x275*x276)*(0.0001*
x275 + 0.0001*x276) + 0.002*exp(-x276*x277)*(0.0001*x276 + 0.0001*x277) +
0.002*exp(-x277*x278)*(0.0001*x277 + 0.0001*x278) + 0.002*exp(-x278*x279)*
(0.0001*x278 + 0.0001*x279) + 0.002*exp(-x279*x280)*(0.0001*x279 + 0.0001*
x280) + 0.002*exp(-x280*x281)*(0.0001*x280 + 0.0001*x281) + 0.002*exp(-
x281*x282)*(0.0001*x281 + 0.0001*x282) + 0.002*exp(-x282*x283)*(0.0001*
x282 + 0.0001*x283) + 0.002*exp(-x283*x284)*(0.0001*x283 + 0.0001*x284) +
0.002*exp(-x284*x285)*(0.0001*x284 + 0.0001*x285) + 0.002*exp(-x285*x286)*
(0.0001*x285 + 0.0001*x286) + 0.002*exp(-x286*x287)*(0.0001*x286 + 0.0001*
x287) + 0.002*exp(-x287*x288)*(0.0001*x287 + 0.0001*x288) + 0.002*exp(-
x288*x289)*(0.0001*x288 + 0.0001*x289) + 0.002*exp(-x289*x290)*(0.0001*
x289 + 0.0001*x290) + 0.002*exp(-x290*x291)*(0.0001*x290 + 0.0001*x291) +
0.002*exp(-x291*x292)*(0.0001*x291 + 0.0001*x292) + 0.002*exp(-x292*x293)*
(0.0001*x292 + 0.0001*x293) + 0.002*exp(-x293*x294)*(0.0001*x293 + 0.0001*
x294) + 0.002*exp(-x294*x295)*(0.0001*x294 + 0.0001*x295) + 0.002*exp(-
x295*x296)*(0.0001*x295 + 0.0001*x296) + 0.002*exp(-x296*x297)*(0.0001*
x296 + 0.0001*x297) + 0.002*exp(-x297*x298)*(0.0001*x297 + 0.0001*x298) +
0.002*exp(-x298*x299)*(0.0001*x298 + 0.0001*x299) + 0.002*exp(-x299*x300)*
(0.0001*x299 + 0.0001*x300) + 0.002*exp(-x300*x301)*(0.0001*x300 + 0.0001*
x301) + 0.002*exp(-x301*x302)*(0.0001*x301 + 0.0001*x302) + 0.002*exp(-
x302*x303)*(0.0001*x302 + 0.0001*x303) + 0.002*exp(-x303*x304)*(0.0001*
x303 + 0.0001*x304) + 0.002*exp(-x304*x305)*(0.0001*x304 + 0.0001*x305) +
0.002*exp(-x305*x306)*(0.0001*x305 + 0.0001*x306) + 0.002*exp(-x306*x307)*
(0.0001*x306 + 0.0001*x307) + 0.002*exp(-x307*x308)*(0.0001*x307 + 0.0001*
x308) + 0.002*exp(-x308*x309)*(0.0001*x308 + 0.0001*x309) + 0.002*exp(-
x309*x310)*(0.0001*x309 + 0.0001*x310) + 0.002*exp(-x310*x311)*(0.0001*
x310 + 0.0001*x311) + 0.002*exp(-x311*x312)*(0.0001*x311 + 0.0001*x312) +
0.002*exp(-x312*x313)*(0.0001*x312 + 0.0001*x313) + 0.002*exp(-x313*x314)*
(0.0001*x313 + 0.0001*x314) + 0.002*exp(-x314*x315)*(0.0001*x314 + 0.0001*
x315) + 0.002*exp(-x315*x316)*(0.0001*x315 + 0.0001*x316) + 0.002*exp(-
x316*x317)*(0.0001*x316 + 0.0001*x317) + 0.002*exp(-x317*x318)*(0.0001*
x317 + 0.0001*x318) + 0.002*exp(-x318*x319)*(0.0001*x318 + 0.0001*x319) +
0.002*exp(-x319*x320)*(0.0001*x319 + 0.0001*x320) + 0.002*exp(-x320*x321)*
(0.0001*x320 + 0.0001*x321) + 0.002*exp(-x321*x322)*(0.0001*x321 + 0.0001*
x322) + 0.002*exp(-x322*x323)*(0.0001*x322 + 0.0001*x323) + 0.002*exp(-
x323*x324)*(0.0001*x323 + 0.0001*x324) + 0.002*exp(-x324*x325)*(0.0001*
x324 + 0.0001*x325) + 0.002*exp(-x325*x326)*(0.0001*x325 + 0.0001*x326) +
0.002*exp(-x326*x327)*(0.0001*x326 + 0.0001*x327) + 0.002*exp(-x327*x328)*
(0.0001*x327 + 0.0001*x328) + 0.002*exp(-x328*x329)*(0.0001*x328 + 0.0001*
x329) + 0.002*exp(-x329*x330)*(0.0001*x329 + 0.0001*x330) + 0.002*exp(-
x330*x331)*(0.0001*x330 + 0.0001*x331) + 0.002*exp(-x331*x332)*(0.0001*
x331 + 0.0001*x332) + 0.002*exp(-x332*x333)*(0.0001*x332 + 0.0001*x333) +
0.002*exp(-x333*x334)*(0.0001*x333 + 0.0001*x334) + 0.002*exp(-x334*x335)*
(0.0001*x334 + 0.0001*x335) + 0.002*exp(-x335*x336)*(0.0001*x335 + 0.0001*
x336) + 0.002*exp(-x336*x337)*(0.0001*x336 + 0.0001*x337) + 0.002*exp(-
x337*x338)*(0.0001*x337 + 0.0001*x338) + 0.002*exp(-x338*x339)*(0.0001*
x338 + 0.0001*x339) + 0.002*exp(-x339*x340)*(0.0001*x339 + 0.0001*x340) +
0.002*exp(-x340*x341)*(0.0001*x340 + 0.0001*x341) + 0.002*exp(-x341*x342)*
(0.0001*x341 + 0.0001*x342) + 0.002*exp(-x342*x343)*(0.0001*x342 + 0.0001*
x343) + 0.002*exp(-x343*x344)*(0.0001*x343 + 0.0001*x344) + 0.002*exp(-
x344*x345)*(0.0001*x344 + 0.0001*x345) + 0.002*exp(-x345*x346)*(0.0001*
x345 + 0.0001*x346) + 0.002*exp(-x346*x347)*(0.0001*x346 + 0.0001*x347) +
0.002*exp(-x347*x348)*(0.0001*x347 + 0.0001*x348) + 0.002*exp(-x348*x349)*
(0.0001*x348 + 0.0001*x349) + 0.002*exp(-x349*x350)*(0.0001*x349 + 0.0001*
x350) + 0.002*exp(-x350*x351)*(0.0001*x350 + 0.0001*x351) + 0.002*exp(-
x351*x352)*(0.0001*x351 + 0.0001*x352) + 0.002*exp(-x352*x353)*(0.0001*
x352 + 0.0001*x353) + 0.002*exp(-x353*x354)*(0.0001*x353 + 0.0001*x354) +
0.002*exp(-x354*x355)*(0.0001*x354 + 0.0001*x355) + 0.002*exp(-x355*x356)*
(0.0001*x355 + 0.0001*x356) + 0.002*exp(-x356*x357)*(0.0001*x356 + 0.0001*
x357) + 0.002*exp(-x357*x358)*(0.0001*x357 + 0.0001*x358) + 0.002*exp(-
x358*x359)*(0.0001*x358 + 0.0001*x359) + 0.002*exp(-x359*x360)*(0.0001*
x359 + 0.0001*x360) + 0.002*exp(-x360*x361)*(0.0001*x360 + 0.0001*x361) +
0.002*exp(-x361*x362)*(0.0001*x361 + 0.0001*x362) + 0.002*exp(-x362*x363)*
(0.0001*x362 + 0.0001*x363) + 0.002*exp(-x363*x364)*(0.0001*x363 + 0.0001*
x364) + 0.002*exp(-x364*x365)*(0.0001*x364 + 0.0001*x365) + 0.002*exp(-
x365*x366)*(0.0001*x365 + 0.0001*x366) + 0.002*exp(-x366*x367)*(0.0001*
x366 + 0.0001*x367) + 0.002*exp(-x367*x368)*(0.0001*x367 + 0.0001*x368) +
0.002*exp(-x368*x369)*(0.0001*x368 + 0.0001*x369) + 0.002*exp(-x369*x370)*
(0.0001*x369 + 0.0001*x370) + 0.002*exp(-x370*x371)*(0.0001*x370 + 0.0001*
x371) + 0.002*exp(-x371*x372)*(0.0001*x371 + 0.0001*x372) + 0.002*exp(-
x372*x373)*(0.0001*x372 + 0.0001*x373) + 0.002*exp(-x373*x374)*(0.0001*
x373 + 0.0001*x374) + 0.002*exp(-x374*x375)*(0.0001*x374 + 0.0001*x375) +
0.002*exp(-x375*x376)*(0.0001*x375 + 0.0001*x376) + 0.002*exp(-x376*x377)*
(0.0001*x376 + 0.0001*x377) + 0.002*exp(-x377*x378)*(0.0001*x377 + 0.0001*
x378) + 0.002*exp(-x378*x379)*(0.0001*x378 + 0.0001*x379) + 0.002*exp(-
x379*x380)*(0.0001*x379 + 0.0001*x380) + 0.002*exp(-x380*x381)*(0.0001*
x380 + 0.0001*x381) + 0.002*exp(-x381*x382)*(0.0001*x381 + 0.0001*x382) +
0.002*exp(-x382*x383)*(0.0001*x382 + 0.0001*x383) + 0.002*exp(-x383*x384)*
(0.0001*x383 + 0.0001*x384) + 0.002*exp(-x384*x385)*(0.0001*x384 + 0.0001*
x385) + 0.002*exp(-x385*x386)*(0.0001*x385 + 0.0001*x386) + 0.002*exp(-
x386*x387)*(0.0001*x386 + 0.0001*x387) + 0.002*exp(-x387*x388)*(0.0001*
x387 + 0.0001*x388) + 0.002*exp(-x388*x389)*(0.0001*x388 + 0.0001*x389) +
0.002*exp(-x389*x390)*(0.0001*x389 + 0.0001*x390) + 0.002*exp(-x390*x391)*
(0.0001*x390 + 0.0001*x391) + 0.002*exp(-x391*x392)*(0.0001*x391 + 0.0001*
x392) + 0.002*exp(-x392*x393)*(0.0001*x392 + 0.0001*x393) + 0.002*exp(-
x393*x394)*(0.0001*x393 + 0.0001*x394) + 0.002*exp(-x394*x395)*(0.0001*
x394 + 0.0001*x395) + 0.002*exp(-x395*x396)*(0.0001*x395 + 0.0001*x396) +
0.002*exp(-x396*x397)*(0.0001*x396 + 0.0001*x397) + 0.002*exp(-x397*x398)*
(0.0001*x397 + 0.0001*x398) + 0.002*exp(-x398*x399)*(0.0001*x398 + 0.0001*
x399) + 0.002*exp(-x399*x400)*(0.0001*x399 + 0.0001*x400) + 0.002*exp(-
x400*x401)*(0.0001*x400 + 0.0001*x401) + 0.002*exp(-x401*x402)*(0.0001*
x401 + 0.0001*x402) + 0.002*exp(-x402*x403)*(0.0001*x402 + 0.0001*x403) +
0.002*exp(-x403*x404)*(0.0001*x403 + 0.0001*x404) + 0.002*exp(-x404*x405)*
(0.0001*x404 + 0.0001*x405) + 0.002*exp(-x405*x406)*(0.0001*x405 + 0.0001*
x406) + 0.002*exp(-x406*x407)*(0.0001*x406 + 0.0001*x407) + 0.002*exp(-
x407*x408)*(0.0001*x407 + 0.0001*x408) + 0.002*exp(-x408*x409)*(0.0001*
x408 + 0.0001*x409) + 0.002*exp(-x409*x410)*(0.0001*x409 + 0.0001*x410) +
0.002*exp(-x410*x411)*(0.0001*x410 + 0.0001*x411) + 0.002*exp(-x411*x412)*
(0.0001*x411 + 0.0001*x412) + 0.002*exp(-x412*x413)*(0.0001*x412 + 0.0001*
x413) + 0.002*exp(-x413*x414)*(0.0001*x413 + 0.0001*x414) + 0.002*exp(-
x414*x415)*(0.0001*x414 + 0.0001*x415) + 0.002*exp(-x415*x416)*(0.0001*
x415 + 0.0001*x416) + 0.002*exp(-x416*x417)*(0.0001*x416 + 0.0001*x417) +
0.002*exp(-x417*x418)*(0.0001*x417 + 0.0001*x418) + 0.002*exp(-x418*x419)*
(0.0001*x418 + 0.0001*x419) + 0.002*exp(-x419*x420)*(0.0001*x419 + 0.0001*
x420) + 0.002*exp(-x420*x421)*(0.0001*x420 + 0.0001*x421) + 0.002*exp(-
x421*x422)*(0.0001*x421 + 0.0001*x422) + 0.002*exp(-x422*x423)*(0.0001*
x422 + 0.0001*x423) + 0.002*exp(-x423*x424)*(0.0001*x423 + 0.0001*x424) +
0.002*exp(-x424*x425)*(0.0001*x424 + 0.0001*x425) + 0.002*exp(-x425*x426)*
(0.0001*x425 + 0.0001*x426) + 0.002*exp(-x426*x427)*(0.0001*x426 + 0.0001*
x427) + 0.002*exp(-x427*x428)*(0.0001*x427 + 0.0001*x428) + 0.002*exp(-
x428*x429)*(0.0001*x428 + 0.0001*x429) + 0.002*exp(-x429*x430)*(0.0001*
x429 + 0.0001*x430) + 0.002*exp(-x430*x431)*(0.0001*x430 + 0.0001*x431) +
0.002*exp(-x431*x432)*(0.0001*x431 + 0.0001*x432) + 0.002*exp(-x432*x433)*
(0.0001*x432 + 0.0001*x433) + 0.002*exp(-x433*x434)*(0.0001*x433 + 0.0001*
x434) + 0.002*exp(-x434*x435)*(0.0001*x434 + 0.0001*x435) + 0.002*exp(-
x435*x436)*(0.0001*x435 + 0.0001*x436) + 0.002*exp(-x436*x437)*(0.0001*
x436 + 0.0001*x437) + 0.002*exp(-x437*x438)*(0.0001*x437 + 0.0001*x438) +
0.002*exp(-x438*x439)*(0.0001*x438 + 0.0001*x439) + 0.002*exp(-x439*x440)*
(0.0001*x439 + 0.0001*x440) + 0.002*exp(-x440*x441)*(0.0001*x440 + 0.0001*
x441) + 0.002*exp(-x441*x442)*(0.0001*x441 + 0.0001*x442) + 0.002*exp(-
x442*x443)*(0.0001*x442 + 0.0001*x443) + 0.002*exp(-x443*x444)*(0.0001*
x443 + 0.0001*x444) + 0.002*exp(-x444*x445)*(0.0001*x444 + 0.0001*x445) +
0.002*exp(-x445*x446)*(0.0001*x445 + 0.0001*x446) + 0.002*exp(-x446*x447)*
(0.0001*x446 + 0.0001*x447) + 0.002*exp(-x447*x448)*(0.0001*x447 + 0.0001*
x448) + 0.002*exp(-x448*x449)*(0.0001*x448 + 0.0001*x449) + 0.002*exp(-
x449*x450)*(0.0001*x449 + 0.0001*x450) + 0.002*exp(-x450*x451)*(0.0001*
x450 + 0.0001*x451) + 0.002*exp(-x451*x452)*(0.0001*x451 + 0.0001*x452) +
0.002*exp(-x452*x453)*(0.0001*x452 + 0.0001*x453) + 0.002*exp(-x453*x454)*
(0.0001*x453 + 0.0001*x454) + 0.002*exp(-x454*x455)*(0.0001*x454 + 0.0001*
x455) + 0.002*exp(-x455*x456)*(0.0001*x455 + 0.0001*x456) + 0.002*exp(-
x456*x457)*(0.0001*x456 + 0.0001*x457) + 0.002*exp(-x457*x458)*(0.0001*
x457 + 0.0001*x458) + 0.002*exp(-x458*x459)*(0.0001*x458 + 0.0001*x459) +
0.002*exp(-x459*x460)*(0.0001*x459 + 0.0001*x460) + 0.002*exp(-x460*x461)*
(0.0001*x460 + 0.0001*x461) + 0.002*exp(-x461*x462)*(0.0001*x461 + 0.0001*
x462) + 0.002*exp(-x462*x463)*(0.0001*x462 + 0.0001*x463) + 0.002*exp(-
x463*x464)*(0.0001*x463 + 0.0001*x464) + 0.002*exp(-x464*x465)*(0.0001*
x464 + 0.0001*x465) + 0.002*exp(-x465*x466)*(0.0001*x465 + 0.0001*x466) +
0.002*exp(-x466*x467)*(0.0001*x466 + 0.0001*x467) + 0.002*exp(-x467*x468)*
(0.0001*x467 + 0.0001*x468) + 0.002*exp(-x468*x469)*(0.0001*x468 + 0.0001*
x469) + 0.002*exp(-x469*x470)*(0.0001*x469 + 0.0001*x470) + 0.002*exp(-
x470*x471)*(0.0001*x470 + 0.0001*x471) + 0.002*exp(-x471*x472)*(0.0001*
x471 + 0.0001*x472) + 0.002*exp(-x472*x473)*(0.0001*x472 + 0.0001*x473) +
0.002*exp(-x473*x474)*(0.0001*x473 + 0.0001*x474) + 0.002*exp(-x474*x475)*
(0.0001*x474 + 0.0001*x475) + 0.002*exp(-x475*x476)*(0.0001*x475 + 0.0001*
x476) + 0.002*exp(-x476*x477)*(0.0001*x476 + 0.0001*x477) + 0.002*exp(-
x477*x478)*(0.0001*x477 + 0.0001*x478) + 0.002*exp(-x478*x479)*(0.0001*
x478 + 0.0001*x479) + 0.002*exp(-x479*x480)*(0.0001*x479 + 0.0001*x480) +
0.002*exp(-x480*x481)*(0.0001*x480 + 0.0001*x481) + 0.002*exp(-x481*x482)*
(0.0001*x481 + 0.0001*x482) + 0.002*exp(-x482*x483)*(0.0001*x482 + 0.0001*
x483) + 0.002*exp(-x483*x484)*(0.0001*x483 + 0.0001*x484) + 0.002*exp(-
x484*x485)*(0.0001*x484 + 0.0001*x485) + 0.002*exp(-x485*x486)*(0.0001*
x485 + 0.0001*x486) + 0.002*exp(-x486*x487)*(0.0001*x486 + 0.0001*x487) +
0.002*exp(-x487*x488)*(0.0001*x487 + 0.0001*x488) + 0.002*exp(-x488*x489)*
(0.0001*x488 + 0.0001*x489) + 0.002*exp(-x489*x490)*(0.0001*x489 + 0.0001*
x490) + 0.002*exp(-x490*x491)*(0.0001*x490 + 0.0001*x491) + 0.002*exp(-
x491*x492)*(0.0001*x491 + 0.0001*x492) + 0.002*exp(-x492*x493)*(0.0001*
x492 + 0.0001*x493) + 0.002*exp(-x493*x494)*(0.0001*x493 + 0.0001*x494) +
0.002*exp(-x494*x495)*(0.0001*x494 + 0.0001*x495) + 0.002*exp(-x495*x496)*
(0.0001*x495 + 0.0001*x496) + 0.002*exp(-x496*x497)*(0.0001*x496 + 0.0001*
x497) + 0.002*exp(-x497*x498)*(0.0001*x497 + 0.0001*x498) + 0.002*exp(-
x498*x499)*(0.0001*x498 + 0.0001*x499) + 0.002*exp(-x499*x500)*(0.0001*
x499 + 0.0001*x500) + 100*sqr((-1) + 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)) + objvar =E= 0;
* set non default bounds
x1.up = 5;
x2.up = 5;
x3.up = 5;
x4.up = 5;
x5.up = 5;
x6.up = 5;
x7.up = 5;
x8.up = 5;
x9.up = 5;
x10.up = 5;
x11.up = 5;
x12.up = 5;
x13.up = 5;
x14.up = 5;
x15.up = 5;
x16.up = 5;
x17.up = 5;
x18.up = 5;
x19.up = 5;
x20.up = 5;
x21.up = 5;
x22.up = 5;
x23.up = 5;
x24.up = 5;
x25.up = 5;
x26.up = 5;
x27.up = 5;
x28.up = 5;
x29.up = 5;
x30.up = 5;
x31.up = 5;
x32.up = 5;
x33.up = 5;
x34.up = 5;
x35.up = 5;
x36.up = 5;
x37.up = 5;
x38.up = 5;
x39.up = 5;
x40.up = 5;
x41.up = 5;
x42.up = 5;
x43.up = 5;
x44.up = 5;
x45.up = 5;
x46.up = 5;
x47.up = 5;
x48.up = 5;
x49.up = 5;
x50.up = 5;
x51.up = 5;
x52.up = 5;
x53.up = 5;
x54.up = 5;
x55.up = 5;
x56.up = 5;
x57.up = 5;
x58.up = 5;
x59.up = 5;
x60.up = 5;
x61.up = 5;
x62.up = 5;
x63.up = 5;
x64.up = 5;
x65.up = 5;
x66.up = 5;
x67.up = 5;
x68.up = 5;
x69.up = 5;
x70.up = 5;
x71.up = 5;
x72.up = 5;
x73.up = 5;
x74.up = 5;
x75.up = 5;
x76.up = 5;
x77.up = 5;
x78.up = 5;
x79.up = 5;
x80.up = 5;
x81.up = 5;
x82.up = 5;
x83.up = 5;
x84.up = 5;
x85.up = 5;
x86.up = 5;
x87.up = 5;
x88.up = 5;
x89.up = 5;
x90.up = 5;
x91.up = 5;
x92.up = 5;
x93.up = 5;
x94.up = 5;
x95.up = 5;
x96.up = 5;
x97.up = 5;
x98.up = 5;
x99.up = 5;
x100.up = 5;
x101.up = 5;
x102.up = 5;
x103.up = 5;
x104.up = 5;
x105.up = 5;
x106.up = 5;
x107.up = 5;
x108.up = 5;
x109.up = 5;
x110.up = 5;
x111.up = 5;
x112.up = 5;
x113.up = 5;
x114.up = 5;
x115.up = 5;
x116.up = 5;
x117.up = 5;
x118.up = 5;
x119.up = 5;
x120.up = 5;
x121.up = 5;
x122.up = 5;
x123.up = 5;
x124.up = 5;
x125.up = 5;
x126.up = 5;
x127.up = 5;
x128.up = 5;
x129.up = 5;
x130.up = 5;
x131.up = 5;
x132.up = 5;
x133.up = 5;
x134.up = 5;
x135.up = 5;
x136.up = 5;
x137.up = 5;
x138.up = 5;
x139.up = 5;
x140.up = 5;
x141.up = 5;
x142.up = 5;
x143.up = 5;
x144.up = 5;
x145.up = 5;
x146.up = 5;
x147.up = 5;
x148.up = 5;
x149.up = 5;
x150.up = 5;
x151.up = 5;
x152.up = 5;
x153.up = 5;
x154.up = 5;
x155.up = 5;
x156.up = 5;
x157.up = 5;
x158.up = 5;
x159.up = 5;
x160.up = 5;
x161.up = 5;
x162.up = 5;
x163.up = 5;
x164.up = 5;
x165.up = 5;
x166.up = 5;
x167.up = 5;
x168.up = 5;
x169.up = 5;
x170.up = 5;
x171.up = 5;
x172.up = 5;
x173.up = 5;
x174.up = 5;
x175.up = 5;
x176.up = 5;
x177.up = 5;
x178.up = 5;
x179.up = 5;
x180.up = 5;
x181.up = 5;
x182.up = 5;
x183.up = 5;
x184.up = 5;
x185.up = 5;
x186.up = 5;
x187.up = 5;
x188.up = 5;
x189.up = 5;
x190.up = 5;
x191.up = 5;
x192.up = 5;
x193.up = 5;
x194.up = 5;
x195.up = 5;
x196.up = 5;
x197.up = 5;
x198.up = 5;
x199.up = 5;
x200.up = 5;
x201.up = 5;
x202.up = 5;
x203.up = 5;
x204.up = 5;
x205.up = 5;
x206.up = 5;
x207.up = 5;
x208.up = 5;
x209.up = 5;
x210.up = 5;
x211.up = 5;
x212.up = 5;
x213.up = 5;
x214.up = 5;
x215.up = 5;
x216.up = 5;
x217.up = 5;
x218.up = 5;
x219.up = 5;
x220.up = 5;
x221.up = 5;
x222.up = 5;
x223.up = 5;
x224.up = 5;
x225.up = 5;
x226.up = 5;
x227.up = 5;
x228.up = 5;
x229.up = 5;
x230.up = 5;
x231.up = 5;
x232.up = 5;
x233.up = 5;
x234.up = 5;
x235.up = 5;
x236.up = 5;
x237.up = 5;
x238.up = 5;
x239.up = 5;
x240.up = 5;
x241.up = 5;
x242.up = 5;
x243.up = 5;
x244.up = 5;
x245.up = 5;
x246.up = 5;
x247.up = 5;
x248.up = 5;
x249.up = 5;
x250.up = 5;
x251.up = 5;
x252.up = 5;
x253.up = 5;
x254.up = 5;
x255.up = 5;
x256.up = 5;
x257.up = 5;
x258.up = 5;
x259.up = 5;
x260.up = 5;
x261.up = 5;
x262.up = 5;
x263.up = 5;
x264.up = 5;
x265.up = 5;
x266.up = 5;
x267.up = 5;
x268.up = 5;
x269.up = 5;
x270.up = 5;
x271.up = 5;
x272.up = 5;
x273.up = 5;
x274.up = 5;
x275.up = 5;
x276.up = 5;
x277.up = 5;
x278.up = 5;
x279.up = 5;
x280.up = 5;
x281.up = 5;
x282.up = 5;
x283.up = 5;
x284.up = 5;
x285.up = 5;
x286.up = 5;
x287.up = 5;
x288.up = 5;
x289.up = 5;
x290.up = 5;
x291.up = 5;
x292.up = 5;
x293.up = 5;
x294.up = 5;
x295.up = 5;
x296.up = 5;
x297.up = 5;
x298.up = 5;
x299.up = 5;
x300.up = 5;
x301.up = 5;
x302.up = 5;
x303.up = 5;
x304.up = 5;
x305.up = 5;
x306.up = 5;
x307.up = 5;
x308.up = 5;
x309.up = 5;
x310.up = 5;
x311.up = 5;
x312.up = 5;
x313.up = 5;
x314.up = 5;
x315.up = 5;
x316.up = 5;
x317.up = 5;
x318.up = 5;
x319.up = 5;
x320.up = 5;
x321.up = 5;
x322.up = 5;
x323.up = 5;
x324.up = 5;
x325.up = 5;
x326.up = 5;
x327.up = 5;
x328.up = 5;
x329.up = 5;
x330.up = 5;
x331.up = 5;
x332.up = 5;
x333.up = 5;
x334.up = 5;
x335.up = 5;
x336.up = 5;
x337.up = 5;
x338.up = 5;
x339.up = 5;
x340.up = 5;
x341.up = 5;
x342.up = 5;
x343.up = 5;
x344.up = 5;
x345.up = 5;
x346.up = 5;
x347.up = 5;
x348.up = 5;
x349.up = 5;
x350.up = 5;
x351.up = 5;
x352.up = 5;
x353.up = 5;
x354.up = 5;
x355.up = 5;
x356.up = 5;
x357.up = 5;
x358.up = 5;
x359.up = 5;
x360.up = 5;
x361.up = 5;
x362.up = 5;
x363.up = 5;
x364.up = 5;
x365.up = 5;
x366.up = 5;
x367.up = 5;
x368.up = 5;
x369.up = 5;
x370.up = 5;
x371.up = 5;
x372.up = 5;
x373.up = 5;
x374.up = 5;
x375.up = 5;
x376.up = 5;
x377.up = 5;
x378.up = 5;
x379.up = 5;
x380.up = 5;
x381.up = 5;
x382.up = 5;
x383.up = 5;
x384.up = 5;
x385.up = 5;
x386.up = 5;
x387.up = 5;
x388.up = 5;
x389.up = 5;
x390.up = 5;
x391.up = 5;
x392.up = 5;
x393.up = 5;
x394.up = 5;
x395.up = 5;
x396.up = 5;
x397.up = 5;
x398.up = 5;
x399.up = 5;
x400.up = 5;
x401.up = 5;
x402.up = 5;
x403.up = 5;
x404.up = 5;
x405.up = 5;
x406.up = 5;
x407.up = 5;
x408.up = 5;
x409.up = 5;
x410.up = 5;
x411.up = 5;
x412.up = 5;
x413.up = 5;
x414.up = 5;
x415.up = 5;
x416.up = 5;
x417.up = 5;
x418.up = 5;
x419.up = 5;
x420.up = 5;
x421.up = 5;
x422.up = 5;
x423.up = 5;
x424.up = 5;
x425.up = 5;
x426.up = 5;
x427.up = 5;
x428.up = 5;
x429.up = 5;
x430.up = 5;
x431.up = 5;
x432.up = 5;
x433.up = 5;
x434.up = 5;
x435.up = 5;
x436.up = 5;
x437.up = 5;
x438.up = 5;
x439.up = 5;
x440.up = 5;
x441.up = 5;
x442.up = 5;
x443.up = 5;
x444.up = 5;
x445.up = 5;
x446.up = 5;
x447.up = 5;
x448.up = 5;
x449.up = 5;
x450.up = 5;
x451.up = 5;
x452.up = 5;
x453.up = 5;
x454.up = 5;
x455.up = 5;
x456.up = 5;
x457.up = 5;
x458.up = 5;
x459.up = 5;
x460.up = 5;
x461.up = 5;
x462.up = 5;
x463.up = 5;
x464.up = 5;
x465.up = 5;
x466.up = 5;
x467.up = 5;
x468.up = 5;
x469.up = 5;
x470.up = 5;
x471.up = 5;
x472.up = 5;
x473.up = 5;
x474.up = 5;
x475.up = 5;
x476.up = 5;
x477.up = 5;
x478.up = 5;
x479.up = 5;
x480.up = 5;
x481.up = 5;
x482.up = 5;
x483.up = 5;
x484.up = 5;
x485.up = 5;
x486.up = 5;
x487.up = 5;
x488.up = 5;
x489.up = 5;
x490.up = 5;
x491.up = 5;
x492.up = 5;
x493.up = 5;
x494.up = 5;
x495.up = 5;
x496.up = 5;
x497.up = 5;
x498.up = 5;
x499.up = 5;
x500.up = 5;
* set non default levels
x1.l = 0.5;
x2.l = 0.5;
x3.l = 0.5;
x4.l = 0.5;
x5.l = 0.5;
x6.l = 0.5;
x7.l = 0.5;
x8.l = 0.5;
x9.l = 0.5;
x10.l = 0.5;
x11.l = 0.5;
x12.l = 0.5;
x13.l = 0.5;
x14.l = 0.5;
x15.l = 0.5;
x16.l = 0.5;
x17.l = 0.5;
x18.l = 0.5;
x19.l = 0.5;
x20.l = 0.5;
x21.l = 0.5;
x22.l = 0.5;
x23.l = 0.5;
x24.l = 0.5;
x25.l = 0.5;
x26.l = 0.5;
x27.l = 0.5;
x28.l = 0.5;
x29.l = 0.5;
x30.l = 0.5;
x31.l = 0.5;
x32.l = 0.5;
x33.l = 0.5;
x34.l = 0.5;
x35.l = 0.5;
x36.l = 0.5;
x37.l = 0.5;
x38.l = 0.5;
x39.l = 0.5;
x40.l = 0.5;
x41.l = 0.5;
x42.l = 0.5;
x43.l = 0.5;
x44.l = 0.5;
x45.l = 0.5;
x46.l = 0.5;
x47.l = 0.5;
x48.l = 0.5;
x49.l = 0.5;
x50.l = 0.5;
x51.l = 0.5;
x52.l = 0.5;
x53.l = 0.5;
x54.l = 0.5;
x55.l = 0.5;
x56.l = 0.5;
x57.l = 0.5;
x58.l = 0.5;
x59.l = 0.5;
x60.l = 0.5;
x61.l = 0.5;
x62.l = 0.5;
x63.l = 0.5;
x64.l = 0.5;
x65.l = 0.5;
x66.l = 0.5;
x67.l = 0.5;
x68.l = 0.5;
x69.l = 0.5;
x70.l = 0.5;
x71.l = 0.5;
x72.l = 0.5;
x73.l = 0.5;
x74.l = 0.5;
x75.l = 0.5;
x76.l = 0.5;
x77.l = 0.5;
x78.l = 0.5;
x79.l = 0.5;
x80.l = 0.5;
x81.l = 0.5;
x82.l = 0.5;
x83.l = 0.5;
x84.l = 0.5;
x85.l = 0.5;
x86.l = 0.5;
x87.l = 0.5;
x88.l = 0.5;
x89.l = 0.5;
x90.l = 0.5;
x91.l = 0.5;
x92.l = 0.5;
x93.l = 0.5;
x94.l = 0.5;
x95.l = 0.5;
x96.l = 0.5;
x97.l = 0.5;
x98.l = 0.5;
x99.l = 0.5;
x100.l = 0.5;
x101.l = 0.5;
x102.l = 0.5;
x103.l = 0.5;
x104.l = 0.5;
x105.l = 0.5;
x106.l = 0.5;
x107.l = 0.5;
x108.l = 0.5;
x109.l = 0.5;
x110.l = 0.5;
x111.l = 0.5;
x112.l = 0.5;
x113.l = 0.5;
x114.l = 0.5;
x115.l = 0.5;
x116.l = 0.5;
x117.l = 0.5;
x118.l = 0.5;
x119.l = 0.5;
x120.l = 0.5;
x121.l = 0.5;
x122.l = 0.5;
x123.l = 0.5;
x124.l = 0.5;
x125.l = 0.5;
x126.l = 0.5;
x127.l = 0.5;
x128.l = 0.5;
x129.l = 0.5;
x130.l = 0.5;
x131.l = 0.5;
x132.l = 0.5;
x133.l = 0.5;
x134.l = 0.5;
x135.l = 0.5;
x136.l = 0.5;
x137.l = 0.5;
x138.l = 0.5;
x139.l = 0.5;
x140.l = 0.5;
x141.l = 0.5;
x142.l = 0.5;
x143.l = 0.5;
x144.l = 0.5;
x145.l = 0.5;
x146.l = 0.5;
x147.l = 0.5;
x148.l = 0.5;
x149.l = 0.5;
x150.l = 0.5;
x151.l = 0.5;
x152.l = 0.5;
x153.l = 0.5;
x154.l = 0.5;
x155.l = 0.5;
x156.l = 0.5;
x157.l = 0.5;
x158.l = 0.5;
x159.l = 0.5;
x160.l = 0.5;
x161.l = 0.5;
x162.l = 0.5;
x163.l = 0.5;
x164.l = 0.5;
x165.l = 0.5;
x166.l = 0.5;
x167.l = 0.5;
x168.l = 0.5;
x169.l = 0.5;
x170.l = 0.5;
x171.l = 0.5;
x172.l = 0.5;
x173.l = 0.5;
x174.l = 0.5;
x175.l = 0.5;
x176.l = 0.5;
x177.l = 0.5;
x178.l = 0.5;
x179.l = 0.5;
x180.l = 0.5;
x181.l = 0.5;
x182.l = 0.5;
x183.l = 0.5;
x184.l = 0.5;
x185.l = 0.5;
x186.l = 0.5;
x187.l = 0.5;
x188.l = 0.5;
x189.l = 0.5;
x190.l = 0.5;
x191.l = 0.5;
x192.l = 0.5;
x193.l = 0.5;
x194.l = 0.5;
x195.l = 0.5;
x196.l = 0.5;
x197.l = 0.5;
x198.l = 0.5;
x199.l = 0.5;
x200.l = 0.5;
x201.l = 0.5;
x202.l = 0.5;
x203.l = 0.5;
x204.l = 0.5;
x205.l = 0.5;
x206.l = 0.5;
x207.l = 0.5;
x208.l = 0.5;
x209.l = 0.5;
x210.l = 0.5;
x211.l = 0.5;
x212.l = 0.5;
x213.l = 0.5;
x214.l = 0.5;
x215.l = 0.5;
x216.l = 0.5;
x217.l = 0.5;
x218.l = 0.5;
x219.l = 0.5;
x220.l = 0.5;
x221.l = 0.5;
x222.l = 0.5;
x223.l = 0.5;
x224.l = 0.5;
x225.l = 0.5;
x226.l = 0.5;
x227.l = 0.5;
x228.l = 0.5;
x229.l = 0.5;
x230.l = 0.5;
x231.l = 0.5;
x232.l = 0.5;
x233.l = 0.5;
x234.l = 0.5;
x235.l = 0.5;
x236.l = 0.5;
x237.l = 0.5;
x238.l = 0.5;
x239.l = 0.5;
x240.l = 0.5;
x241.l = 0.5;
x242.l = 0.5;
x243.l = 0.5;
x244.l = 0.5;
x245.l = 0.5;
x246.l = 0.5;
x247.l = 0.5;
x248.l = 0.5;
x249.l = 0.5;
x250.l = 0.5;
x251.l = 0.5;
x252.l = 0.5;
x253.l = 0.5;
x254.l = 0.5;
x255.l = 0.5;
x256.l = 0.5;
x257.l = 0.5;
x258.l = 0.5;
x259.l = 0.5;
x260.l = 0.5;
x261.l = 0.5;
x262.l = 0.5;
x263.l = 0.5;
x264.l = 0.5;
x265.l = 0.5;
x266.l = 0.5;
x267.l = 0.5;
x268.l = 0.5;
x269.l = 0.5;
x270.l = 0.5;
x271.l = 0.5;
x272.l = 0.5;
x273.l = 0.5;
x274.l = 0.5;
x275.l = 0.5;
x276.l = 0.5;
x277.l = 0.5;
x278.l = 0.5;
x279.l = 0.5;
x280.l = 0.5;
x281.l = 0.5;
x282.l = 0.5;
x283.l = 0.5;
x284.l = 0.5;
x285.l = 0.5;
x286.l = 0.5;
x287.l = 0.5;
x288.l = 0.5;
x289.l = 0.5;
x290.l = 0.5;
x291.l = 0.5;
x292.l = 0.5;
x293.l = 0.5;
x294.l = 0.5;
x295.l = 0.5;
x296.l = 0.5;
x297.l = 0.5;
x298.l = 0.5;
x299.l = 0.5;
x300.l = 0.5;
x301.l = 0.5;
x302.l = 0.5;
x303.l = 0.5;
x304.l = 0.5;
x305.l = 0.5;
x306.l = 0.5;
x307.l = 0.5;
x308.l = 0.5;
x309.l = 0.5;
x310.l = 0.5;
x311.l = 0.5;
x312.l = 0.5;
x313.l = 0.5;
x314.l = 0.5;
x315.l = 0.5;
x316.l = 0.5;
x317.l = 0.5;
x318.l = 0.5;
x319.l = 0.5;
x320.l = 0.5;
x321.l = 0.5;
x322.l = 0.5;
x323.l = 0.5;
x324.l = 0.5;
x325.l = 0.5;
x326.l = 0.5;
x327.l = 0.5;
x328.l = 0.5;
x329.l = 0.5;
x330.l = 0.5;
x331.l = 0.5;
x332.l = 0.5;
x333.l = 0.5;
x334.l = 0.5;
x335.l = 0.5;
x336.l = 0.5;
x337.l = 0.5;
x338.l = 0.5;
x339.l = 0.5;
x340.l = 0.5;
x341.l = 0.5;
x342.l = 0.5;
x343.l = 0.5;
x344.l = 0.5;
x345.l = 0.5;
x346.l = 0.5;
x347.l = 0.5;
x348.l = 0.5;
x349.l = 0.5;
x350.l = 0.5;
x351.l = 0.5;
x352.l = 0.5;
x353.l = 0.5;
x354.l = 0.5;
x355.l = 0.5;
x356.l = 0.5;
x357.l = 0.5;
x358.l = 0.5;
x359.l = 0.5;
x360.l = 0.5;
x361.l = 0.5;
x362.l = 0.5;
x363.l = 0.5;
x364.l = 0.5;
x365.l = 0.5;
x366.l = 0.5;
x367.l = 0.5;
x368.l = 0.5;
x369.l = 0.5;
x370.l = 0.5;
x371.l = 0.5;
x372.l = 0.5;
x373.l = 0.5;
x374.l = 0.5;
x375.l = 0.5;
x376.l = 0.5;
x377.l = 0.5;
x378.l = 0.5;
x379.l = 0.5;
x380.l = 0.5;
x381.l = 0.5;
x382.l = 0.5;
x383.l = 0.5;
x384.l = 0.5;
x385.l = 0.5;
x386.l = 0.5;
x387.l = 0.5;
x388.l = 0.5;
x389.l = 0.5;
x390.l = 0.5;
x391.l = 0.5;
x392.l = 0.5;
x393.l = 0.5;
x394.l = 0.5;
x395.l = 0.5;
x396.l = 0.5;
x397.l = 0.5;
x398.l = 0.5;
x399.l = 0.5;
x400.l = 0.5;
x401.l = 0.5;
x402.l = 0.5;
x403.l = 0.5;
x404.l = 0.5;
x405.l = 0.5;
x406.l = 0.5;
x407.l = 0.5;
x408.l = 0.5;
x409.l = 0.5;
x410.l = 0.5;
x411.l = 0.5;
x412.l = 0.5;
x413.l = 0.5;
x414.l = 0.5;
x415.l = 0.5;
x416.l = 0.5;
x417.l = 0.5;
x418.l = 0.5;
x419.l = 0.5;
x420.l = 0.5;
x421.l = 0.5;
x422.l = 0.5;
x423.l = 0.5;
x424.l = 0.5;
x425.l = 0.5;
x426.l = 0.5;
x427.l = 0.5;
x428.l = 0.5;
x429.l = 0.5;
x430.l = 0.5;
x431.l = 0.5;
x432.l = 0.5;
x433.l = 0.5;
x434.l = 0.5;
x435.l = 0.5;
x436.l = 0.5;
x437.l = 0.5;
x438.l = 0.5;
x439.l = 0.5;
x440.l = 0.5;
x441.l = 0.5;
x442.l = 0.5;
x443.l = 0.5;
x444.l = 0.5;
x445.l = 0.5;
x446.l = 0.5;
x447.l = 0.5;
x448.l = 0.5;
x449.l = 0.5;
x450.l = 0.5;
x451.l = 0.5;
x452.l = 0.5;
x453.l = 0.5;
x454.l = 0.5;
x455.l = 0.5;
x456.l = 0.5;
x457.l = 0.5;
x458.l = 0.5;
x459.l = 0.5;
x460.l = 0.5;
x461.l = 0.5;
x462.l = 0.5;
x463.l = 0.5;
x464.l = 0.5;
x465.l = 0.5;
x466.l = 0.5;
x467.l = 0.5;
x468.l = 0.5;
x469.l = 0.5;
x470.l = 0.5;
x471.l = 0.5;
x472.l = 0.5;
x473.l = 0.5;
x474.l = 0.5;
x475.l = 0.5;
x476.l = 0.5;
x477.l = 0.5;
x478.l = 0.5;
x479.l = 0.5;
x480.l = 0.5;
x481.l = 0.5;
x482.l = 0.5;
x483.l = 0.5;
x484.l = 0.5;
x485.l = 0.5;
x486.l = 0.5;
x487.l = 0.5;
x488.l = 0.5;
x489.l = 0.5;
x490.l = 0.5;
x491.l = 0.5;
x492.l = 0.5;
x493.l = 0.5;
x494.l = 0.5;
x495.l = 0.5;
x496.l = 0.5;
x497.l = 0.5;
x498.l = 0.5;
x499.l = 0.5;
x500.l = 0.5;
* 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;