* NLP written by GAMS Convert at 10/06/06 11:33:48 * * Equation counts * Total E G L N X C * 151 151 0 0 0 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 203 203 0 0 0 0 0 0 * FX 1 1 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 602 152 450 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,objvar; Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19 ,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36 ,e37,e38,e39,e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,e50,e51,e52,e53 ,e54,e55,e56,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; e1.. 0.00437*cos(x5)/((0.0162079 + 0.486237*x153*x153)*x4*x4) + x1 - x3 =E= 0; e2.. (-1/x4) - cos(x5)/((0.0162079 + 0.486237*x153*x153)*x4*x4*x4) =E= 0; e3.. 0.00437*x153/((0.01 + 0.3*x153*x153)*x4*x4) - 0.00437*cos(x5)/((0.0162079 + 0.486237*x153*x153)*x4*x4*x4*x4) - 0.1*x5 + 0.1*x2 =E= 0; e4.. 0.00437*cos(x8)/((0.0162079 + 0.486237*x154*x154)*x7*x7) + x3 - x6 =E= 0; e5.. (-1/x7) - cos(x8)/((0.0162079 + 0.486237*x154*x154)*x7*x7*x7) =E= 0; e6.. 0.00437*x154/((0.01 + 0.3*x154*x154)*x7*x7) - 0.00437*cos(x8)/((0.0162079 + 0.486237*x154*x154)*x7*x7*x7*x7) - 0.1*x8 + 0.1*x5 =E= 0; e7.. 0.00437*cos(x11)/((0.0162079 + 0.486237*x155*x155)*x10*x10) + x6 - x9 =E= 0; e8.. (-1/x10) - cos(x11)/((0.0162079 + 0.486237*x155*x155)*x10*x10*x10) =E= 0; e9.. 0.00437*x155/((0.01 + 0.3*x155*x155)*x10*x10) - 0.00437*cos(x11)/(( 0.0162079 + 0.486237*x155*x155)*x10*x10*x10*x10) - 0.1*x11 + 0.1*x8 =E= 0; e10.. 0.00437*cos(x14)/((0.0162079 + 0.486237*x156*x156)*x13*x13) + x9 - x12 =E= 0; e11.. (-1/x13) - cos(x14)/((0.0162079 + 0.486237*x156*x156)*x13*x13*x13) =E= 0; e12.. 0.00437*x156/((0.01 + 0.3*x156*x156)*x13*x13) - 0.00437*cos(x14)/(( 0.0162079 + 0.486237*x156*x156)*x13*x13*x13*x13) - 0.1*x14 + 0.1*x11 =E= 0; e13.. 0.00437*cos(x17)/((0.0162079 + 0.486237*x157*x157)*x16*x16) + x12 - x15 =E= 0; e14.. (-1/x16) - cos(x17)/((0.0162079 + 0.486237*x157*x157)*x16*x16*x16) =E= 0; e15.. 0.00437*x157/((0.01 + 0.3*x157*x157)*x16*x16) - 0.00437*cos(x17)/(( 0.0162079 + 0.486237*x157*x157)*x16*x16*x16*x16) - 0.1*x17 + 0.1*x14 =E= 0; e16.. 0.00437*cos(x20)/((0.0162079 + 0.486237*x158*x158)*x19*x19) + x15 - x18 =E= 0; e17.. (-1/x19) - cos(x20)/((0.0162079 + 0.486237*x158*x158)*x19*x19*x19) =E= 0; e18.. 0.00437*x158/((0.01 + 0.3*x158*x158)*x19*x19) - 0.00437*cos(x20)/(( 0.0162079 + 0.486237*x158*x158)*x19*x19*x19*x19) - 0.1*x20 + 0.1*x17 =E= 0; e19.. 0.00437*cos(x23)/((0.0162079 + 0.486237*x159*x159)*x22*x22) + x18 - x21 =E= 0; e20.. (-1/x22) - cos(x23)/((0.0162079 + 0.486237*x159*x159)*x22*x22*x22) =E= 0; e21.. 0.00437*x159/((0.01 + 0.3*x159*x159)*x22*x22) - 0.00437*cos(x23)/(( 0.0162079 + 0.486237*x159*x159)*x22*x22*x22*x22) - 0.1*x23 + 0.1*x20 =E= 0; e22.. 0.00437*cos(x26)/((0.0162079 + 0.486237*x160*x160)*x25*x25) + x21 - x24 =E= 0; e23.. (-1/x25) - cos(x26)/((0.0162079 + 0.486237*x160*x160)*x25*x25*x25) =E= 0; e24.. 0.00437*x160/((0.01 + 0.3*x160*x160)*x25*x25) - 0.00437*cos(x26)/(( 0.0162079 + 0.486237*x160*x160)*x25*x25*x25*x25) - 0.1*x26 + 0.1*x23 =E= 0; e25.. 0.00437*cos(x29)/((0.0162079 + 0.486237*x161*x161)*x28*x28) + x24 - x27 =E= 0; e26.. (-1/x28) - cos(x29)/((0.0162079 + 0.486237*x161*x161)*x28*x28*x28) =E= 0; e27.. 0.00437*x161/((0.01 + 0.3*x161*x161)*x28*x28) - 0.00437*cos(x29)/(( 0.0162079 + 0.486237*x161*x161)*x28*x28*x28*x28) - 0.1*x29 + 0.1*x26 =E= 0; e28.. 0.00437*cos(x32)/((0.0162079 + 0.486237*x162*x162)*x31*x31) + x27 - x30 =E= 0; e29.. (-1/x31) - cos(x32)/((0.0162079 + 0.486237*x162*x162)*x31*x31*x31) =E= 0; e30.. 0.00437*x162/((0.01 + 0.3*x162*x162)*x31*x31) - 0.00437*cos(x32)/(( 0.0162079 + 0.486237*x162*x162)*x31*x31*x31*x31) - 0.1*x32 + 0.1*x29 =E= 0; e31.. 0.00437*cos(x35)/((0.0162079 + 0.486237*x163*x163)*x34*x34) + x30 - x33 =E= 0; e32.. (-1/x34) - cos(x35)/((0.0162079 + 0.486237*x163*x163)*x34*x34*x34) =E= 0; e33.. 0.00437*x163/((0.01 + 0.3*x163*x163)*x34*x34) - 0.00437*cos(x35)/(( 0.0162079 + 0.486237*x163*x163)*x34*x34*x34*x34) - 0.1*x35 + 0.1*x32 =E= 0; e34.. 0.00437*cos(x38)/((0.0162079 + 0.486237*x164*x164)*x37*x37) + x33 - x36 =E= 0; e35.. (-1/x37) - cos(x38)/((0.0162079 + 0.486237*x164*x164)*x37*x37*x37) =E= 0; e36.. 0.00437*x164/((0.01 + 0.3*x164*x164)*x37*x37) - 0.00437*cos(x38)/(( 0.0162079 + 0.486237*x164*x164)*x37*x37*x37*x37) - 0.1*x38 + 0.1*x35 =E= 0; e37.. 0.00437*cos(x41)/((0.0162079 + 0.486237*x165*x165)*x40*x40) + x36 - x39 =E= 0; e38.. (-1/x40) - cos(x41)/((0.0162079 + 0.486237*x165*x165)*x40*x40*x40) =E= 0; e39.. 0.00437*x165/((0.01 + 0.3*x165*x165)*x40*x40) - 0.00437*cos(x41)/(( 0.0162079 + 0.486237*x165*x165)*x40*x40*x40*x40) - 0.1*x41 + 0.1*x38 =E= 0; e40.. 0.00437*cos(x44)/((0.0162079 + 0.486237*x166*x166)*x43*x43) + x39 - x42 =E= 0; e41.. (-1/x43) - cos(x44)/((0.0162079 + 0.486237*x166*x166)*x43*x43*x43) =E= 0; e42.. 0.00437*x166/((0.01 + 0.3*x166*x166)*x43*x43) - 0.00437*cos(x44)/(( 0.0162079 + 0.486237*x166*x166)*x43*x43*x43*x43) - 0.1*x44 + 0.1*x41 =E= 0; e43.. 0.00437*cos(x47)/((0.0162079 + 0.486237*x167*x167)*x46*x46) + x42 - x45 =E= 0; e44.. (-1/x46) - cos(x47)/((0.0162079 + 0.486237*x167*x167)*x46*x46*x46) =E= 0; e45.. 0.00437*x167/((0.01 + 0.3*x167*x167)*x46*x46) - 0.00437*cos(x47)/(( 0.0162079 + 0.486237*x167*x167)*x46*x46*x46*x46) - 0.1*x47 + 0.1*x44 =E= 0; e46.. 0.00437*cos(x50)/((0.0162079 + 0.486237*x168*x168)*x49*x49) + x45 - x48 =E= 0; e47.. (-1/x49) - cos(x50)/((0.0162079 + 0.486237*x168*x168)*x49*x49*x49) =E= 0; e48.. 0.00437*x168/((0.01 + 0.3*x168*x168)*x49*x49) - 0.00437*cos(x50)/(( 0.0162079 + 0.486237*x168*x168)*x49*x49*x49*x49) - 0.1*x50 + 0.1*x47 =E= 0; e49.. 0.00437*cos(x53)/((0.0162079 + 0.486237*x169*x169)*x52*x52) + x48 - x51 =E= 0; e50.. (-1/x52) - cos(x53)/((0.0162079 + 0.486237*x169*x169)*x52*x52*x52) =E= 0; e51.. 0.00437*x169/((0.01 + 0.3*x169*x169)*x52*x52) - 0.00437*cos(x53)/(( 0.0162079 + 0.486237*x169*x169)*x52*x52*x52*x52) - 0.1*x53 + 0.1*x50 =E= 0; e52.. 0.00437*cos(x56)/((0.0162079 + 0.486237*x170*x170)*x55*x55) + x51 - x54 =E= 0; e53.. (-1/x55) - cos(x56)/((0.0162079 + 0.486237*x170*x170)*x55*x55*x55) =E= 0; e54.. 0.00437*x170/((0.01 + 0.3*x170*x170)*x55*x55) - 0.00437*cos(x56)/(( 0.0162079 + 0.486237*x170*x170)*x55*x55*x55*x55) - 0.1*x56 + 0.1*x53 =E= 0; e55.. 0.00437*cos(x59)/((0.0162079 + 0.486237*x171*x171)*x58*x58) + x54 - x57 =E= 0; e56.. (-1/x58) - cos(x59)/((0.0162079 + 0.486237*x171*x171)*x58*x58*x58) =E= 0; e57.. 0.00437*x171/((0.01 + 0.3*x171*x171)*x58*x58) - 0.00437*cos(x59)/(( 0.0162079 + 0.486237*x171*x171)*x58*x58*x58*x58) - 0.1*x59 + 0.1*x56 =E= 0; e58.. 0.00437*cos(x62)/((0.0162079 + 0.486237*x172*x172)*x61*x61) + x57 - x60 =E= 0; e59.. (-1/x61) - cos(x62)/((0.0162079 + 0.486237*x172*x172)*x61*x61*x61) =E= 0; e60.. 0.00437*x172/((0.01 + 0.3*x172*x172)*x61*x61) - 0.00437*cos(x62)/(( 0.0162079 + 0.486237*x172*x172)*x61*x61*x61*x61) - 0.1*x62 + 0.1*x59 =E= 0; e61.. 0.00437*cos(x65)/((0.0162079 + 0.486237*x173*x173)*x64*x64) + x60 - x63 =E= 0; e62.. (-1/x64) - cos(x65)/((0.0162079 + 0.486237*x173*x173)*x64*x64*x64) =E= 0; e63.. 0.00437*x173/((0.01 + 0.3*x173*x173)*x64*x64) - 0.00437*cos(x65)/(( 0.0162079 + 0.486237*x173*x173)*x64*x64*x64*x64) - 0.1*x65 + 0.1*x62 =E= 0; e64.. 0.00437*cos(x68)/((0.0162079 + 0.486237*x174*x174)*x67*x67) + x63 - x66 =E= 0; e65.. (-1/x67) - cos(x68)/((0.0162079 + 0.486237*x174*x174)*x67*x67*x67) =E= 0; e66.. 0.00437*x174/((0.01 + 0.3*x174*x174)*x67*x67) - 0.00437*cos(x68)/(( 0.0162079 + 0.486237*x174*x174)*x67*x67*x67*x67) - 0.1*x68 + 0.1*x65 =E= 0; e67.. 0.00437*cos(x71)/((0.0162079 + 0.486237*x175*x175)*x70*x70) + x66 - x69 =E= 0; e68.. (-1/x70) - cos(x71)/((0.0162079 + 0.486237*x175*x175)*x70*x70*x70) =E= 0; e69.. 0.00437*x175/((0.01 + 0.3*x175*x175)*x70*x70) - 0.00437*cos(x71)/(( 0.0162079 + 0.486237*x175*x175)*x70*x70*x70*x70) - 0.1*x71 + 0.1*x68 =E= 0; e70.. 0.00437*cos(x74)/((0.0162079 + 0.486237*x176*x176)*x73*x73) + x69 - x72 =E= 0; e71.. (-1/x73) - cos(x74)/((0.0162079 + 0.486237*x176*x176)*x73*x73*x73) =E= 0; e72.. 0.00437*x176/((0.01 + 0.3*x176*x176)*x73*x73) - 0.00437*cos(x74)/(( 0.0162079 + 0.486237*x176*x176)*x73*x73*x73*x73) - 0.1*x74 + 0.1*x71 =E= 0; e73.. 0.00437*cos(x77)/((0.0162079 + 0.486237*x177*x177)*x76*x76) + x72 - x75 =E= 0; e74.. (-1/x76) - cos(x77)/((0.0162079 + 0.486237*x177*x177)*x76*x76*x76) =E= 0; e75.. 0.00437*x177/((0.01 + 0.3*x177*x177)*x76*x76) - 0.00437*cos(x77)/(( 0.0162079 + 0.486237*x177*x177)*x76*x76*x76*x76) - 0.1*x77 + 0.1*x74 =E= 0; e76.. 0.00437*cos(x80)/((0.0162079 + 0.486237*x178*x178)*x79*x79) + x75 - x78 =E= 0; e77.. (-1/x79) - cos(x80)/((0.0162079 + 0.486237*x178*x178)*x79*x79*x79) =E= 0; e78.. 0.00437*x178/((0.01 + 0.3*x178*x178)*x79*x79) - 0.00437*cos(x80)/(( 0.0162079 + 0.486237*x178*x178)*x79*x79*x79*x79) - 0.1*x80 + 0.1*x77 =E= 0; e79.. 0.00437*cos(x83)/((0.0162079 + 0.486237*x179*x179)*x82*x82) + x78 - x81 =E= 0; e80.. (-1/x82) - cos(x83)/((0.0162079 + 0.486237*x179*x179)*x82*x82*x82) =E= 0; e81.. 0.00437*x179/((0.01 + 0.3*x179*x179)*x82*x82) - 0.00437*cos(x83)/(( 0.0162079 + 0.486237*x179*x179)*x82*x82*x82*x82) - 0.1*x83 + 0.1*x80 =E= 0; e82.. 0.00437*cos(x86)/((0.0162079 + 0.486237*x180*x180)*x85*x85) + x81 - x84 =E= 0; e83.. (-1/x85) - cos(x86)/((0.0162079 + 0.486237*x180*x180)*x85*x85*x85) =E= 0; e84.. 0.00437*x180/((0.01 + 0.3*x180*x180)*x85*x85) - 0.00437*cos(x86)/(( 0.0162079 + 0.486237*x180*x180)*x85*x85*x85*x85) - 0.1*x86 + 0.1*x83 =E= 0; e85.. 0.00437*cos(x89)/((0.0162079 + 0.486237*x181*x181)*x88*x88) + x84 - x87 =E= 0; e86.. (-1/x88) - cos(x89)/((0.0162079 + 0.486237*x181*x181)*x88*x88*x88) =E= 0; e87.. 0.00437*x181/((0.01 + 0.3*x181*x181)*x88*x88) - 0.00437*cos(x89)/(( 0.0162079 + 0.486237*x181*x181)*x88*x88*x88*x88) - 0.1*x89 + 0.1*x86 =E= 0; e88.. 0.00437*cos(x92)/((0.0162079 + 0.486237*x182*x182)*x91*x91) + x87 - x90 =E= 0; e89.. (-1/x91) - cos(x92)/((0.0162079 + 0.486237*x182*x182)*x91*x91*x91) =E= 0; e90.. 0.00437*x182/((0.01 + 0.3*x182*x182)*x91*x91) - 0.00437*cos(x92)/(( 0.0162079 + 0.486237*x182*x182)*x91*x91*x91*x91) - 0.1*x92 + 0.1*x89 =E= 0; e91.. 0.00437*cos(x95)/((0.0162079 + 0.486237*x183*x183)*x94*x94) + x90 - x93 =E= 0; e92.. (-1/x94) - cos(x95)/((0.0162079 + 0.486237*x183*x183)*x94*x94*x94) =E= 0; e93.. 0.00437*x183/((0.01 + 0.3*x183*x183)*x94*x94) - 0.00437*cos(x95)/(( 0.0162079 + 0.486237*x183*x183)*x94*x94*x94*x94) - 0.1*x95 + 0.1*x92 =E= 0; e94.. 0.00437*cos(x98)/((0.0162079 + 0.486237*x184*x184)*x97*x97) + x93 - x96 =E= 0; e95.. (-1/x97) - cos(x98)/((0.0162079 + 0.486237*x184*x184)*x97*x97*x97) =E= 0; e96.. 0.00437*x184/((0.01 + 0.3*x184*x184)*x97*x97) - 0.00437*cos(x98)/(( 0.0162079 + 0.486237*x184*x184)*x97*x97*x97*x97) - 0.1*x98 + 0.1*x95 =E= 0; e97.. 0.00437*cos(x101)/((0.0162079 + 0.486237*x185*x185)*x100*x100) + x96 - x99 =E= 0; e98.. (-1/x100) - cos(x101)/((0.0162079 + 0.486237*x185*x185)*x100*x100*x100) =E= 0; e99.. 0.00437*x185/((0.01 + 0.3*x185*x185)*x100*x100) - 0.00437*cos(x101)/(( 0.0162079 + 0.486237*x185*x185)*x100*x100*x100*x100) - 0.1*x101 + 0.1*x98 =E= 0; e100.. 0.00437*cos(x104)/((0.0162079 + 0.486237*x186*x186)*x103*x103) + x99 - x102 =E= 0; e101.. (-1/x103) - cos(x104)/((0.0162079 + 0.486237*x186*x186)*x103*x103*x103) =E= 0; e102.. 0.00437*x186/((0.01 + 0.3*x186*x186)*x103*x103) - 0.00437*cos(x104)/(( 0.0162079 + 0.486237*x186*x186)*x103*x103*x103*x103) - 0.1*x104 + 0.1*x101 =E= 0; e103.. 0.00437*cos(x107)/((0.0162079 + 0.486237*x187*x187)*x106*x106) + x102 - x105 =E= 0; e104.. (-1/x106) - cos(x107)/((0.0162079 + 0.486237*x187*x187)*x106*x106*x106) =E= 0; e105.. 0.00437*x187/((0.01 + 0.3*x187*x187)*x106*x106) - 0.00437*cos(x107)/(( 0.0162079 + 0.486237*x187*x187)*x106*x106*x106*x106) - 0.1*x107 + 0.1*x104 =E= 0; e106.. 0.00437*cos(x110)/((0.0162079 + 0.486237*x188*x188)*x109*x109) + x105 - x108 =E= 0; e107.. (-1/x109) - cos(x110)/((0.0162079 + 0.486237*x188*x188)*x109*x109*x109) =E= 0; e108.. 0.00437*x188/((0.01 + 0.3*x188*x188)*x109*x109) - 0.00437*cos(x110)/(( 0.0162079 + 0.486237*x188*x188)*x109*x109*x109*x109) - 0.1*x110 + 0.1*x107 =E= 0; e109.. 0.00437*cos(x113)/((0.0162079 + 0.486237*x189*x189)*x112*x112) + x108 - x111 =E= 0; e110.. (-1/x112) - cos(x113)/((0.0162079 + 0.486237*x189*x189)*x112*x112*x112) =E= 0; e111.. 0.00437*x189/((0.01 + 0.3*x189*x189)*x112*x112) - 0.00437*cos(x113)/(( 0.0162079 + 0.486237*x189*x189)*x112*x112*x112*x112) - 0.1*x113 + 0.1*x110 =E= 0; e112.. 0.00437*cos(x116)/((0.0162079 + 0.486237*x190*x190)*x115*x115) + x111 - x114 =E= 0; e113.. (-1/x115) - cos(x116)/((0.0162079 + 0.486237*x190*x190)*x115*x115*x115) =E= 0; e114.. 0.00437*x190/((0.01 + 0.3*x190*x190)*x115*x115) - 0.00437*cos(x116)/(( 0.0162079 + 0.486237*x190*x190)*x115*x115*x115*x115) - 0.1*x116 + 0.1*x113 =E= 0; e115.. 0.00437*cos(x119)/((0.0162079 + 0.486237*x191*x191)*x118*x118) + x114 - x117 =E= 0; e116.. (-1/x118) - cos(x119)/((0.0162079 + 0.486237*x191*x191)*x118*x118*x118) =E= 0; e117.. 0.00437*x191/((0.01 + 0.3*x191*x191)*x118*x118) - 0.00437*cos(x119)/(( 0.0162079 + 0.486237*x191*x191)*x118*x118*x118*x118) - 0.1*x119 + 0.1*x116 =E= 0; e118.. 0.00437*cos(x122)/((0.0162079 + 0.486237*x192*x192)*x121*x121) + x117 - x120 =E= 0; e119.. (-1/x121) - cos(x122)/((0.0162079 + 0.486237*x192*x192)*x121*x121*x121) =E= 0; e120.. 0.00437*x192/((0.01 + 0.3*x192*x192)*x121*x121) - 0.00437*cos(x122)/(( 0.0162079 + 0.486237*x192*x192)*x121*x121*x121*x121) - 0.1*x122 + 0.1*x119 =E= 0; e121.. 0.00437*cos(x125)/((0.0162079 + 0.486237*x193*x193)*x124*x124) + x120 - x123 =E= 0; e122.. (-1/x124) - cos(x125)/((0.0162079 + 0.486237*x193*x193)*x124*x124*x124) =E= 0; e123.. 0.00437*x193/((0.01 + 0.3*x193*x193)*x124*x124) - 0.00437*cos(x125)/(( 0.0162079 + 0.486237*x193*x193)*x124*x124*x124*x124) - 0.1*x125 + 0.1*x122 =E= 0; e124.. 0.00437*cos(x128)/((0.0162079 + 0.486237*x194*x194)*x127*x127) + x123 - x126 =E= 0; e125.. (-1/x127) - cos(x128)/((0.0162079 + 0.486237*x194*x194)*x127*x127*x127) =E= 0; e126.. 0.00437*x194/((0.01 + 0.3*x194*x194)*x127*x127) - 0.00437*cos(x128)/(( 0.0162079 + 0.486237*x194*x194)*x127*x127*x127*x127) - 0.1*x128 + 0.1*x125 =E= 0; e127.. 0.00437*cos(x131)/((0.0162079 + 0.486237*x195*x195)*x130*x130) + x126 - x129 =E= 0; e128.. (-1/x130) - cos(x131)/((0.0162079 + 0.486237*x195*x195)*x130*x130*x130) =E= 0; e129.. 0.00437*x195/((0.01 + 0.3*x195*x195)*x130*x130) - 0.00437*cos(x131)/(( 0.0162079 + 0.486237*x195*x195)*x130*x130*x130*x130) - 0.1*x131 + 0.1*x128 =E= 0; e130.. 0.00437*cos(x134)/((0.0162079 + 0.486237*x196*x196)*x133*x133) + x129 - x132 =E= 0; e131.. (-1/x133) - cos(x134)/((0.0162079 + 0.486237*x196*x196)*x133*x133*x133) =E= 0; e132.. 0.00437*x196/((0.01 + 0.3*x196*x196)*x133*x133) - 0.00437*cos(x134)/(( 0.0162079 + 0.486237*x196*x196)*x133*x133*x133*x133) - 0.1*x134 + 0.1*x131 =E= 0; e133.. 0.00437*cos(x137)/((0.0162079 + 0.486237*x197*x197)*x136*x136) + x132 - x135 =E= 0; e134.. (-1/x136) - cos(x137)/((0.0162079 + 0.486237*x197*x197)*x136*x136*x136) =E= 0; e135.. 0.00437*x197/((0.01 + 0.3*x197*x197)*x136*x136) - 0.00437*cos(x137)/(( 0.0162079 + 0.486237*x197*x197)*x136*x136*x136*x136) - 0.1*x137 + 0.1*x134 =E= 0; e136.. 0.00437*cos(x140)/((0.0162079 + 0.486237*x198*x198)*x139*x139) + x135 - x138 =E= 0; e137.. (-1/x139) - cos(x140)/((0.0162079 + 0.486237*x198*x198)*x139*x139*x139) =E= 0; e138.. 0.00437*x198/((0.01 + 0.3*x198*x198)*x139*x139) - 0.00437*cos(x140)/(( 0.0162079 + 0.486237*x198*x198)*x139*x139*x139*x139) - 0.1*x140 + 0.1*x137 =E= 0; e139.. 0.00437*cos(x143)/((0.0162079 + 0.486237*x199*x199)*x142*x142) + x138 - x141 =E= 0; e140.. (-1/x142) - cos(x143)/((0.0162079 + 0.486237*x199*x199)*x142*x142*x142) =E= 0; e141.. 0.00437*x199/((0.01 + 0.3*x199*x199)*x142*x142) - 0.00437*cos(x143)/(( 0.0162079 + 0.486237*x199*x199)*x142*x142*x142*x142) - 0.1*x143 + 0.1*x140 =E= 0; e142.. 0.00437*cos(x146)/((0.0162079 + 0.486237*x200*x200)*x145*x145) + x141 - x144 =E= 0; e143.. (-1/x145) - cos(x146)/((0.0162079 + 0.486237*x200*x200)*x145*x145*x145) =E= 0; e144.. 0.00437*x200/((0.01 + 0.3*x200*x200)*x145*x145) - 0.00437*cos(x146)/(( 0.0162079 + 0.486237*x200*x200)*x145*x145*x145*x145) - 0.1*x146 + 0.1*x143 =E= 0; e145.. 0.00437*cos(x149)/((0.0162079 + 0.486237*x201*x201)*x148*x148) + x144 - x147 =E= 0; e146.. (-1/x148) - cos(x149)/((0.0162079 + 0.486237*x201*x201)*x148*x148*x148) =E= 0; e147.. 0.00437*x201/((0.01 + 0.3*x201*x201)*x148*x148) - 0.00437*cos(x149)/(( 0.0162079 + 0.486237*x201*x201)*x148*x148*x148*x148) - 0.1*x149 + 0.1*x146 =E= 0; e148.. 0.00437*cos(x152)/((0.0162079 + 0.486237*x202*x202)*x151*x151) + x147 - x150 =E= 0; e149.. (-1/x151) - cos(x152)/((0.0162079 + 0.486237*x202*x202)*x151*x151*x151) =E= 0; e150.. 0.00437*x202/((0.01 + 0.3*x202*x202)*x151*x151) - 0.00437*cos(x152)/(( 0.0162079 + 0.486237*x202*x202)*x151*x151*x151*x151) - 0.1*x152 + 0.1*x149 =E= 0; e151.. - x150 + objvar =E= 0; * set non default bounds x1.fx = 0.005; x2.lo = 0.005; x2.up = 6.2881854; x5.lo = 0.005; x5.up = 6.2881854; x8.lo = 0.005; x8.up = 6.2881854; x11.lo = 0.005; x11.up = 6.2881854; x14.lo = 0.005; x14.up = 6.2881854; x17.lo = 0.005; x17.up = 6.2881854; x20.lo = 0.005; x20.up = 6.2881854; x23.lo = 0.005; x23.up = 6.2881854; x26.lo = 0.005; x26.up = 6.2881854; x29.lo = 0.005; x29.up = 6.2881854; x32.lo = 0.005; x32.up = 6.2881854; x35.lo = 0.005; x35.up = 6.2881854; x38.lo = 0.005; x38.up = 6.2881854; x41.lo = 0.005; x41.up = 6.2881854; x44.lo = 0.005; x44.up = 6.2881854; x47.lo = 0.005; x47.up = 6.2881854; x50.lo = 0.005; x50.up = 6.2881854; x53.lo = 0.005; x53.up = 6.2881854; x56.lo = 0.005; x56.up = 6.2881854; x59.lo = 0.005; x59.up = 6.2881854; x62.lo = 0.005; x62.up = 6.2881854; x65.lo = 0.005; x65.up = 6.2881854; x68.lo = 0.005; x68.up = 6.2881854; x71.lo = 0.005; x71.up = 6.2881854; x74.lo = 0.005; x74.up = 6.2881854; x77.lo = 0.005; x77.up = 6.2881854; x80.lo = 0.005; x80.up = 6.2881854; x83.lo = 0.005; x83.up = 6.2881854; x86.lo = 0.005; x86.up = 6.2881854; x89.lo = 0.005; x89.up = 6.2881854; x92.lo = 0.005; x92.up = 6.2881854; x95.lo = 0.005; x95.up = 6.2881854; x98.lo = 0.005; x98.up = 6.2881854; x101.lo = 0.005; x101.up = 6.2881854; x104.lo = 0.005; x104.up = 6.2881854; x107.lo = 0.005; x107.up = 6.2881854; x110.lo = 0.005; x110.up = 6.2881854; x113.lo = 0.005; x113.up = 6.2881854; x116.lo = 0.005; x116.up = 6.2881854; x119.lo = 0.005; x119.up = 6.2881854; x122.lo = 0.005; x122.up = 6.2881854; x125.lo = 0.005; x125.up = 6.2881854; x128.lo = 0.005; x128.up = 6.2881854; x131.lo = 0.005; x131.up = 6.2881854; x134.lo = 0.005; x134.up = 6.2881854; x137.lo = 0.005; x137.up = 6.2881854; x140.lo = 0.005; x140.up = 6.2881854; x143.lo = 0.005; x143.up = 6.2881854; x146.lo = 0.005; x146.up = 6.2881854; x149.lo = 0.005; x149.up = 6.2881854; x152.lo = 0.005; x152.up = 6.2881854; x153.lo = 0.08; x153.up = 0.417; x154.lo = 0.08; x154.up = 0.417; x155.lo = 0.08; x155.up = 0.417; x156.lo = 0.08; x156.up = 0.417; x157.lo = 0.08; x157.up = 0.417; x158.lo = 0.08; x158.up = 0.417; x159.lo = 0.08; x159.up = 0.417; x160.lo = 0.08; x160.up = 0.417; x161.lo = 0.08; x161.up = 0.417; x162.lo = 0.08; x162.up = 0.417; x163.lo = 0.08; x163.up = 0.417; x164.lo = 0.08; x164.up = 0.417; x165.lo = 0.08; x165.up = 0.417; x166.lo = 0.08; x166.up = 0.417; x167.lo = 0.08; x167.up = 0.417; x168.lo = 0.08; x168.up = 0.417; x169.lo = 0.08; x169.up = 0.417; x170.lo = 0.08; x170.up = 0.417; x171.lo = 0.08; x171.up = 0.417; x172.lo = 0.08; x172.up = 0.417; x173.lo = 0.08; x173.up = 0.417; x174.lo = 0.08; x174.up = 0.417; x175.lo = 0.08; x175.up = 0.417; x176.lo = 0.08; x176.up = 0.417; x177.lo = 0.08; x177.up = 0.417; x178.lo = 0.08; x178.up = 0.417; x179.lo = 0.08; x179.up = 0.417; x180.lo = 0.08; x180.up = 0.417; x181.lo = 0.08; x181.up = 0.417; x182.lo = 0.08; x182.up = 0.417; x183.lo = 0.08; x183.up = 0.417; x184.lo = 0.08; x184.up = 0.417; x185.lo = 0.08; x185.up = 0.417; x186.lo = 0.08; x186.up = 0.417; x187.lo = 0.08; x187.up = 0.417; x188.lo = 0.08; x188.up = 0.417; x189.lo = 0.08; x189.up = 0.417; x190.lo = 0.08; x190.up = 0.417; x191.lo = 0.08; x191.up = 0.417; x192.lo = 0.08; x192.up = 0.417; x193.lo = 0.08; x193.up = 0.417; x194.lo = 0.08; x194.up = 0.417; x195.lo = 0.08; x195.up = 0.417; x196.lo = 0.08; x196.up = 0.417; x197.lo = 0.08; x197.up = 0.417; x198.lo = 0.08; x198.up = 0.417; x199.lo = 0.08; x199.up = 0.417; x200.lo = 0.08; x200.up = 0.417; x201.lo = 0.08; x201.up = 0.417; x202.lo = 0.08; x202.up = 0.417; * set non default levels x4.l = 1.5; x7.l = 1.5; x10.l = 1.5; x13.l = 1.5; x16.l = 1.5; x19.l = 1.5; x22.l = 1.5; x25.l = 1.5; x28.l = 1.5; x31.l = 1.5; x34.l = 1.5; x37.l = 1.5; x40.l = 1.5; x43.l = 1.5; x46.l = 1.5; x49.l = 1.5; x52.l = 1.5; x55.l = 1.5; x58.l = 1.5; x61.l = 1.5; x64.l = 1.5; x67.l = 1.5; x70.l = 1.5; x73.l = 1.5; x76.l = 1.5; x79.l = 1.5; x82.l = 1.5; x85.l = 1.5; x88.l = 1.5; x91.l = 1.5; x94.l = 1.5; x97.l = 1.5; x100.l = 1.5; x103.l = 1.5; x106.l = 1.5; x109.l = 1.5; x112.l = 1.5; x115.l = 1.5; x118.l = 1.5; x121.l = 1.5; x124.l = 1.5; x127.l = 1.5; x130.l = 1.5; x133.l = 1.5; x136.l = 1.5; x139.l = 1.5; x142.l = 1.5; x145.l = 1.5; x148.l = 1.5; x151.l = 1.09905; * 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;