catenary.gms
<![CDATA[
* NLP written by GAMS Convert at 11/29/04 16:00:03
*
* Equation counts
* Total E G L N X C
* 102 1 0 101 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 4 4 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 504 102 402 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;
Positive Variables x1,x102,x202;
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;
e1.. sqr(x1) + sqr(x102) =L= 0.0004;
e2.. sqr(x2 - x1) + sqr(x103 - x102) =L= 0.0004;
e3.. sqr(x3 - x2) + sqr(x104 - x103) =L= 0.0004;
e4.. sqr(x4 - x3) + sqr(x105 - x104) =L= 0.0004;
e5.. sqr(x5 - x4) + sqr(x106 - x105) =L= 0.0004;
e6.. sqr(x6 - x5) + sqr(x107 - x106) =L= 0.0004;
e7.. sqr(x7 - x6) + sqr(x108 - x107) =L= 0.0004;
e8.. sqr(x8 - x7) + sqr(x109 - x108) =L= 0.0004;
e9.. sqr(x9 - x8) + sqr(x110 - x109) =L= 0.0004;
e10.. sqr(x10 - x9) + sqr(x111 - x110) =L= 0.0004;
e11.. sqr(x11 - x10) + sqr(x112 - x111) =L= 0.0004;
e12.. sqr(x12 - x11) + sqr(x113 - x112) =L= 0.0004;
e13.. sqr(x13 - x12) + sqr(x114 - x113) =L= 0.0004;
e14.. sqr(x14 - x13) + sqr(x115 - x114) =L= 0.0004;
e15.. sqr(x15 - x14) + sqr(x116 - x115) =L= 0.0004;
e16.. sqr(x16 - x15) + sqr(x117 - x116) =L= 0.0004;
e17.. sqr(x17 - x16) + sqr(x118 - x117) =L= 0.0004;
e18.. sqr(x18 - x17) + sqr(x119 - x118) =L= 0.0004;
e19.. sqr(x19 - x18) + sqr(x120 - x119) =L= 0.0004;
e20.. sqr(x20 - x19) + sqr(x121 - x120) =L= 0.0004;
e21.. sqr(x21 - x20) + sqr(x122 - x121) =L= 0.0004;
e22.. sqr(x22 - x21) + sqr(x123 - x122) =L= 0.0004;
e23.. sqr(x23 - x22) + sqr(x124 - x123) =L= 0.0004;
e24.. sqr(x24 - x23) + sqr(x125 - x124) =L= 0.0004;
e25.. sqr(x25 - x24) + sqr(x126 - x125) =L= 0.0004;
e26.. sqr(x26 - x25) + sqr(x127 - x126) =L= 0.0004;
e27.. sqr(x27 - x26) + sqr(x128 - x127) =L= 0.0004;
e28.. sqr(x28 - x27) + sqr(x129 - x128) =L= 0.0004;
e29.. sqr(x29 - x28) + sqr(x130 - x129) =L= 0.0004;
e30.. sqr(x30 - x29) + sqr(x131 - x130) =L= 0.0004;
e31.. sqr(x31 - x30) + sqr(x132 - x131) =L= 0.0004;
e32.. sqr(x32 - x31) + sqr(x133 - x132) =L= 0.0004;
e33.. sqr(x33 - x32) + sqr(x134 - x133) =L= 0.0004;
e34.. sqr(x34 - x33) + sqr(x135 - x134) =L= 0.0004;
e35.. sqr(x35 - x34) + sqr(x136 - x135) =L= 0.0004;
e36.. sqr(x36 - x35) + sqr(x137 - x136) =L= 0.0004;
e37.. sqr(x37 - x36) + sqr(x138 - x137) =L= 0.0004;
e38.. sqr(x38 - x37) + sqr(x139 - x138) =L= 0.0004;
e39.. sqr(x39 - x38) + sqr(x140 - x139) =L= 0.0004;
e40.. sqr(x40 - x39) + sqr(x141 - x140) =L= 0.0004;
e41.. sqr(x41 - x40) + sqr(x142 - x141) =L= 0.0004;
e42.. sqr(x42 - x41) + sqr(x143 - x142) =L= 0.0004;
e43.. sqr(x43 - x42) + sqr(x144 - x143) =L= 0.0004;
e44.. sqr(x44 - x43) + sqr(x145 - x144) =L= 0.0004;
e45.. sqr(x45 - x44) + sqr(x146 - x145) =L= 0.0004;
e46.. sqr(x46 - x45) + sqr(x147 - x146) =L= 0.0004;
e47.. sqr(x47 - x46) + sqr(x148 - x147) =L= 0.0004;
e48.. sqr(x48 - x47) + sqr(x149 - x148) =L= 0.0004;
e49.. sqr(x49 - x48) + sqr(x150 - x149) =L= 0.0004;
e50.. sqr(x50 - x49) + sqr(x151 - x150) =L= 0.0004;
e51.. sqr(x51 - x50) + sqr(x152 - x151) =L= 0.0004;
e52.. sqr(x52 - x51) + sqr(x153 - x152) =L= 0.0004;
e53.. sqr(x53 - x52) + sqr(x154 - x153) =L= 0.0004;
e54.. sqr(x54 - x53) + sqr(x155 - x154) =L= 0.0004;
e55.. sqr(x55 - x54) + sqr(x156 - x155) =L= 0.0004;
e56.. sqr(x56 - x55) + sqr(x157 - x156) =L= 0.0004;
e57.. sqr(x57 - x56) + sqr(x158 - x157) =L= 0.0004;
e58.. sqr(x58 - x57) + sqr(x159 - x158) =L= 0.0004;
e59.. sqr(x59 - x58) + sqr(x160 - x159) =L= 0.0004;
e60.. sqr(x60 - x59) + sqr(x161 - x160) =L= 0.0004;
e61.. sqr(x61 - x60) + sqr(x162 - x161) =L= 0.0004;
e62.. sqr(x62 - x61) + sqr(x163 - x162) =L= 0.0004;
e63.. sqr(x63 - x62) + sqr(x164 - x163) =L= 0.0004;
e64.. sqr(x64 - x63) + sqr(x165 - x164) =L= 0.0004;
e65.. sqr(x65 - x64) + sqr(x166 - x165) =L= 0.0004;
e66.. sqr(x66 - x65) + sqr(x167 - x166) =L= 0.0004;
e67.. sqr(x67 - x66) + sqr(x168 - x167) =L= 0.0004;
e68.. sqr(x68 - x67) + sqr(x169 - x168) =L= 0.0004;
e69.. sqr(x69 - x68) + sqr(x170 - x169) =L= 0.0004;
e70.. sqr(x70 - x69) + sqr(x171 - x170) =L= 0.0004;
e71.. sqr(x71 - x70) + sqr(x172 - x171) =L= 0.0004;
e72.. sqr(x72 - x71) + sqr(x173 - x172) =L= 0.0004;
e73.. sqr(x73 - x72) + sqr(x174 - x173) =L= 0.0004;
e74.. sqr(x74 - x73) + sqr(x175 - x174) =L= 0.0004;
e75.. sqr(x75 - x74) + sqr(x176 - x175) =L= 0.0004;
e76.. sqr(x76 - x75) + sqr(x177 - x176) =L= 0.0004;
e77.. sqr(x77 - x76) + sqr(x178 - x177) =L= 0.0004;
e78.. sqr(x78 - x77) + sqr(x179 - x178) =L= 0.0004;
e79.. sqr(x79 - x78) + sqr(x180 - x179) =L= 0.0004;
e80.. sqr(x80 - x79) + sqr(x181 - x180) =L= 0.0004;
e81.. sqr(x81 - x80) + sqr(x182 - x181) =L= 0.0004;
e82.. sqr(x82 - x81) + sqr(x183 - x182) =L= 0.0004;
e83.. sqr(x83 - x82) + sqr(x184 - x183) =L= 0.0004;
e84.. sqr(x84 - x83) + sqr(x185 - x184) =L= 0.0004;
e85.. sqr(x85 - x84) + sqr(x186 - x185) =L= 0.0004;
e86.. sqr(x86 - x85) + sqr(x187 - x186) =L= 0.0004;
e87.. sqr(x87 - x86) + sqr(x188 - x187) =L= 0.0004;
e88.. sqr(x88 - x87) + sqr(x189 - x188) =L= 0.0004;
e89.. sqr(x89 - x88) + sqr(x190 - x189) =L= 0.0004;
e90.. sqr(x90 - x89) + sqr(x191 - x190) =L= 0.0004;
e91.. sqr(x91 - x90) + sqr(x192 - x191) =L= 0.0004;
e92.. sqr(x92 - x91) + sqr(x193 - x192) =L= 0.0004;
e93.. sqr(x93 - x92) + sqr(x194 - x193) =L= 0.0004;
e94.. sqr(x94 - x93) + sqr(x195 - x194) =L= 0.0004;
e95.. sqr(x95 - x94) + sqr(x196 - x195) =L= 0.0004;
e96.. sqr(x96 - x95) + sqr(x197 - x196) =L= 0.0004;
e97.. sqr(x97 - x96) + sqr(x198 - x197) =L= 0.0004;
e98.. sqr(x98 - x97) + sqr(x199 - x198) =L= 0.0004;
e99.. sqr(x99 - x98) + sqr(x200 - x199) =L= 0.0004;
e100.. sqr(x100 - x99) + sqr(x201 - x200) =L= 0.0004;
e101.. sqr(x101 - x100) + sqr(x202 - x201) =L= 0.0004;
e102.. - 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
- 0.5*x202 + objvar =E= 0;
* set non default bounds
x1.fx = 0;
x101.fx = 1;
x102.fx = 0;
x202.fx = 0;
* set non default levels
x2.l = 0.01;
x3.l = 0.02;
x4.l = 0.03;
x5.l = 0.04;
x6.l = 0.05;
x7.l = 0.06;
x8.l = 0.07;
x9.l = 0.08;
x10.l = 0.09;
x11.l = 0.1;
x12.l = 0.11;
x13.l = 0.12;
x14.l = 0.13;
x15.l = 0.14;
x16.l = 0.15;
x17.l = 0.16;
x18.l = 0.17;
x19.l = 0.18;
x20.l = 0.19;
x21.l = 0.2;
x22.l = 0.21;
x23.l = 0.22;
x24.l = 0.23;
x25.l = 0.24;
x26.l = 0.25;
x27.l = 0.26;
x28.l = 0.27;
x29.l = 0.28;
x30.l = 0.29;
x31.l = 0.3;
x32.l = 0.31;
x33.l = 0.32;
x34.l = 0.33;
x35.l = 0.34;
x36.l = 0.35;
x37.l = 0.36;
x38.l = 0.37;
x39.l = 0.38;
x40.l = 0.39;
x41.l = 0.4;
x42.l = 0.41;
x43.l = 0.42;
x44.l = 0.43;
x45.l = 0.44;
x46.l = 0.45;
x47.l = 0.46;
x48.l = 0.47;
x49.l = 0.48;
x50.l = 0.49;
x51.l = 0.5;
x52.l = 0.51;
x53.l = 0.52;
x54.l = 0.53;
x55.l = 0.54;
x56.l = 0.55;
x57.l = 0.56;
x58.l = 0.57;
x59.l = 0.58;
x60.l = 0.59;
x61.l = 0.6;
x62.l = 0.61;
x63.l = 0.62;
x64.l = 0.63;
x65.l = 0.64;
x66.l = 0.65;
x67.l = 0.66;
x68.l = 0.67;
x69.l = 0.68;
x70.l = 0.69;
x71.l = 0.7;
x72.l = 0.71;
x73.l = 0.72;
x74.l = 0.73;
x75.l = 0.74;
x76.l = 0.75;
x77.l = 0.76;
x78.l = 0.77;
x79.l = 0.78;
x80.l = 0.79;
x81.l = 0.8;
x82.l = 0.81;
x83.l = 0.82;
x84.l = 0.83;
x85.l = 0.84;
x86.l = 0.85;
x87.l = 0.86;
x88.l = 0.87;
x89.l = 0.88;
x90.l = 0.89;
x91.l = 0.9;
x92.l = 0.91;
x93.l = 0.92;
x94.l = 0.93;
x95.l = 0.94;
x96.l = 0.95;
x97.l = 0.96;
x98.l = 0.97;
x99.l = 0.98;
x100.l = 0.99;
* set non default marginals
Model m / all /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include '%gams.u1%'
Solve m using NLP minimizing objvar;
]]></plaintext></body>