bc7.gms
<![CDATA[
* NLP written by GAMS Convert at 11/29/04 15:59:52
*
* Equation counts
* Total E G L N X C
* 129 129 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 194 194 0 0 0 0 0 0
* FX 2 2 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 385 321 64 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,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,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;
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;
e1.. 68719476736*x1 - 68719476736*x2 + x130 =E= 0;
e2.. 1090785345.01587*x2 - 1090785345.01587*x3 + x131 =E= 0;
e3.. 103337559.001504*x3 - 103337559.001504*x4 + x132 =E= 0;
e4.. 20409705.000297*x4 - 20409705.000297*x5 + x133 =E= 0;
e5.. 5960575.65582444*x5 - 5960575.65582444*x6 + x134 =E= 0;
e6.. 2214542.77129322*x6 - 2214542.77129322*x7 + x135 =E= 0;
e7.. 967975.388221374*x7 - 967975.388221374*x8 + x136 =E= 0;
e8.. 475583.76923769*x8 - 475583.76923769*x9 + x137 =E= 0;
e9.. 255180.996208647*x9 - 255180.996208647*x10 + x138 =E= 0;
e10.. 146661.309965234*x10 - 146661.309965234*x11 + x139 =E= 0;
e11.. 89065.5135964622*x11 - 89065.5135964622*x12 + x140 =E= 0;
e12.. 56586.1126938472*x12 - 56586.1126938472*x13 + x141 =E= 0;
e13.. 37330.80370812*x13 - 37330.80370812*x14 + x142 =E= 0;
e14.. 25425.9778127795*x14 - 25425.9778127795*x15 + x143 =E= 0;
e15.. 17797.9520119842*x15 - 17797.9520119842*x16 + x144 =E= 0;
e16.. 12757.5078070713*x16 - 12757.5078070713*x17 + x145 =E= 0;
e17.. 9336.43763227117*x17 - 9336.43763227117*x18 + x146 =E= 0;
e18.. 6959.17748377032*x18 - 6959.17748377032*x19 + x147 =E= 0;
e19.. 5272.46318788349*x19 - 5272.46318788349*x20 + x148 =E= 0;
e20.. 4053.26143670456*x20 - 4053.26143670456*x21 + x149 =E= 0;
e21.. 3157.17608736991*x21 - 3157.17608736991*x22 + x150 =E= 0;
e22.. 2488.59334977754*x22 - 2488.59334977754*x23 + x151 =E= 0;
e23.. 1982.90358031953*x23 - 1982.90358031953*x24 + x152 =E= 0;
e24.. 1595.63791105723*x24 - 1595.63791105723*x25 + x153 =E= 0;
e25.. 1295.67350800183*x25 - 1295.67350800183*x26 + x154 =E= 0;
e26.. 1060.89257493201*x26 - 1060.89257493201*x27 + x155 =E= 0;
e27.. 875.354792214832*x27 - 875.354792214832*x28 + x156 =E= 0;
e28.. 727.42258186914*x28 - 727.42258186914*x29 + x157 =E= 0;
e29.. 608.497661370368*x29 - 608.497661370368*x30 + x158 =E= 0;
e30.. 512.156637413868*x30 - 512.156637413868*x31 + x159 =E= 0;
e31.. 433.55129863514*x31 - 433.55129863514*x32 + x160 =E= 0;
e32.. 368.987123846053*x32 - 368.987123846053*x33 + x161 =E= 0;
e33.. 315.62344860329*x33 - 315.62344860329*x34 + x162 =E= 0;
e34.. 271.257758406946*x34 - 271.257758406946*x35 + x163 =E= 0;
e35.. 234.168859898618*x35 - 234.168859898618*x36 + x164 =E= 0;
e36.. 203.001726363813*x36 - 203.001726363813*x37 + x165 =E= 0;
e37.. 176.682154341506*x37 - 176.682154341506*x38 + x166 =E= 0;
e38.. 154.352958367566*x38 - 154.352958367566*x39 + x167 =E= 0;
e39.. 135.325873251345*x39 - 135.325873251345*x40 + x168 =E= 0;
e40.. 119.045013450257*x40 - 119.045013450257*x41 + x169 =E= 0;
e41.. 105.058907171953*x41 - 105.058907171953*x42 + x170 =E= 0;
e42.. 92.9989430045616*x42 - 92.9989430045616*x43 + x171 =E= 0;
e43.. 82.562648218548*x43 - 82.562648218548*x44 + x172 =E= 0;
e44.. 73.5006336392199*x44 - 73.5006336392199*x45 + x173 =E= 0;
e45.. 65.6063398523889*x45 - 65.6063398523889*x46 + x174 =E= 0;
e46.. 58.7079375310427*x46 - 58.7079375310427*x47 + x175 =E= 0;
e47.. 52.6618944128288*x47 - 52.6618944128288*x48 + x176 =E= 0;
e48.. 47.347839355123*x48 - 47.347839355123*x49 + x177 =E= 0;
e49.. 42.6644415161561*x49 - 42.6644415161561*x50 + x178 =E= 0;
e50.. 38.5260882662983*x50 - 38.5260882662983*x51 + x179 =E= 0;
e51.. 34.8601947930382*x51 - 34.8601947930382*x52 + x180 =E= 0;
e52.. 31.6050157547443*x52 - 31.6050157547443*x53 + x181 =E= 0;
e53.. 28.7078578293424*x53 - 28.7078578293424*x54 + x182 =E= 0;
e54.. 26.1236138348849*x54 - 26.1236138348849*x55 + x183 =E= 0;
e55.. 23.8135559164912*x55 - 23.8135559164912*x56 + x184 =E= 0;
e56.. 21.7443383164886*x56 - 21.7443383164886*x57 + x185 =E= 0;
e57.. 19.8871703782655*x57 - 19.8871703782655*x58 + x186 =E= 0;
e58.. 18.2171283577676*x58 - 18.2171283577676*x59 + x187 =E= 0;
e59.. 16.7125808402805*x59 - 16.7125808402805*x60 + x188 =E= 0;
e60.. 15.3547074703863*x60 - 15.3547074703863*x61 + x189 =E= 0;
e61.. 14.1270945934911*x61 - 14.1270945934911*x62 + x190 =E= 0;
e62.. 13.0153945025384*x62 - 13.0153945025384*x63 + x191 =E= 0;
e63.. 12.0070374558166*x63 - 12.0070374558166*x64 + x192 =E= 0;
e64.. 11.0909876140619*x64 - 11.0909876140619*x65 + x193 =E= 0;
e65.. - sqrt(137438953472 + 137438953472*sqr(x130)) + x66 =E= 0;
e66.. - sqrt(2114445438.03077 + 2114445438.03077*sqr(x131)) + x67 =E= 0;
e67.. - sqrt(173315199.838588 + 173315199.838588*sqr(x132)) + x68 =E= 0;
e68.. - sqrt(28484757.1962694 + 28484757.1962694*sqr(x133)) + x69 =E= 0;
e69.. - sqrt(6969167.56107702 + 6969167.56107702*sqr(x134)) + x70 =E= 0;
e70.. - sqrt(2206755.72762159 + 2206755.72762159*sqr(x135)) + x71 =E= 0;
e71.. - sqrt(836486.737908159 + 836486.737908159*sqr(x136)) + x72 =E= 0;
e72.. - sqrt(361878.585102938 + 361878.585102938*sqr(x137)) + x73 =E= 0;
e73.. - sqrt(173187.438613381 + 173187.438613381*sqr(x138)) + x74 =E= 0;
e74.. - sqrt(89744.8569497617 + 89744.8569497617*sqr(x139)) + x75 =E= 0;
e75.. - sqrt(49589.0054276272 + 49589.0054276272*sqr(x140)) + x76 =E= 0;
e76.. - sqrt(28888.6292135965 + 28888.6292135965*sqr(x141)) + x77 =E= 0;
e77.. - sqrt(17591.526291814 + 17591.526291814*sqr(x142)) + x78 =E= 0;
e78.. - sqrt(11122.9456179801 + 11122.9456179801*sqr(x143)) + x79 =E= 0;
e79.. - sqrt(7264.15348537468 + 7264.15348537468*sqr(x144)) + x80 =E= 0;
e80.. - sqrt(4879.28604368365 + 4879.28604368365*sqr(x145)) + x81 =E= 0;
e81.. - sqrt(3359.15130611098 + 3359.15130611098*sqr(x146)) + x82 =E= 0;
e82.. - sqrt(2363.53297890502 + 2363.53297890502*sqr(x147)) + x83 =E= 0;
e83.. - sqrt(1695.56090007286 + 1695.56090007286*sqr(x148)) + x84 =E= 0;
e84.. - sqrt(1237.67718563105 + 1237.67718563105*sqr(x149)) + x85 =E= 0;
e85.. - sqrt(917.690546795961 + 917.690546795961*sqr(x150)) + x86 =E= 0;
e86.. - sqrt(690.141585663083 + 690.141585663083*sqr(x151)) + x87 =E= 0;
e87.. - sqrt(525.748471026768 + 525.748471026768*sqr(x152)) + x88 =E= 0;
e88.. - sqrt(405.258634901665 + 405.258634901665*sqr(x153)) + x89 =E= 0;
e89.. - sqrt(315.774782572852 + 315.774782572852*sqr(x154)) + x90 =E= 0;
e90.. - sqrt(248.508024178894 + 248.508024178894*sqr(x155)) + x91 =E= 0;
e91.. - sqrt(197.374401392121 + 197.374401392121*sqr(x156)) + x92 =E= 0;
e92.. - sqrt(158.101054972177 + 158.101054972177*sqr(x157)) + x93 =E= 0;
e93.. - sqrt(127.646711512544 + 127.646711512544*sqr(x158)) + x94 =E= 0;
e94.. - sqrt(103.819710146955 + 103.819710146955*sqr(x159)) + x95 =E= 0;
e95.. - sqrt(85.0223572562344 + 85.0223572562344*sqr(x160)) + x96 =E= 0;
e96.. - sqrt(70.0773835410269 + 70.0773835410269*sqr(x161)) + x97 =E= 0;
e97.. - sqrt(58.1085677383714 + 58.1085677383714*sqr(x162)) + x98 =E= 0;
e98.. - sqrt(48.4576002639464 + 48.4576002639464*sqr(x163)) + x99 =E= 0;
e99.. - sqrt(40.6255122732766 + 40.6255122732766*sqr(x164)) + x100 =E= 0;
e100.. - sqrt(34.2309618233724 + 34.2309618233724*sqr(x165)) + x101 =E= 0;
e101.. - sqrt(28.9802214106408 + 28.9802214106408*sqr(x166)) + x102 =E= 0;
e102.. - sqrt(24.6453763789551 + 24.6453763789551*sqr(x167)) + x103 =E= 0;
e103.. - sqrt(21.0483439341576 + 21.0483439341576*sqr(x168)) + x104 =E= 0;
e104.. - sqrt(18.0490582199119 + 18.0490582199119*sqr(x169)) + x105 =E= 0;
e105.. - sqrt(15.5366644714627 + 15.5366644714627*sqr(x170)) + x106 =E= 0;
e106.. - sqrt(13.4229053773037 + 13.4229053773037*sqr(x171)) + x107 =E= 0;
e107.. - sqrt(11.6371176307722 + 11.6371176307722*sqr(x172)) + x108 =E= 0;
e108.. - sqrt(10.1224203841077 + 10.1224203841077*sqr(x173)) + x109 =E= 0;
e109.. - sqrt(8.8327925085359 + 8.8327925085359*sqr(x174)) + x110 =E= 0;
e110.. - sqrt(7.730817309797 + 7.730817309797*sqr(x175)) + x111 =E= 0;
e111.. - sqrt(6.78593183963183 + 6.78593183963183*sqr(x176)) + x112 =E= 0;
e112.. - sqrt(5.9730601252537 + 5.9730601252537*sqr(x177)) + x113 =E= 0;
e113.. - sqrt(5.27154028712339 + 5.27154028712339*sqr(x178)) + x114 =E= 0;
e114.. - sqrt(4.66427794362012 + 4.66427794362012*sqr(x179)) + x115 =E= 0;
e115.. - sqrt(4.13707482669781 + 4.13707482669781*sqr(x180)) + x116 =E= 0;
e116.. - sqrt(3.67809378878012 + 3.67809378878012*sqr(x181)) + x117 =E= 0;
e117.. - sqrt(3.27743052810095 + 3.27743052810095*sqr(x182)) + x118 =E= 0;
e118.. - sqrt(2.92676922734499 + 2.92676922734499*sqr(x183)) + x119 =E= 0;
e119.. - sqrt(2.61910448657899 + 2.61910448657899*sqr(x184)) + x120 =E= 0;
e120.. - sqrt(2.34851586954992 + 2.34851586954992*sqr(x185)) + x121 =E= 0;
e121.. - sqrt(2.10998438876693 + 2.10998438876693*sqr(x186)) + x122 =E= 0;
e122.. - sqrt(1.89924256161383 + 1.89924256161383*sqr(x187)) + x123 =E= 0;
e123.. - sqrt(1.71265144856549 + 1.71265144856549*sqr(x188)) + x124 =E= 0;
e124.. - sqrt(1.54709946278868 + 1.54709946278868*sqr(x189)) + x125 =E= 0;
e125.. - sqrt(1.39991881312126 + 1.39991881312126*sqr(x190)) + x126 =E= 0;
e126.. - sqrt(1.26881628105478 + 1.26881628105478*sqr(x191)) + x127 =E= 0;
e127.. - sqrt(1.15181569077584 + 1.15181569077584*sqr(x192)) + x128 =E= 0;
e128.. - sqrt(1.04720995040502 + 1.04720995040502*sqr(x193)) + x129 =E= 0;
e129.. - 1.45519152283669E-11*x66 - 9.16770659387112E-10*x67
- 9.67702362686396E-9*x68 - 4.89962985739112E-8*x69
- 1.67769030667841E-7*x70 - 4.51560481451452E-7*x71
- 1.03308411780745E-6*x72 - 2.10267899092287E-6*x73
- 3.91878711525351E-6*x74 - 6.81843084748834E-6*x75
- 1.1227690265514E-5*x76 - 1.7672180547379E-5*x77
- 2.67875293502584E-5*x78 - 3.93298541894183E-5*x79
- 5.61862398171797E-5*x80 - 7.83852156018838E-5*x81
- 0.000107107232906856*x82 - 0.000143695142469369*x83
- 0.00018966467177961*x84 - 0.000246714902459644*x85
- 0.000316738747642376*x86 - 0.000401833429350518*x87
- 0.000504310955875553*x88 - 0.0006267085991567*x89
- 0.000771799372159876*x90 - 0.000942602506256662*x91
- 0.00114239392860327*x92 - 0.0013747167395195*x93
- 0.00164339168986771*x94 - 0.00195252765843179*x95
- 0.0023065321292961*x96 - 0.00271012166922446*x97
- 0.00316833240503911*x98 - 0.00368653050099965*x99
- 0.00427042263618205*x100 - 0.00492606648185756*x101
- 0.00565988117887173*x102 - 0.00647865781502333*x103
- 0.00738956990244333*x104 - 0.00840018385497387*x105
- 0.00951846946554724*x106 - 0.0107528103835648*x107
- 0.012112014592276*x108 - 0.0136053248861572*x109
- 0.0152424293482909*x110 - 0.0170334718277445*x111
- 0.0189890624169493*x112 - 0.0211202879290795*x113
- 0.0234387223754311*x114 - 0.025956437442801*x115
- 0.0286860129708657*x116 - 0.0316405474295607*x117
- 0.034833668396459*x118 - 0.0382795430341503*x119
- 0.04199288856762*x120 - 0.045988982761628*x121
- 0.0502836743980879*x122 - 0.0548933937534457*x123
- 0.059835163076059*x124 - 0.0651266070635756*x125
- 0.0707859633403132*x126 - 0.0768320929346373*x127
- 0.0832844907563413*x128 - 0.0901632960740244*x129 + objvar =E= 0;
* set non default bounds
x1.fx = 0;
x65.fx = 1;
x66.up = 100000000;
x67.up = 100000000;
x68.up = 100000000;
x69.up = 100000000;
x70.up = 100000000;
x71.up = 100000000;
x72.up = 100000000;
x73.up = 100000000;
x74.up = 100000000;
x75.up = 100000000;
x76.up = 100000000;
x77.up = 100000000;
x78.up = 100000000;
x79.up = 100000000;
x80.up = 100000000;
x81.up = 100000000;
x82.up = 100000000;
x83.up = 100000000;
x84.up = 100000000;
x85.up = 100000000;
x86.up = 100000000;
x87.up = 100000000;
x88.up = 100000000;
x89.up = 100000000;
x90.up = 100000000;
x91.up = 100000000;
x92.up = 100000000;
x93.up = 100000000;
x94.up = 100000000;
x95.up = 100000000;
x96.up = 100000000;
x97.up = 100000000;
x98.up = 100000000;
x99.up = 100000000;
x100.up = 100000000;
x101.up = 100000000;
x102.up = 100000000;
x103.up = 100000000;
x104.up = 100000000;
x105.up = 100000000;
x106.up = 100000000;
x107.up = 100000000;
x108.up = 100000000;
x109.up = 100000000;
x110.up = 100000000;
x111.up = 100000000;
x112.up = 100000000;
x113.up = 100000000;
x114.up = 100000000;
x115.up = 100000000;
x116.up = 100000000;
x117.up = 100000000;
x118.up = 100000000;
x119.up = 100000000;
x120.up = 100000000;
x121.up = 100000000;
x122.up = 100000000;
x123.up = 100000000;
x124.up = 100000000;
x125.up = 100000000;
x126.up = 100000000;
x127.up = 100000000;
x128.up = 100000000;
x129.up = 100000000;
* set non default levels
x2.l = 1.45519152283669E-11;
x3.l = 9.31322574615479E-10;
x4.l = 1.06083462014794E-8;
x5.l = 5.96046447753906E-8;
x6.l = 2.27373675443232E-7;
x7.l = 6.78934156894684E-7;
x8.l = 1.71201827470213E-6;
x9.l = 3.814697265625E-6;
x10.l = 7.73348438087851E-6;
x11.l = 1.45519152283669E-5;
x12.l = 2.57796054938808E-5;
x13.l = 4.34517860412598E-5;
x14.l = 7.02393153915182E-5;
x15.l = 0.000109569169580936;
x16.l = 0.000165755409398116;
x17.l = 0.000244140625;
x18.l = 0.000351247857906856;
x19.l = 0.000494943000376225;
x20.l = 0.000684607672155835;
x21.l = 0.000931322574615479;
x22.l = 0.00124806132225785;
x23.l = 0.00164989475160837;
x24.l = 0.00215420570748392;
x25.l = 0.00278091430664062;
x26.l = 0.0035527136788005;
x27.l = 0.00449531618505716;
x28.l = 0.00563771011366043;
x29.l = 0.00701242685317993;
x30.l = 0.00865581854304764;
x31.l = 0.0106083462014794;
x32.l = 0.0129148783307755;
x33.l = 0.015625;
x34.l = 0.0187933324050391;
x35.l = 0.0224798629060388;
x36.l = 0.0267502855422208;
x37.l = 0.0316763520240784;
x38.l = 0.0373362332029501;
x39.l = 0.0438148910179734;
x40.l = 0.0512044609204168;
x41.l = 0.0596046447753906;
x42.l = 0.0691231142409379;
x43.l = 0.0798759246245027;
x44.l = 0.0919879392167786;
x45.l = 0.105593264102936;
x46.l = 0.120835693451227;
x47.l = 0.137869165278971;
x48.l = 0.156858227695921;
x49.l = 0.177978515625;
x50.l = 0.201417238000431;
x51.l = 0.227373675443232;
x52.l = 0.256059688414098;
x53.l = 0.287700235843658;
x54.l = 0.322533904240117;
x55.l = 0.360813447274268;
x56.l = 0.402806335841888;
x57.l = 0.448795318603516;
x58.l = 0.499078993001604;
x59.l = 0.553972386755049;
x60.l = 0.613807549831108;
x61.l = 0.678934156894684;
x62.l = 0.749720120234997;
x63.l = 0.826552213169634;
x64.l = 0.909836703925976;
* 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>