# MINLP written by GAMS Convert at 02/17/22 17:23:52 # # Equation counts # Total E G L N X C B # 137 53 14 70 0 0 0 0 # # Variable counts # x b i s1s s2s sc si # Total cont binary integer sos1 sos2 scont sint # 195 69 28 0 98 0 0 0 # FX 0 # # Nonzero counts # Total const NL # 753 707 46 # # Reformulation has removed 1 variable and 1 equation from pyomo.environ import * model = m = ConcreteModel() m.x1 = Var(within=Reals, bounds=(0,None), initialize=0) m.x2 = Var(within=Reals, bounds=(0,None), initialize=0) m.x3 = Var(within=Reals, bounds=(0,None), initialize=0) m.x4 = Var(within=Reals, bounds=(0,None), initialize=0) m.x5 = Var(within=Reals, bounds=(0,None), initialize=0) m.x6 = Var(within=Reals, bounds=(0,None), initialize=0) m.x7 = Var(within=Reals, bounds=(0,None), initialize=0) m.x8 = Var(within=Reals, bounds=(0,None), initialize=0) m.x9 = Var(within=Reals, bounds=(0,None), initialize=0) m.x10 = Var(within=Reals, bounds=(0,None), initialize=0) m.x11 = Var(within=Reals, bounds=(0,None), initialize=0) m.x12 = Var(within=Reals, bounds=(0,None), initialize=0) m.x13 = Var(within=Reals, bounds=(0,None), initialize=0) m.x14 = Var(within=Reals, bounds=(0,None), initialize=0) m.x15 = Var(within=Reals, bounds=(0,None), initialize=0) m.x16 = Var(within=Reals, bounds=(0,None), initialize=0) m.x17 = Var(within=Reals, bounds=(0,None), initialize=0) m.x18 = Var(within=Reals, bounds=(0,None), initialize=0) m.x19 = Var(within=Reals, bounds=(0,None), initialize=0) m.x20 = Var(within=Reals, bounds=(0,None), initialize=0) m.x21 = Var(within=Reals, bounds=(0,None), initialize=0) m.x22 = Var(within=Reals, bounds=(0,None), initialize=0) m.x23 = Var(within=Reals, bounds=(0,None), initialize=0) m.x24 = Var(within=Reals, bounds=(0,None), initialize=0) m.x25 = Var(within=Reals, bounds=(0,None), initialize=0) m.x26 = Var(within=Reals, bounds=(0,None), initialize=0) m.x27 = Var(within=Reals, bounds=(0,None), initialize=0) m.x28 = Var(within=Reals, bounds=(0,None), initialize=0) m.x29 = Var(within=Reals, bounds=(None,None), initialize=0) m.x30 = Var(within=Reals, bounds=(None,None), initialize=0) m.x31 = Var(within=Reals, bounds=(None,None), initialize=0) m.x32 = Var(within=Reals, bounds=(None,None), initialize=0) m.x33 = Var(within=Reals, bounds=(None,None), initialize=0) m.x34 = Var(within=Reals, bounds=(None,None), initialize=0) m.x35 = Var(within=Reals, bounds=(None,None), initialize=0) m.x36 = Var(within=Reals, bounds=(None,None), initialize=0) m.x37 = Var(within=Reals, bounds=(None,None), initialize=0) m.x38 = Var(within=Reals, bounds=(None,None), initialize=0) m.x39 = Var(within=Reals, bounds=(None,None), initialize=0) m.x40 = Var(within=Reals, bounds=(None,None), initialize=0) m.x41 = Var(within=Reals, bounds=(None,None), initialize=0) m.x42 = Var(within=Reals, bounds=(None,None), initialize=0) m.x43 = Var(within=Reals, bounds=(6.5,None), initialize=11.5) m.x44 = Var(within=Reals, bounds=(3.25,None), initialize=8.25) m.x45 = Var(within=Reals, bounds=(16.58,None), initialize=21.58) m.x46 = Var(within=Reals, bounds=(14.92,None), initialize=19.92) m.x47 = Var(within=Reals, bounds=(12.925,None), initialize=17.925) m.x48 = Var(within=Reals, bounds=(12.26,None), initialize=17.26) m.x49 = Var(within=Reals, bounds=(8.76,None), initialize=13.76) m.x50 = Var(within=Reals, bounds=(16.08,None), initialize=21.08) m.x51 = Var(within=Reals, bounds=(None,None), initialize=0) m.x52 = Var(within=Reals, bounds=(None,None), initialize=0) m.x53 = Var(within=Reals, bounds=(None,None), initialize=0) m.x54 = Var(within=Reals, bounds=(None,None), initialize=0) m.x55 = Var(within=Reals, bounds=(None,None), initialize=0) m.x56 = Var(within=Reals, bounds=(None,None), initialize=0) m.x57 = Var(within=Reals, bounds=(None,None), initialize=0) m.x58 = Var(within=Reals, bounds=(None,None), initialize=0) m.x59 = Var(within=Reals, bounds=(None,None), initialize=0) m.x60 = Var(within=Reals, bounds=(None,None), initialize=0) m.x61 = Var(within=Reals, bounds=(None,None), initialize=0) m.x62 = Var(within=Reals, bounds=(None,None), initialize=0) m.x63 = Var(within=Reals, bounds=(None,None), initialize=0) m.x64 = Var(within=Reals, bounds=(None,None), initialize=0) m.x65 = Var(within=Reals, bounds=(0,2.5), initialize=0.961470588235294) m.x66 = Var(within=Reals, bounds=(0,6), initialize=2.30752941176471) m.x67 = Var(within=Reals, bounds=(None,None), initialize=0) m.x68 = Var(within=Reals, bounds=(None,None), initialize=0) m.x69 = Var(within=Reals, bounds=(None,None), initialize=0) m.b70 = Var(within=Binary, bounds=(0,1), initialize=0) m.b71 = Var(within=Binary, bounds=(0,1), initialize=0) m.b72 = Var(within=Binary, bounds=(0,1), initialize=0) m.b73 = Var(within=Binary, bounds=(0,1), initialize=0) m.b74 = Var(within=Binary, bounds=(0,1), initialize=0) m.b75 = Var(within=Binary, bounds=(0,1), initialize=0) m.b76 = Var(within=Binary, bounds=(0,1), initialize=0) m.b77 = Var(within=Binary, bounds=(0,1), initialize=0) m.b78 = Var(within=Binary, bounds=(0,1), initialize=0) m.b79 = Var(within=Binary, bounds=(0,1), initialize=0) m.b80 = Var(within=Binary, bounds=(0,1), initialize=0) m.b81 = Var(within=Binary, bounds=(0,1), initialize=0) m.b82 = Var(within=Binary, bounds=(0,1), initialize=0) m.b83 = Var(within=Binary, bounds=(0,1), initialize=0) m.b84 = Var(within=Binary, bounds=(0,1), initialize=0) m.b85 = Var(within=Binary, bounds=(0,1), initialize=0) m.b86 = Var(within=Binary, bounds=(0,1), initialize=0) m.b87 = Var(within=Binary, bounds=(0,1), initialize=0) m.b88 = Var(within=Binary, bounds=(0,1), initialize=0) m.b89 = Var(within=Binary, bounds=(0,1), initialize=0) m.b90 = Var(within=Binary, bounds=(0,1), initialize=0) m.b91 = Var(within=Binary, bounds=(0,1), initialize=0) m.b92 = Var(within=Binary, bounds=(0,1), initialize=0) m.b93 = Var(within=Binary, bounds=(0,1), initialize=0) m.b94 = Var(within=Binary, bounds=(0,1), initialize=0) m.b95 = Var(within=Binary, bounds=(0,1), initialize=0) m.b96 = Var(within=Binary, bounds=(0,1), initialize=0) m.b97 = Var(within=Binary, bounds=(0,1), initialize=0) m.s1s1 = Var([98,99,100,101,102,103,104], within=Reals, bounds=(0,None)) m.s1s2 = Var([105,106,107,108,109,110,111], within=Reals, bounds=(0,None)) m.s1s3 = Var([112,113,114,115,116,117,118], within=Reals, bounds=(0,None)) m.s1s4 = Var([119,120,121,122,123,124,125], within=Reals, bounds=(0,None)) m.s1s5 = Var([126,127,128,129,130,131,132], within=Reals, bounds=(0,None)) m.s1s6 = Var([133,134,135,136,137,138,139], within=Reals, bounds=(0,None)) m.s1s7 = Var([140,141,142,143,144,145,146], within=Reals, bounds=(0,None)) m.s1s8 = Var([147,148,149,150,151,152,153], within=Reals, bounds=(0,None)) m.s1s9 = Var([154,155,156,157,158,159,160], within=Reals, bounds=(0,None)) m.s1s10 = Var([161,162,163,164,165,166,167], within=Reals, bounds=(0,None)) m.s1s11 = Var([168,169,170,171,172,173,174], within=Reals, bounds=(0,None)) m.s1s12 = Var([175,176,177,178,179,180,181], within=Reals, bounds=(0,None)) m.s1s13 = Var([182,183,184,185,186,187,188], within=Reals, bounds=(0,None)) m.s1s14 = Var([189,190,191,192,193,194,195], within=Reals, bounds=(0,None)) m.obj = Objective(sense=minimize, expr= 10 * m.x67 + m.x68 + 10 * m.x69) m.es1s1 = SOSConstraint(sos=1, var=m.s1s1) m.es1s2 = SOSConstraint(sos=1, var=m.s1s2) m.es1s3 = SOSConstraint(sos=1, var=m.s1s3) m.es1s4 = SOSConstraint(sos=1, var=m.s1s4) m.es1s5 = SOSConstraint(sos=1, var=m.s1s5) m.es1s6 = SOSConstraint(sos=1, var=m.s1s6) m.es1s7 = SOSConstraint(sos=1, var=m.s1s7) m.es1s8 = SOSConstraint(sos=1, var=m.s1s8) m.es1s9 = SOSConstraint(sos=1, var=m.s1s9) m.es1s10 = SOSConstraint(sos=1, var=m.s1s10) m.es1s11 = SOSConstraint(sos=1, var=m.s1s11) m.es1s12 = SOSConstraint(sos=1, var=m.s1s12) m.es1s13 = SOSConstraint(sos=1, var=m.s1s13) m.es1s14 = SOSConstraint(sos=1, var=m.s1s14) m.e1 = Constraint(expr= -m.x1 - m.x2 - m.x3 + m.x15 + m.x16 + m.x17 + m.x65 == 0) m.e2 = Constraint(expr= -m.x4 - m.x5 - m.x6 - m.x7 + m.x18 + m.x19 + m.x20 + m.x21 + m.x66 == 0) m.e3 = Constraint(expr= m.x1 + m.x4 - m.x8 - m.x9 - m.x10 - m.x11 - m.x15 - m.x18 + m.x22 + m.x23 + m.x24 + m.x25 == 1.212) m.e4 = Constraint(expr= m.x2 + m.x8 + m.x12 - m.x16 - m.x22 - m.x26 == 0.452) m.e5 = Constraint(expr= m.x9 - m.x12 + m.x13 - m.x23 + m.x26 - m.x27 == 0.245) m.e6 = Constraint(expr= m.x5 + m.x10 - m.x13 - m.x14 - m.x19 - m.x24 + m.x27 + m.x28 == 0.652) m.e7 = Constraint(expr= m.x6 + m.x14 - m.x20 - m.x28 == 0.252) m.e8 = Constraint(expr= m.x3 + m.x7 + m.x11 - m.x17 - m.x21 - m.x25 == 0.456) m.e9 = Constraint(expr= m.x29 - 38721.1970117411 * m.s1s1[98] - 2543.8701482414 * m.s1s1[99] - 207.747320703761 * m.s1s1[100] - 23.9314504121258 * m.s1s1[101] - 1.5722267648148 * m.s1s1[102] - 0.181112645550961 * m.s1s1[103] - 0.0390863672545667 * m.s1s1[104] == 0) m.e10 = Constraint(expr= m.x30 - 32510.4890865135 * m.s1s2[105] - 2135.84468132099 * m.s1s2[106] - 174.425573683688 * m.s1s2[107] - 20.0929521164322 * m.s1s2[108] - 1.32004857865156 * m.s1s2[109] - 0.152062982061963 * m.s1s2[110] - 0.0328170876451919 * m.s1s2[111] == 0) m.e11 = Constraint(expr= m.x31 - 63468.4628982673 * m.s1s3[112] - 4169.69361956223 * m.s1s3[113] - 340.521578201805 * m.s1s3[114] - 39.2263796008983 * m.s1s3[115] - 2.57705917665854 * m.s1s3[116] - 0.296864304610023 * m.s1s3[117] - 0.0640670186196026 * m.s1s3[118] == 0) m.e12 = Constraint(expr= m.x32 - 50797.5773435889 * m.s1s4[119] - 3337.25325093014 * m.s1s4[120] - 272.539627020641 * m.s1s4[121] - 31.3951994533022 * m.s1s4[122] - 2.06257339263358 * m.s1s4[123] - 0.237598120158509 * m.s1s4[124] - 0.0512766370081929 * m.s1s4[125] == 0) m.e13 = Constraint(expr= m.x33 - 59165.7349698592 * m.s1s5[126] - 3887.01689524085 * m.s1s5[127] - 317.436542928413 * m.s1s5[128] - 36.5670992066393 * m.s1s5[129] - 2.40235218067626 * m.s1s5[130] - 0.27673893405488 * m.s1s5[131] - 0.0597237127048799 * m.s1s5[132] == 0) m.e14 = Constraint(expr= m.x34 - 32977.2294678044 * m.s1s6[133] - 2166.50816836621 * m.s1s6[134] - 176.929733450444 * m.s1s6[135] - 20.3814187742893 * m.s1s6[136] - 1.339 * m.s1s6[137] - 0.154246090843839 * m.s1s6[138] - 0.0332882297421199 * m.s1s6[139] == 0) m.e15 = Constraint(expr= m.x35 - 33843.9321019273 * m.s1s7[140] - 2223.4480134252 * m.s1s7[141] - 181.579774357788 * m.s1s7[142] - 20.9170801874496 * m.s1s7[143] - 1.37419139860501 * m.s1s7[144] - 0.158299963634093 * m.s1s7[145] - 0.0341631060391402 * m.s1s7[146] == 0) m.e16 = Constraint(expr= m.x36 - 31810.181054648 * m.s1s8[147] - 2089.8364782095 * m.s1s8[148] - 170.668274619734 * m.s1s8[149] - 19.660130090483 * m.s1s8[150] - 1.2916134290104 * m.s1s8[151] - 0.148787395299671 * m.s1s8[152] - 0.0321101751776739 * m.s1s8[153] == 0) m.e17 = Constraint(expr= m.x37 - 39461.9459070343 * m.s1s9[154] - 2592.53519858857 * m.s1s9[155] - 211.721593458417 * m.s1s9[156] - 24.3892667200816 * m.s1s9[157] - 1.60230396616872 * m.s1s9[158] - 0.184577388442944 * m.s1s9[159] - 0.0398341019735132 * m.s1s9[160] == 0) m.e18 = Constraint(expr= m.x38 - 32977.2294678044 * m.s1s10[161] - 2166.50816836621 * m.s1s10[162] - 176.929733450444 * m.s1s10[163] - 20.3814187742893 * m.s1s10[164] - 1.339 * m.s1s10[165] - 0.154246090843839 * m.s1s10[166] - 0.0332882297421199 * m.s1s10[167] == 0) m.e19 = Constraint(expr= m.x39 - 52785.5148814787 * m.s1s11[168] - 3467.85497167945 * m.s1s11[169] - 283.205327698691 * m.s1s11[170] - 32.6238347301504 * m.s1s11[171] - 2.14329116080854 * m.s1s11[172] - 0.246896402610059 * m.s1s11[173] - 0.0532833223041444 * m.s1s11[174] == 0) m.e20 = Constraint(expr= m.x40 - 30677.4142839491 * m.s1s12[175] - 2015.41699236491 * m.s1s12[176] - 164.590743970989 * m.s1s12[177] - 18.9600290116536 * m.s1s12[178] - 1.24561882211213 * m.s1s12[179] - 0.143489047044288 * m.s1s12[180] - 0.0309667255575633 * m.s1s12[181] == 0) m.e21 = Constraint(expr= m.x41 - 28361.2795383154 * m.s1s13[182] - 1863.25366856746 * m.s1s13[183] - 152.164196629274 * m.s1s13[184] - 17.5285530220005 * m.s1s13[185] - 1.15157500841239 * m.s1s13[186] - 0.132655670919396 * m.s1s13[187] - 0.0286287479053886 * m.s1s13[188] == 0) m.e22 = Constraint(expr= m.x42 - 50797.5773435889 * m.s1s14[189] - 3337.25325093014 * m.s1s14[190] - 272.539627020641 * m.s1s14[191] - 31.3951994533022 * m.s1s14[192] - 2.06257339263358 * m.s1s14[193] - 0.237598120158509 * m.s1s14[194] - 0.0512766370081929 * m.s1s14[195] == 0) m.e23 = Constraint(expr= -m.x29 * (m.x1 + m.x15) * (m.x1 - m.x15) + m.x43 - m.x45 - m.x51 == 0) m.e24 = Constraint(expr= -m.x30 * (m.x2 + m.x16) * (m.x2 - m.x16) + m.x43 - m.x46 - m.x52 == 0) m.e25 = Constraint(expr= -m.x31 * (m.x3 + m.x17) * (m.x3 - m.x17) + m.x43 - m.x50 - m.x53 == 0) m.e26 = Constraint(expr= -m.x32 * (m.x4 + m.x18) * (m.x4 - m.x18) + m.x44 - m.x45 - m.x54 == 0) m.e27 = Constraint(expr= -m.x33 * (m.x5 + m.x19) * (m.x5 - m.x19) + m.x44 - m.x48 - m.x55 == 0) m.e28 = Constraint(expr= -m.x34 * (m.x6 + m.x20) * (m.x6 - m.x20) + m.x44 - m.x49 - m.x56 == 0) m.e29 = Constraint(expr= -m.x35 * (m.x7 + m.x21) * (m.x7 - m.x21) + m.x44 - m.x50 - m.x57 == 0) m.e30 = Constraint(expr= -m.x36 * (m.x8 + m.x22) * (m.x8 - m.x22) + m.x45 - m.x46 - m.x58 == 0) m.e31 = Constraint(expr= -m.x37 * (m.x9 + m.x23) * (m.x9 - m.x23) + m.x45 - m.x47 - m.x59 == 0) m.e32 = Constraint(expr= -m.x38 * (m.x10 + m.x24) * (m.x10 - m.x24) + m.x45 - m.x48 - m.x60 == 0) m.e33 = Constraint(expr= -m.x39 * (m.x11 + m.x25) * (m.x11 - m.x25) + m.x45 - m.x50 - m.x61 == 0) m.e34 = Constraint(expr= -m.x40 * (m.x12 + m.x26) * (m.x12 - m.x26) - m.x46 + m.x47 - m.x62 == 0) m.e35 = Constraint(expr= -m.x41 * (m.x13 + m.x27) * (m.x13 - m.x27) - m.x47 + m.x48 - m.x63 == 0) m.e36 = Constraint(expr= -m.x42 * (m.x14 + m.x28) * (m.x14 - m.x28) + m.x48 - m.x49 - m.x64 == 0) m.e37 = Constraint(expr= m.x51 + 12 * m.b84 <= 12) m.e38 = Constraint(expr= m.x52 + 12 * m.b85 <= 12) m.e39 = Constraint(expr= m.x53 + 12 * m.b86 <= 12) m.e40 = Constraint(expr= m.x54 + 12 * m.b87 <= 12) m.e41 = Constraint(expr= m.x55 + 12 * m.b88 <= 12) m.e42 = Constraint(expr= m.x56 + 12 * m.b89 <= 12) m.e43 = Constraint(expr= m.x57 + 12 * m.b90 <= 12) m.e44 = Constraint(expr= m.x58 + 12 * m.b91 <= 12) m.e45 = Constraint(expr= m.x59 + 12 * m.b92 <= 12) m.e46 = Constraint(expr= m.x60 + 12 * m.b93 <= 12) m.e47 = Constraint(expr= m.x61 + 12 * m.b94 <= 12) m.e48 = Constraint(expr= m.x62 + 12 * m.b95 <= 12) m.e49 = Constraint(expr= m.x63 + 12 * m.b96 <= 12) m.e50 = Constraint(expr= m.x64 + 12 * m.b97 <= 12) m.e51 = Constraint(expr= m.x51 - 12 * m.b84 >= -12) m.e52 = Constraint(expr= m.x52 - 12 * m.b85 >= -12) m.e53 = Constraint(expr= m.x53 - 12 * m.b86 >= -12) m.e54 = Constraint(expr= m.x54 - 12 * m.b87 >= -12) m.e55 = Constraint(expr= m.x55 - 12 * m.b88 >= -12) m.e56 = Constraint(expr= m.x56 - 12 * m.b89 >= -12) m.e57 = Constraint(expr= m.x57 - 12 * m.b90 >= -12) m.e58 = Constraint(expr= m.x58 - 12 * m.b91 >= -12) m.e59 = Constraint(expr= m.x59 - 12 * m.b92 >= -12) m.e60 = Constraint(expr= m.x60 - 12 * m.b93 >= -12) m.e61 = Constraint(expr= m.x61 - 12 * m.b94 >= -12) m.e62 = Constraint(expr= m.x62 - 12 * m.b95 >= -12) m.e63 = Constraint(expr= m.x63 - 12 * m.b96 >= -12) m.e64 = Constraint(expr= m.x64 - 12 * m.b97 >= -12) m.e65 = Constraint(expr= -1.02 * m.x65 * (-6.5 + m.x43) - 1.02 * m.x66 * (-3.25 + m.x44) + m.x67 == 0) m.e66 = Constraint(expr= m.x68 - 9.11349113439539 * m.s1s1[98] - 17.6144733325531 * m.s1s1[99] - 32.2986551864818 * m.s1s1[100] - 54.4931814987685 * m.s1s1[101] - 105.323928905069 * m.s1s1[102] - 177.698914733437 * m.s1s1[103] - 257.546555368226 * m.s1s1[104] - 7.65172765642961 * m.s1s2[105] - 14.7891900880288 * m.s1s2[106] - 27.118094428506 * m.s1s2[107] - 45.7527173518919 * m.s1s2[108] - 88.4304387640365 * m.s1s2[109] - 149.196798497086 * m.s1s2[110] - 216.237232413786 * m.s1s2[111] - 14.9380525029139 * m.s1s3[112] - 28.8721329260735 * m.s1s3[113] - 52.941183552398 * m.s1s3[114] - 89.3205462402005 * m.s1s3[115] - 172.637944844116 * m.s1s3[116] - 291.268810037089 * m.s1s3[117] - 422.148209648796 * m.s1s3[118] - 11.9558099050809 * m.s1s4[119] - 23.1080813747994 * m.s1s4[120] - 42.3719709499612 * m.s1s4[121] - 71.4885338137291 * m.s1s4[122] - 138.172392322055 * m.s1s4[123] - 233.119713791557 * m.s1s4[124] - 337.870264236031 * m.s1s4[125] - 13.9253546563734 * m.s1s5[126] - 26.9147996770731 * m.s1s5[127] - 49.3521332015331 * m.s1s5[128] - 83.2652237802191 * m.s1s5[129] - 160.93427229773 * m.s1s5[130] - 271.522775764452 * m.s1s5[131] - 393.529446744536 * m.s1s5[132] - 7.76158051882097 * m.s1s6[133] - 15.0015127080393 * m.s1s6[134] - 27.5074183079396 * m.s1s6[135] - 46.4095712271164 * m.s1s6[136] - 89.7 * m.s1s6[137] - 151.338758602103 * m.s1s6[138] - 219.341665817957 * m.s1s6[139] - 7.96556922221359 * m.s1s7[140] - 15.3957802311063 * m.s1s7[141] - 28.2303641796868 * m.s1s7[142] - 47.6293006671023 * m.s1s7[143] - 92.0574820424717 * m.s1s7[144] - 155.316221319321 * m.s1s7[145] - 225.10637081608 * m.s1s7[146] - 7.48690188831565 * m.s1s8[147] - 14.4706163324673 * m.s1s8[148] - 26.5339439013751 * m.s1s8[149] - 44.7671586494086 * m.s1s8[150] - 86.5255598074927 * m.s1s8[151] - 145.982952158506 * m.s1s8[152] - 211.579268940989 * m.s1s8[153] - 9.28783513744935 * m.s1s9[154] - 17.9514438466182 * m.s1s9[155] - 32.916538800503 * m.s1s9[156] - 55.5356535066454 * m.s1s9[157] - 107.338809384118 * m.s1s9[158] - 181.098351861986 * m.s1s9[159] - 262.473503425068 * m.s1s9[160] - 7.76158051882097 * m.s1s10[161] - 15.0015127080393 * m.s1s10[162] - 27.5074183079396 * m.s1s10[163] - 46.4095712271164 * m.s1s10[164] - 89.7 * m.s1s10[165] - 151.338758602103 * m.s1s10[166] - 219.341665817957 * m.s1s10[167] - 12.4236944883441 * m.s1s11[168] - 24.0124044704238 * m.s1s11[169] - 44.0301766363479 * m.s1s11[170] - 74.2862014846846 * m.s1s11[171] - 143.579699122125 * m.s1s11[172] - 242.242736071415 * m.s1s11[173] - 351.092646411238 * m.s1s11[174] - 7.22029184733547 * m.s1s12[175] - 13.9553148538372 * m.s1s12[176] - 25.5890649679471 * m.s1s12[177] - 43.1729913716576 * m.s1s12[178] - 83.44436769489 * m.s1s12[179] - 140.784470672041 * m.s1s12[180] - 204.044889780639 * m.s1s12[181] - 6.67516217420068 * m.s1s13[182] - 12.9016931463472 * m.s1s13[183] - 23.6570989315674 * m.s1s13[184] - 39.913444642481 * m.s1s13[185] - 77.1443452237428 * m.s1s13[186] - 130.155289178744 * m.s1s13[187] - 188.639567333459 * m.s1s13[188] - 11.9558099050809 * m.s1s14[189] - 23.1080813747994 * m.s1s14[190] - 42.3719709499612 * m.s1s14[191] - 71.4885338137291 * m.s1s14[192] - 138.172392322055 * m.s1s14[193] - 233.119713791557 * m.s1s14[194] - 337.870264236031 * m.s1s14[195] == 0) m.e67 = Constraint(expr= -0.2 * m.x65 - 0.17 * m.x66 + m.x69 == 0) m.e68 = Constraint(expr= m.x1 - 2 * m.b70 <= 0) m.e69 = Constraint(expr= m.x2 - 2 * m.b71 <= 0) m.e70 = Constraint(expr= m.x3 - 2 * m.b72 <= 0) m.e71 = Constraint(expr= m.x4 - 2 * m.b73 <= 0) m.e72 = Constraint(expr= m.x5 - 2 * m.b74 <= 0) m.e73 = Constraint(expr= m.x6 - 2 * m.b75 <= 0) m.e74 = Constraint(expr= m.x7 - 2 * m.b76 <= 0) m.e75 = Constraint(expr= m.x8 - 2 * m.b77 <= 0) m.e76 = Constraint(expr= m.x9 - 2 * m.b78 <= 0) m.e77 = Constraint(expr= m.x10 - 2 * m.b79 <= 0) m.e78 = Constraint(expr= m.x11 - 2 * m.b80 <= 0) m.e79 = Constraint(expr= m.x12 - 2 * m.b81 <= 0) m.e80 = Constraint(expr= m.x13 - 2 * m.b82 <= 0) m.e81 = Constraint(expr= m.x14 - 2 * m.b83 <= 0) m.e82 = Constraint(expr= m.x15 + 2 * m.b70 <= 2) m.e83 = Constraint(expr= m.x16 + 2 * m.b71 <= 2) m.e84 = Constraint(expr= m.x17 + 2 * m.b72 <= 2) m.e85 = Constraint(expr= m.x18 + 2 * m.b73 <= 2) m.e86 = Constraint(expr= m.x19 + 2 * m.b74 <= 2) m.e87 = Constraint(expr= m.x20 + 2 * m.b75 <= 2) m.e88 = Constraint(expr= m.x21 + 2 * m.b76 <= 2) m.e89 = Constraint(expr= m.x22 + 2 * m.b77 <= 2) m.e90 = Constraint(expr= m.x23 + 2 * m.b78 <= 2) m.e91 = Constraint(expr= m.x24 + 2 * m.b79 <= 2) m.e92 = Constraint(expr= m.x25 + 2 * m.b80 <= 2) m.e93 = Constraint(expr= m.x26 + 2 * m.b81 <= 2) m.e94 = Constraint(expr= m.x27 + 2 * m.b82 <= 2) m.e95 = Constraint(expr= m.x28 + 2 * m.b83 <= 2) m.e96 = Constraint(expr= m.x1 + m.x15 - 2 * m.b84 <= 0) m.e97 = Constraint(expr= m.x2 + m.x16 - 2 * m.b85 <= 0) m.e98 = Constraint(expr= m.x3 + m.x17 - 2 * m.b86 <= 0) m.e99 = Constraint(expr= m.x4 + m.x18 - 2 * m.b87 <= 0) m.e100 = Constraint(expr= m.x5 + m.x19 - 2 * m.b88 <= 0) m.e101 = Constraint(expr= m.x6 + m.x20 - 2 * m.b89 <= 0) m.e102 = Constraint(expr= m.x7 + m.x21 - 2 * m.b90 <= 0) m.e103 = Constraint(expr= m.x8 + m.x22 - 2 * m.b91 <= 0) m.e104 = Constraint(expr= m.x9 + m.x23 - 2 * m.b92 <= 0) m.e105 = Constraint(expr= m.x10 + m.x24 - 2 * m.b93 <= 0) m.e106 = Constraint(expr= m.x11 + m.x25 - 2 * m.b94 <= 0) m.e107 = Constraint(expr= m.x12 + m.x26 - 2 * m.b95 <= 0) m.e108 = Constraint(expr= m.x13 + m.x27 - 2 * m.b96 <= 0) m.e109 = Constraint(expr= m.x14 + m.x28 - 2 * m.b97 <= 0) m.e110 = Constraint(expr= -m.b84 + m.s1s1[98] + m.s1s1[99] + m.s1s1[100] + m.s1s1[101] + m.s1s1[102] + m.s1s1[103] + m.s1s1[104] == 0) m.e111 = Constraint(expr= -m.b85 + m.s1s2[105] + m.s1s2[106] + m.s1s2[107] + m.s1s2[108] + m.s1s2[109] + m.s1s2[110] + m.s1s2[111] == 0) m.e112 = Constraint(expr= -m.b86 + m.s1s3[112] + m.s1s3[113] + m.s1s3[114] + m.s1s3[115] + m.s1s3[116] + m.s1s3[117] + m.s1s3[118] == 0) m.e113 = Constraint(expr= -m.b87 + m.s1s4[119] + m.s1s4[120] + m.s1s4[121] + m.s1s4[122] + m.s1s4[123] + m.s1s4[124] + m.s1s4[125] == 0) m.e114 = Constraint(expr= -m.b88 + m.s1s5[126] + m.s1s5[127] + m.s1s5[128] + m.s1s5[129] + m.s1s5[130] + m.s1s5[131] + m.s1s5[132] == 0) m.e115 = Constraint(expr= -m.b89 + m.s1s6[133] + m.s1s6[134] + m.s1s6[135] + m.s1s6[136] + m.s1s6[137] + m.s1s6[138] + m.s1s6[139] == 0) m.e116 = Constraint(expr= -m.b90 + m.s1s7[140] + m.s1s7[141] + m.s1s7[142] + m.s1s7[143] + m.s1s7[144] + m.s1s7[145] + m.s1s7[146] == 0) m.e117 = Constraint(expr= -m.b91 + m.s1s8[147] + m.s1s8[148] + m.s1s8[149] + m.s1s8[150] + m.s1s8[151] + m.s1s8[152] + m.s1s8[153] == 0) m.e118 = Constraint(expr= -m.b92 + m.s1s9[154] + m.s1s9[155] + m.s1s9[156] + m.s1s9[157] + m.s1s9[158] + m.s1s9[159] + m.s1s9[160] == 0) m.e119 = Constraint(expr= -m.b93 + m.s1s10[161] + m.s1s10[162] + m.s1s10[163] + m.s1s10[164] + m.s1s10[165] + m.s1s10[166] + m.s1s10[167] == 0) m.e120 = Constraint(expr= -m.b94 + m.s1s11[168] + m.s1s11[169] + m.s1s11[170] + m.s1s11[171] + m.s1s11[172] + m.s1s11[173] + m.s1s11[174] == 0) m.e121 = Constraint(expr= -m.b95 + m.s1s12[175] + m.s1s12[176] + m.s1s12[177] + m.s1s12[178] + m.s1s12[179] + m.s1s12[180] + m.s1s12[181] == 0) m.e122 = Constraint(expr= -m.b96 + m.s1s13[182] + m.s1s13[183] + m.s1s13[184] + m.s1s13[185] + m.s1s13[186] + m.s1s13[187] + m.s1s13[188] == 0) m.e123 = Constraint(expr= -m.b97 + m.s1s14[189] + m.s1s14[190] + m.s1s14[191] + m.s1s14[192] + m.s1s14[193] + m.s1s14[194] + m.s1s14[195] == 0) m.e124 = Constraint(expr= m.x1 + m.x15 - 0.0176041976445841 * m.s1s1[98] - 0.0686820348432157 * m.s1s1[99] - 0.240338257044582 * m.s1s1[100] - 0.708118780382974 * m.s1s1[101] - 2.76269556439349 * m.s1s1[102] - 8.13984688781727 * m.s1s1[103] - 17.521769741878 * m.s1s1[104] <= 0) m.e125 = Constraint(expr= m.x2 + m.x16 - 0.0192122784105588 * m.s1s2[105] - 0.0749558941482069 * m.s1s2[106] - 0.262292300976835 * m.s1s2[107] - 0.772802909347502 * m.s1s2[108] - 3.01505796619323 * m.s1s2[109] - 8.88339291488115 * m.s1s2[110] - 19.1223210124677 * m.s1s2[111] <= 0) m.e126 = Constraint(expr= m.x3 + m.x17 - 0.0137502828767635 * m.s1s3[112] - 0.0536461488738445 * m.s1s3[113] - 0.187723353667753 * m.s1s3[114] - 0.553097263345606 * m.s1s3[115] - 2.15788565203237 * m.s1s3[116] - 6.35786984108667 * m.s1s3[117] - 13.6859001084016 * m.s1s3[118] <= 0) m.e127 = Constraint(expr= m.x4 + m.x18 - 0.0153698320860398 * m.s1s4[119] - 0.0599647518268192 * m.s1s4[120] - 0.209833968534382 * m.s1s4[121] - 0.618242703881818 * m.s1s4[122] - 2.41204784147821 * m.s1s4[123] - 7.10671865867823 * m.s1s4[124] - 15.297866123752 * m.s1s4[125] <= 0) m.e128 = Constraint(expr= m.x5 + m.x19 - 0.0142414920290718 * m.s1s5[126] - 0.0555625806701283 * m.s1s5[127] - 0.194429501479406 * m.s1s5[128] - 0.572855870518057 * m.s1s5[129] - 2.23497302480958 * m.s1s5[130] - 6.58499563065165 * m.s1s5[131] - 14.1748092785675 * m.s1s5[132] <= 0) m.e129 = Constraint(expr= m.x6 + m.x20 - 0.0190758342372385 * m.s1s6[133] - 0.0744235629590588 * m.s1s6[134] - 0.260429520550158 * m.s1s6[135] - 0.767314520523847 * m.s1s6[136] - 2.99364524861134 * m.s1s6[137] - 8.82030371865732 * m.s1s6[138] - 18.9865157099025 * m.s1s6[139] <= 0) m.e130 = Constraint(expr= m.x7 + m.x21 - 0.0188299954674205 * m.s1s7[140] - 0.0734644333642121 * m.s1s7[141] - 0.257073249355929 * m.s1s7[142] - 0.757425796631457 * m.s1s7[143] - 2.95506480929543 * m.s1s7[144] - 8.70663253716935 * m.s1s7[145] - 18.7418280277177 * m.s1s7[146] <= 0) m.e131 = Constraint(expr= m.x8 + m.x22 - 0.0194226083350049 * m.s1s8[147] - 0.0757764874800376 * m.s1s8[148] - 0.265163793814297 * m.s1s8[149] - 0.781263310246409 * m.s1s8[150] - 3.04806586357418 * m.s1s8[151] - 8.98064548017735 * m.s1s8[152] - 19.3316661118795 * m.s1s8[153] <= 0) m.e132 = Constraint(expr= m.x9 + m.x23 - 0.0174381887671401 * m.s1s9[154] - 0.0680343582075014 * m.s1s9[155] - 0.238071849619242 * m.s1s9[156] - 0.701441168247406 * m.s1s9[157] - 2.73664314220282 * m.s1s9[158] - 8.06308753350335 * m.s1s9[159] - 17.3565381656137 * m.s1s9[160] <= 0) m.e133 = Constraint(expr= m.x10 + m.x24 - 0.0190758342372385 * m.s1s10[161] - 0.0744235629590588 * m.s1s10[162] - 0.260429520550158 * m.s1s10[163] - 0.767314520523847 * m.s1s10[164] - 2.99364524861134 * m.s1s10[165] - 8.82030371865732 * m.s1s10[166] - 18.9865157099025 * m.s1s10[167] <= 0) m.e134 = Constraint(expr= m.x11 + m.x25 - 0.0150776355652448 * m.s1s11[168] - 0.0588247594211735 * m.s1s11[169] - 0.205844806180028 * m.s1s11[170] - 0.606489265973719 * m.s1s11[171] - 2.36619229905421 * m.s1s11[172] - 6.97161253294371 * m.s1s11[173] - 15.0070377508768 * m.s1s11[174] <= 0) m.e135 = Constraint(expr= m.x12 + m.x26 - 0.0197779487583483 * m.s1s12[175] - 0.0771628331590627 * m.s1s12[176] - 0.270015017353593 * m.s1s12[177] - 0.795556675515238 * m.s1s12[178] - 3.10383082550202 * m.s1s12[179] - 9.14494814806731 * m.s1s12[180] - 19.6853427294401 * m.s1s12[181] <= 0) m.e136 = Constraint(expr= m.x13 + m.x27 - 0.02056968839856 * m.s1s13[182] - 0.0802517719822704 * m.s1s13[183] - 0.280824105561038 * m.s1s13[184] - 0.827403949655566 * m.s1s13[185] - 3.22808162274528 * m.s1s13[186] - 9.5110335315907 * m.s1s13[187] - 20.4733752175652 * m.s1s13[188] <= 0) m.e137 = Constraint(expr= m.x14 + m.x28 - 0.0153698320860398 * m.s1s14[189] - 0.0599647518268192 * m.s1s14[190] - 0.209833968534382 * m.s1s14[191] - 0.618242703881818 * m.s1s14[192] - 2.41204784147821 * m.s1s14[193] - 7.10671865867823 * m.s1s14[194] - 15.297866123752 * m.s1s14[195] <= 0)