hvycrash.gms
<![CDATA[
* 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;
]]></plaintext></body>