spanhyd.gms
<![CDATA[
* NLP written by GAMS Convert at 10/06/06 11:52:25
*
* 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
* 249 195 54 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.. - ((438.00995511036 + 1.0862222572*((58.71308 + 0.1265*x1 - 0.00025225*
x1*x1)*x11 + (228.89241 + 2.57*(0.031*x2 - 7.2e-6*x2*x2 - (0.558*x3 -
0.0014052*x3*x3)))*x12 + 2.2146*(0.558*x3 - 0.0014052*x3*x3 - (0.0262*x9
- 2.4e-6*x9*x9)) + (122.3956928 + 2.5128*(0.7584*x4 - 0.014353*x4*x4 - (
0.0262*x9 - 2.4e-6*x9*x9)))*x14 + (84.6622463000001 + 2.5707*(0.0541*x6
- 2.01e-5*x6*x6 - (0.0014139*x7*x7 + 0.9038*x7)))*x16 + 2.7601*(
0.0014139*x7*x7 + 0.9038*x7 - (0.0017195*x8*x8 + 0.1557*x8)) + (
77.2387110000001 + 0.0017195*x8*x8 + 0.1557*x8)*x18 + (107.811243 +
0.0262*x9 - 2.4e-6*x9*x9)*x19 + 0.0002372*x10*x10 - 0.3111*x10) -
6.79942546339484*x13 - 9.69572848999292*x17 - 223.024240070561*x20)*(
403.2415578 + (58.71308 + 0.1265*x1 - 0.00025225*x1*x1)*x11 + (228.89241
+ 2.57*(0.031*x2 - 7.2e-6*x2*x2 - (0.558*x3 - 0.0014052*x3*x3)))*x12 +
2.2146*(0.558*x3 - 0.0014052*x3*x3 - (0.0262*x9 - 2.4e-6*x9*x9)) + (
122.3956928 + 2.5128*(0.7584*x4 - 0.014353*x4*x4 - (0.0262*x9 - 2.4e-6*x9
*x9)))*x14 + (84.6622463000001 + 2.5707*(0.0541*x6 - 2.01e-5*x6*x6 - (
0.0014139*x7*x7 + 0.9038*x7)))*x16 + 2.7601*(0.0014139*x7*x7 + 0.9038*x7
- (0.0017195*x8*x8 + 0.1557*x8)) + (77.2387110000001 + 0.0017195*x8*x8
+ 0.1557*x8)*x18 + (107.811243 + 0.0262*x9 - 2.4e-6*x9*x9)*x19 +
0.0002372*x10*x10 - 0.3111*x10 - 6.2597*x13 - 8.9261*x17 - 205.321*x20)
- 6.255270637*((58.71308 + 0.1265*x1 - 0.00025225*x1*x1)*x11 + (
228.89241 + 2.57*(0.031*x2 - 7.2e-6*x2*x2 - (0.558*x3 - 0.0014052*x3*x3))
)*x12 + 2.2146*(0.558*x3 - 0.0014052*x3*x3 - (0.0262*x9 - 2.4e-6*x9*x9))
+ (122.3956928 + 2.5128*(0.7584*x4 - 0.014353*x4*x4 - (0.0262*x9 -
2.4e-6*x9*x9)))*x14 + (84.6622463000001 + 2.5707*(0.0541*x6 - 2.01e-5*x6*
x6 - (0.0014139*x7*x7 + 0.9038*x7)))*x16 + 2.7601*(0.0014139*x7*x7 +
0.9038*x7 - (0.0017195*x8*x8 + 0.1557*x8)) + (77.2387110000001 +
0.0017195*x8*x8 + 0.1557*x8)*x18 + (107.811243 + 0.0262*x9 - 2.4e-6*x9*x9
)*x19 + 0.0002372*x10*x10 - 0.3111*x10) + 39.1561176064289*x13 +
55.8351712329257*x17 + 1284.33842245948*x20 - 6.021899007*((58.71308 +
0.1265*x30 - 0.00025225*x30*x30)*x40 + (228.89241 + 2.57*(0.031*x31 -
7.2e-6*x31*x31 - (0.558*x32 - 0.0014052*x32*x32)))*x41 + 2.2146*(0.558*
x32 - 0.0014052*x32*x32 - (0.0262*x38 - 2.4e-6*x38*x38)) + (122.3956928
+ 2.5128*(0.7584*x33 - 0.014353*x33*x33 - (0.0262*x38 - 2.4e-6*x38*x38))
)*x43 + (84.6622463000001 + 2.5707*(0.0541*x35 - 2.01e-5*x35*x35 - (
0.0014139*x36*x36 + 0.9038*x36)))*x45 + 2.7601*(0.0014139*x36*x36 +
0.9038*x36 - (0.0017195*x37*x37 + 0.1557*x37)) + (77.2387110000001 +
0.0017195*x37*x37 + 0.1557*x37)*x47 + (107.811243 + 0.0262*x38 - 2.4e-6*
x38*x38)*x48 + 0.0002372*x39*x39 - 0.3111*x39) + (3398.46649455427 +
8.4278676858*((58.71308 + 0.1265*x30 - 0.00025225*x30*x30)*x40 + (
228.89241 + 2.57*(0.031*x31 - 7.2e-6*x31*x31 - (0.558*x32 - 0.0014052*x32
*x32)))*x41 + 2.2146*(0.558*x32 - 0.0014052*x32*x32 - (0.0262*x38 -
2.4e-6*x38*x38)) + (122.3956928 + 2.5128*(0.7584*x33 - 0.014353*x33*x33
- (0.0262*x38 - 2.4e-6*x38*x38)))*x43 + (84.6622463000001 + 2.5707*(
0.0541*x35 - 2.01e-5*x35*x35 - (0.0014139*x36*x36 + 0.9038*x36)))*x45 +
2.7601*(0.0014139*x36*x36 + 0.9038*x36 - (0.0017195*x37*x37 + 0.1557*x37)
) + (77.2387110000001 + 0.0017195*x37*x37 + 0.1557*x37)*x47 + (107.811243
+ 0.0262*x38 - 2.4e-6*x38*x38)*x48 + 0.0002372*x39*x39 - 0.3111*x39) -
52.7559233528023*x42 - 75.2279897502194*x46 - 1730.41822111614*x49)*(
403.2415578 + (58.71308 + 0.1265*x30 - 0.00025225*x30*x30)*x40 + (
228.89241 + 2.57*(0.031*x31 - 7.2e-6*x31*x31 - (0.558*x32 - 0.0014052*x32
*x32)))*x41 + 2.2146*(0.558*x32 - 0.0014052*x32*x32 - (0.0262*x38 -
2.4e-6*x38*x38)) + (122.3956928 + 2.5128*(0.7584*x33 - 0.014353*x33*x33
- (0.0262*x38 - 2.4e-6*x38*x38)))*x43 + (84.6622463000001 + 2.5707*(
0.0541*x35 - 2.01e-5*x35*x35 - (0.0014139*x36*x36 + 0.9038*x36)))*x45 +
2.7601*(0.0014139*x36*x36 + 0.9038*x36 - (0.0017195*x37*x37 + 0.1557*x37)
) + (77.2387110000001 + 0.0017195*x37*x37 + 0.1557*x37)*x47 + (107.811243
+ 0.0262*x38 - 2.4e-6*x38*x38)*x48 + 0.0002372*x39*x39 - 0.3111*x39 -
6.2597*x42 - 8.9261*x46 - 205.321*x49) + 37.6952812141179*x42 +
53.7520727263827*x46 + 1236.42232601625*x49 - 6.08910295*((58.71308 +
0.1265*x59 - 0.00025225*x59*x59)*x69 + (228.89241 + 2.57*(0.031*x60 -
7.2e-6*x60*x60 - (0.558*x61 - 0.0014052*x61*x61)))*x70 + 2.2146*(0.558*
x61 - 0.0014052*x61*x61 - (0.0262*x67 - 2.4e-6*x67*x67)) + (122.3956928
+ 2.5128*(0.7584*x62 - 0.014353*x62*x62 - (0.0262*x67 - 2.4e-6*x67*x67))
)*x72 + (84.6622463000001 + 2.5707*(0.0541*x64 - 2.01e-5*x64*x64 - (
0.0014139*x65*x65 + 0.9038*x65)))*x74 + 2.7601*(0.0014139*x65*x65 +
0.9038*x65 - (0.0017195*x66*x66 + 0.1557*x66)) + (77.2387110000001 +
0.0017195*x66*x66 + 0.1557*x66)*x76 + (107.811243 + 0.0262*x67 - 2.4e-6*
x67*x67)*x77 + 0.0002372*x68*x68 - 0.3111*x68) + (3597.3022819255 +
8.9209611766*((58.71308 + 0.1265*x59 - 0.00025225*x59*x59)*x69 + (
228.89241 + 2.57*(0.031*x60 - 7.2e-6*x60*x60 - (0.558*x61 - 0.0014052*x61
*x61)))*x70 + 2.2146*(0.558*x61 - 0.0014052*x61*x61 - (0.0262*x67 -
2.4e-6*x67*x67)) + (122.3956928 + 2.5128*(0.7584*x62 - 0.014353*x62*x62
- (0.0262*x67 - 2.4e-6*x67*x67)))*x72 + (84.6622463000001 + 2.5707*(
0.0541*x64 - 2.01e-5*x64*x64 - (0.0014139*x65*x65 + 0.9038*x65)))*x74 +
2.7601*(0.0014139*x65*x65 + 0.9038*x65 - (0.0017195*x66*x66 + 0.1557*x66)
) + (77.2387110000001 + 0.0017195*x66*x66 + 0.1557*x66)*x76 + (107.811243
+ 0.0262*x67 - 2.4e-6*x67*x67)*x77 + 0.0002372*x68*x68 - 0.3111*x68) -
55.842540677163*x71 - 79.6293915584493*x75 - 1831.66066974069*x78)*(
403.2415578 + (58.71308 + 0.1265*x59 - 0.00025225*x59*x59)*x69 + (
228.89241 + 2.57*(0.031*x60 - 7.2e-6*x60*x60 - (0.558*x61 - 0.0014052*x61
*x61)))*x70 + 2.2146*(0.558*x61 - 0.0014052*x61*x61 - (0.0262*x67 -
2.4e-6*x67*x67)) + (122.3956928 + 2.5128*(0.7584*x62 - 0.014353*x62*x62
- (0.0262*x67 - 2.4e-6*x67*x67)))*x72 + (84.6622463000001 + 2.5707*(
0.0541*x64 - 2.01e-5*x64*x64 - (0.0014139*x65*x65 + 0.9038*x65)))*x74 +
2.7601*(0.0014139*x65*x65 + 0.9038*x65 - (0.0017195*x66*x66 + 0.1557*x66)
) + (77.2387110000001 + 0.0017195*x66*x66 + 0.1557*x66)*x76 + (107.811243
+ 0.0262*x67 - 2.4e-6*x67*x67)*x77 + 0.0002372*x68*x68 - 0.3111*x68 -
6.2597*x71 - 8.9261*x75 - 205.321*x78) + 38.115957736115*x71 +
54.351941841995*x75 + 1250.22070679695*x78) + objvar
=E= -7155.18601141601;
* 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;
]]></plaintext></body>