linspanh.gms
* NLP written by GAMS Convert at 10/06/06 11:33:51
*
* Equation counts
* Total E G L N X C
* 34 34 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 98 98 0 0 0 0 0 0
* FX 16 16 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 196 196 0 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,objvar;
Positive Variables x5,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24
,x25,x26,x27,x28,x29,x34,x40,x41,x42,x43,x44,x45,x46,x47,x48,x49,x50
,x51,x52,x53,x54,x55,x56,x57,x58,x63,x69,x70,x71,x72,x73,x74,x75,x76
,x77,x78,x79,x80,x81,x82,x83,x84,x85,x86,x87,x92;
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;
e1.. x1 - x11 - x21 - x30 =E= -51.38;
e2.. x2 + x11 - x12 + x21 - x22 - x31 =E= -13.84;
e3.. x3 + x12 - x13 + x22 - x23 - x32 =E= -2.58;
e4.. x4 - x14 - x24 - x33 =E= -21.91;
e5.. x5 + x13 + x14 - x15 - x34 =E= 0;
e6.. x6 - x16 - x25 - x35 =E= -12.97;
e7.. x7 + x16 - x17 + x25 - x26 - x36 =E= 0;
e8.. x8 + x17 - x18 + x26 - x27 - x37 =E= -2.89;
e9.. x9 + x15 + x18 - x19 + x23 + x24 + x27 - x28 - x38 =E= -20.84;
e10.. x10 + x19 - x20 + x28 - x29 - x39 =E= -17.14;
e11.. x30 - x40 - x50 - x59 =E= -32.06;
e12.. x31 + x40 - x41 + x50 - x51 - x60 =E= -0.28;
e13.. x32 + x41 - x42 + x51 - x52 - x61 =E= -4.2;
e14.. x33 - x43 - x53 - x62 =E= -48.37;
e15.. x34 + x42 + x43 - x44 - x63 =E= 0;
e16.. x35 - x45 - x54 - x64 =E= -18.13;
e17.. x36 + x45 - x46 + x54 - x55 - x65 =E= 0;
e18.. x37 + x46 - x47 + x55 - x56 - x66 =E= 1.61;
e19.. x38 + x44 + x47 - x48 + x52 + x53 + x56 - x57 - x67 =E= -26.6;
e20.. x39 + x48 - x49 + x57 - x58 - x68 =E= -18.76;
e21.. x59 - x69 - x79 - x88 =E= -18.13;
e22.. x60 + x69 - x70 + x79 - x80 - x89 =E= 0;
e23.. x61 + x70 - x71 + x80 - x81 - x90 =E= 0;
e24.. x62 - x72 - x82 - x91 =E= -18.13;
e25.. x63 + x71 + x72 - x73 - x92 =E= 0;
e26.. x64 - x74 - x83 - x93 =E= -9.1;
e27.. x65 + x74 - x75 + x83 - x84 - x94 =E= 0;
e28.. x66 + x75 - x76 + x84 - x85 - x95 =E= 5.81;
e29.. x67 + x73 + x76 - x77 + x81 + x82 + x85 - x86 - x96 =E= -9.1;
e30.. x68 + x77 - x78 + x86 - x87 - x97 =E= -6.02;
e31.. x20 + x29 + x49 + x58 + x78 + x87 =E= 608.35;
e32.. - x1 - x2 - x3 - x4 - x5 - x6 - x7 - x8 - x9 - x10 =E= -4626.34;
e33.. x88 + x89 + x90 + x91 + x92 + x93 + x94 + x95 + x96 + x97 =E= 4363;
e34.. x1 + objvar =E= 0;
* set non default bounds
x1.lo = 77; x1.up = 77.01;
x2.lo = 1124.52; x2.up = 1124.53;
x3.lo = 158; x3.up = 158.01;
x4.lo = 16; x4.up = 16.01;
x5.fx = 0;
x6.lo = 783.65; x6.up = 783.66;
x7.lo = 11; x7.up = 11.01;
x8.lo = 49; x8.up = 49.01;
x9.lo = 2155.17; x9.up = 2155.18;
x10.lo = 252; x10.up = 252.01;
x11.up = 397.84;
x12.up = 222.32;
x13.up = 205.63;
x14.up = 205.63;
x15.up = 205.63;
x16.up = 124.83;
x17.up = 127.01;
x18.up = 61.08;
x19.up = 614.84;
x20.up = 778.08;
x21.up = 3024;
x22.up = 3024;
x23.up = 3024;
x24.up = 3024;
x25.up = 7257.6;
x26.up = 1209.6;
x27.up = 907.2;
x28.up = 7257.6;
x29.up = 7257.6;
x30.fx = 77;
x31.lo = 403.4; x31.up = 1312;
x32.fx = 158;
x33.fx = 16;
x34.fx = 0;
x35.lo = 502; x35.up = 928.46;
x36.fx = 11;
x37.fx = 49;
x38.lo = 915.3; x38.up = 2611.6;
x39.fx = 252;
x40.up = 397.84;
x41.up = 222.32;
x42.up = 205.63;
x43.up = 205.63;
x44.up = 205.63;
x45.up = 124.83;
x46.up = 127.01;
x47.up = 61.08;
x48.up = 614.84;
x49.up = 778.08;
x50.up = 3024;
x51.up = 3024;
x52.up = 3024;
x53.up = 3024;
x54.up = 7257.6;
x55.up = 1209.6;
x56.up = 907.2;
x57.up = 7257.6;
x58.up = 7257.6;
x59.fx = 77;
x60.lo = 403.4; x60.up = 1312;
x61.fx = 158;
x62.fx = 16;
x63.fx = 0;
x64.lo = 505.64; x64.up = 928.46;
x65.fx = 11;
x66.fx = 49;
x67.lo = 915.3; x67.up = 2611.6;
x68.fx = 252;
x69.up = 397.84;
x70.up = 222.32;
x71.up = 205.63;
x72.up = 205.63;
x73.up = 205.63;
x74.up = 124.83;
x75.up = 127.01;
x76.up = 61.08;
x77.up = 614.84;
x78.up = 778.08;
x79.up = 3024;
x80.up = 3024;
x81.up = 3024;
x82.up = 3024;
x83.up = 7257.6;
x84.up = 1209.6;
x85.up = 907.2;
x86.up = 7257.6;
x87.up = 7257.6;
x88.lo = 77; x88.up = 77.01;
x89.lo = 1100; x89.up = 1100.01;
x90.lo = 158; x90.up = 158.01;
x91.lo = 16; x91.up = 16.01;
x92.fx = 0;
x93.lo = 700; x93.up = 700.01;
x94.lo = 11; x94.up = 11.01;
x95.lo = 49; x95.up = 49.01;
x96.lo = 2000; x96.up = 2000.01;
x97.lo = 252; x97.up = 252.01;
* set non default levels
x11.l = 51.38;
x12.l = 140.21;
x13.l = 142.79;
x14.l = 21.91;
x15.l = 164.7;
x16.l = 58.19;
x17.l = 58.19;
x18.l = 61.08;
x19.l = 566.43;
x20.l = 583.57;
x31.l = 1049.53;
x35.l = 738.43;
x38.l = 1835.36;
x40.l = 32.06;
x42.l = 4.2;
x43.l = 48.37;
x44.l = 52.57;
x45.l = 59.85;
x46.l = 59.85;
x47.l = 58.24;
x49.l = 18.76;
x60.l = 1081.87;
x64.l = 696.71;
x67.l = 1972.77;
x69.l = 18.13;
x72.l = 18.13;
x73.l = 18.13;
x74.l = 5.81;
x75.l = 5.81;
x78.l = 6.02;
* 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;