# MINLP written by GAMS Convert at 03/03/22 14:44:52 # # Equation counts # Total E G L N X C B # 467 198 40 229 0 0 0 0 # # Variable counts # x b i s1s s2s sc si # Total cont binary integer sos1 sos2 scont sint # 302 242 60 0 0 0 0 0 # FX 0 # # Nonzero counts # Total const NL # 1011 945 66 # # Reformulation has removed 1 variable and 1 equation from pyomo.environ import * model = m = ConcreteModel() m.x1 = Var(within=Reals, bounds=(None,None), initialize=0) m.x2 = Var(within=Reals, bounds=(None,None), initialize=0) m.x3 = Var(within=Reals, bounds=(None,None), initialize=0) m.x4 = Var(within=Reals, bounds=(None,None), initialize=0) m.x5 = Var(within=Reals, bounds=(None,None), initialize=0) m.x6 = Var(within=Reals, bounds=(None,None), initialize=0) m.x7 = Var(within=Reals, bounds=(None,None), initialize=0) m.x8 = Var(within=Reals, bounds=(None,None), initialize=0) m.x9 = Var(within=Reals, bounds=(None,None), initialize=0) m.x10 = Var(within=Reals, bounds=(None,None), initialize=0) m.x11 = Var(within=Reals, bounds=(None,None), initialize=0) m.x12 = Var(within=Reals, bounds=(None,None), initialize=0) m.x13 = Var(within=Reals, bounds=(None,None), initialize=0) m.x14 = Var(within=Reals, bounds=(None,None), initialize=0) m.x15 = Var(within=Reals, bounds=(None,None), initialize=0) m.x16 = Var(within=Reals, bounds=(None,None), initialize=0) m.x17 = Var(within=Reals, bounds=(None,None), initialize=0) m.x18 = Var(within=Reals, bounds=(None,None), initialize=0) m.x19 = Var(within=Reals, bounds=(None,None), initialize=0) m.x20 = Var(within=Reals, bounds=(None,None), initialize=0) m.x21 = Var(within=Reals, bounds=(None,None), initialize=0) m.x22 = Var(within=Reals, bounds=(None,None), initialize=0) m.x23 = Var(within=Reals, bounds=(None,None), initialize=0) m.x24 = Var(within=Reals, bounds=(None,None), initialize=0) m.x25 = Var(within=Reals, bounds=(None,None), initialize=0) m.x26 = Var(within=Reals, bounds=(None,None), initialize=0) m.x27 = Var(within=Reals, bounds=(None,None), initialize=0) m.x28 = Var(within=Reals, bounds=(None,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=(0,40), initialize=0) m.x32 = Var(within=Reals, bounds=(0,40), initialize=0) m.x33 = Var(within=Reals, bounds=(0,None), initialize=0) m.x34 = Var(within=Reals, bounds=(0,None), initialize=0) m.x35 = Var(within=Reals, bounds=(0,None), initialize=0) m.x36 = Var(within=Reals, bounds=(0,None), initialize=0) m.x37 = Var(within=Reals, bounds=(0,None), initialize=0) m.x38 = Var(within=Reals, bounds=(0,None), initialize=0) m.x39 = Var(within=Reals, bounds=(0,None), initialize=0) m.x40 = Var(within=Reals, bounds=(0,None), initialize=0) m.x41 = Var(within=Reals, bounds=(0,None), initialize=0) m.x42 = Var(within=Reals, bounds=(0,None), initialize=0) m.x43 = Var(within=Reals, bounds=(0,None), initialize=0) m.x44 = Var(within=Reals, bounds=(0,None), initialize=0) m.x45 = Var(within=Reals, bounds=(0,None), initialize=0) m.x46 = Var(within=Reals, bounds=(0,None), initialize=0) m.x47 = Var(within=Reals, bounds=(0,None), initialize=0) m.x48 = Var(within=Reals, bounds=(0,None), initialize=0) m.x49 = Var(within=Reals, bounds=(0,None), initialize=0) m.x50 = Var(within=Reals, bounds=(0,None), initialize=0) m.x51 = Var(within=Reals, bounds=(0,None), initialize=0) m.x52 = Var(within=Reals, bounds=(0,None), initialize=0) m.x53 = Var(within=Reals, bounds=(0,30), initialize=0) m.x54 = Var(within=Reals, bounds=(0,30), initialize=0) m.x55 = Var(within=Reals, bounds=(0,None), initialize=0) m.x56 = Var(within=Reals, bounds=(0,None), initialize=0) m.x57 = Var(within=Reals, bounds=(0,None), initialize=0) m.x58 = Var(within=Reals, bounds=(0,None), initialize=0) m.x59 = Var(within=Reals, bounds=(0,None), initialize=0) m.x60 = Var(within=Reals, bounds=(0,None), initialize=0) m.x61 = Var(within=Reals, bounds=(0,None), initialize=0) m.x62 = Var(within=Reals, bounds=(0,None), initialize=0) m.x63 = Var(within=Reals, bounds=(0,None), initialize=0) m.x64 = Var(within=Reals, bounds=(0,None), initialize=0) m.x65 = Var(within=Reals, bounds=(0,None), initialize=0) m.x66 = Var(within=Reals, bounds=(0,None), initialize=0) m.x67 = Var(within=Reals, bounds=(0,None), initialize=0) m.x68 = Var(within=Reals, bounds=(0,None), initialize=0) m.x69 = Var(within=Reals, bounds=(0,None), initialize=0) m.x70 = Var(within=Reals, bounds=(0,None), initialize=0) m.x71 = Var(within=Reals, bounds=(0,None), initialize=0) m.x72 = Var(within=Reals, bounds=(0,None), initialize=0) m.x73 = Var(within=Reals, bounds=(0,None), initialize=0) m.x74 = Var(within=Reals, bounds=(0,None), initialize=0) m.x75 = Var(within=Reals, bounds=(0,None), initialize=0) m.x76 = Var(within=Reals, bounds=(0,None), initialize=0) m.x77 = Var(within=Reals, bounds=(0,None), initialize=0) m.x78 = Var(within=Reals, bounds=(0,None), initialize=0) m.x79 = Var(within=Reals, bounds=(0,None), initialize=0) m.x80 = Var(within=Reals, bounds=(0,None), initialize=0) m.x81 = Var(within=Reals, bounds=(0,None), initialize=0) m.x82 = Var(within=Reals, bounds=(0,None), initialize=0) m.x83 = Var(within=Reals, bounds=(0,None), initialize=0) m.x84 = Var(within=Reals, bounds=(0,None), initialize=0) m.x85 = Var(within=Reals, bounds=(0,None), initialize=0) m.x86 = Var(within=Reals, bounds=(0,None), initialize=0) m.x87 = Var(within=Reals, bounds=(0,20), initialize=0) m.x88 = Var(within=Reals, bounds=(0,20), initialize=0) m.x89 = Var(within=Reals, bounds=(0,20), initialize=0) m.x90 = Var(within=Reals, bounds=(0,20), initialize=0) m.x91 = Var(within=Reals, bounds=(0,None), initialize=0) m.x92 = Var(within=Reals, bounds=(0,None), initialize=0) m.x93 = Var(within=Reals, bounds=(0,None), initialize=0) m.x94 = Var(within=Reals, bounds=(0,None), initialize=0) m.x95 = Var(within=Reals, bounds=(0,None), initialize=0) m.x96 = Var(within=Reals, bounds=(0,None), initialize=0) m.x97 = Var(within=Reals, bounds=(0,None), initialize=0) m.x98 = Var(within=Reals, bounds=(0,None), initialize=0) m.x99 = Var(within=Reals, bounds=(0,None), initialize=0) m.x100 = Var(within=Reals, bounds=(0,None), initialize=0) m.x101 = Var(within=Reals, bounds=(0,None), initialize=0) m.x102 = Var(within=Reals, bounds=(0,None), initialize=0) m.x103 = Var(within=Reals, bounds=(0,None), initialize=0) m.x104 = Var(within=Reals, bounds=(0,None), initialize=0) m.x105 = Var(within=Reals, bounds=(0,None), initialize=0) m.x106 = Var(within=Reals, bounds=(0,None), initialize=0) m.x107 = Var(within=Reals, bounds=(0,None), initialize=0) m.x108 = Var(within=Reals, bounds=(0,None), initialize=0) m.x109 = Var(within=Reals, bounds=(0,None), initialize=0) m.x110 = Var(within=Reals, bounds=(0,None), initialize=0) m.x111 = Var(within=Reals, bounds=(0,None), initialize=0) m.x112 = Var(within=Reals, bounds=(0,None), initialize=0) m.x113 = Var(within=Reals, bounds=(0,None), initialize=0) m.x114 = Var(within=Reals, bounds=(0,None), initialize=0) m.x115 = Var(within=Reals, bounds=(0,None), initialize=0) m.x116 = Var(within=Reals, bounds=(0,None), initialize=0) m.x117 = Var(within=Reals, bounds=(0,None), initialize=0) m.x118 = Var(within=Reals, bounds=(0,None), initialize=0) m.x119 = Var(within=Reals, bounds=(0,None), initialize=0) m.x120 = Var(within=Reals, bounds=(0,None), initialize=0) m.x121 = Var(within=Reals, bounds=(0,None), initialize=0) m.x122 = Var(within=Reals, bounds=(0,None), initialize=0) m.x123 = Var(within=Reals, bounds=(0,None), initialize=0) m.x124 = Var(within=Reals, bounds=(0,None), initialize=0) m.x125 = Var(within=Reals, bounds=(0,None), initialize=0) m.x126 = Var(within=Reals, bounds=(0,None), initialize=0) m.x127 = Var(within=Reals, bounds=(0,None), initialize=0) m.x128 = Var(within=Reals, bounds=(0,None), initialize=0) m.x129 = Var(within=Reals, bounds=(0,None), initialize=0) m.x130 = Var(within=Reals, bounds=(0,None), initialize=0) m.x131 = Var(within=Reals, bounds=(0,None), initialize=0) m.x132 = Var(within=Reals, bounds=(0,None), initialize=0) m.x133 = Var(within=Reals, bounds=(0,None), initialize=0) m.x134 = Var(within=Reals, bounds=(0,None), initialize=0) m.x135 = Var(within=Reals, bounds=(0,None), initialize=0) m.x136 = Var(within=Reals, bounds=(0,None), initialize=0) m.x137 = Var(within=Reals, bounds=(0,None), initialize=0) m.x138 = Var(within=Reals, bounds=(0,None), initialize=0) m.x139 = Var(within=Reals, bounds=(0,None), initialize=0) m.x140 = Var(within=Reals, bounds=(0,None), initialize=0) m.x141 = Var(within=Reals, bounds=(0,None), initialize=0) m.x142 = Var(within=Reals, bounds=(0,None), initialize=0) m.x143 = Var(within=Reals, bounds=(0,None), initialize=0) m.x144 = Var(within=Reals, bounds=(0,None), initialize=0) m.x145 = Var(within=Reals, bounds=(0,None), initialize=0) m.x146 = Var(within=Reals, bounds=(0,None), initialize=0) m.x147 = Var(within=Reals, bounds=(0,None), initialize=0) m.x148 = Var(within=Reals, bounds=(0,None), initialize=0) m.x149 = Var(within=Reals, bounds=(0,None), initialize=0) m.x150 = Var(within=Reals, bounds=(0,None), initialize=0) m.x151 = Var(within=Reals, bounds=(0,None), initialize=0) m.x152 = Var(within=Reals, bounds=(0,None), initialize=0) m.x153 = Var(within=Reals, bounds=(0,None), initialize=0) m.x154 = Var(within=Reals, bounds=(0,None), initialize=0) m.x155 = Var(within=Reals, bounds=(0,None), initialize=0) m.x156 = Var(within=Reals, bounds=(0,None), initialize=0) m.x157 = Var(within=Reals, bounds=(0,None), initialize=0) m.x158 = Var(within=Reals, bounds=(0,None), initialize=0) m.x159 = Var(within=Reals, bounds=(0,None), initialize=0) m.x160 = Var(within=Reals, bounds=(0,None), initialize=0) m.x161 = Var(within=Reals, bounds=(0,None), initialize=0) m.x162 = Var(within=Reals, bounds=(0,None), initialize=0) m.x163 = Var(within=Reals, bounds=(0,None), initialize=0) m.x164 = Var(within=Reals, bounds=(0,None), initialize=0) m.x165 = Var(within=Reals, bounds=(0,None), initialize=0) m.x166 = Var(within=Reals, bounds=(0,None), initialize=0) m.x167 = Var(within=Reals, bounds=(0,None), initialize=0) m.x168 = Var(within=Reals, bounds=(0,None), initialize=0) m.x169 = Var(within=Reals, bounds=(0,None), initialize=0) m.x170 = Var(within=Reals, bounds=(0,None), initialize=0) m.x171 = Var(within=Reals, bounds=(0,None), initialize=0) m.x172 = Var(within=Reals, bounds=(0,None), initialize=0) m.x173 = Var(within=Reals, bounds=(0,None), initialize=0) m.x174 = Var(within=Reals, bounds=(0,None), initialize=0) m.x175 = Var(within=Reals, bounds=(0,None), initialize=0) m.x176 = Var(within=Reals, bounds=(0,None), initialize=0) m.x177 = Var(within=Reals, bounds=(0,None), initialize=0) m.x178 = Var(within=Reals, bounds=(0,None), initialize=0) m.x179 = Var(within=Reals, bounds=(0,None), initialize=0) m.x180 = Var(within=Reals, bounds=(0,None), initialize=0) m.x181 = Var(within=Reals, bounds=(0,None), initialize=0) m.x182 = Var(within=Reals, bounds=(0,None), initialize=0) m.x183 = Var(within=Reals, bounds=(0,None), initialize=0) m.x184 = Var(within=Reals, bounds=(0,None), initialize=0) m.x185 = Var(within=Reals, bounds=(0,None), initialize=0) m.x186 = Var(within=Reals, bounds=(0,None), initialize=0) m.x187 = Var(within=Reals, bounds=(0,None), initialize=0) m.x188 = Var(within=Reals, bounds=(0,None), initialize=0) m.x189 = Var(within=Reals, bounds=(0,None), initialize=0) m.x190 = Var(within=Reals, bounds=(0,None), initialize=0) m.x191 = Var(within=Reals, bounds=(0,None), initialize=0) m.x192 = Var(within=Reals, bounds=(0,None), initialize=0) m.x193 = Var(within=Reals, bounds=(0,None), initialize=0) m.x194 = Var(within=Reals, bounds=(0,None), initialize=0) m.x195 = Var(within=Reals, bounds=(0,None), initialize=0) m.x196 = Var(within=Reals, bounds=(0,None), initialize=0) m.x197 = Var(within=Reals, bounds=(0,None), initialize=0) m.x198 = Var(within=Reals, bounds=(0,None), initialize=0) m.x199 = Var(within=Reals, bounds=(0,None), initialize=0) m.x200 = Var(within=Reals, bounds=(0,None), initialize=0) m.x201 = Var(within=Reals, bounds=(0,None), initialize=0) m.x202 = Var(within=Reals, bounds=(0,None), initialize=0) m.x203 = Var(within=Reals, bounds=(0,None), initialize=0) m.x204 = Var(within=Reals, bounds=(0,None), initialize=0) m.x205 = Var(within=Reals, bounds=(0,None), initialize=0) m.x206 = Var(within=Reals, bounds=(0,None), initialize=0) m.x207 = Var(within=Reals, bounds=(0,None), initialize=0) m.x208 = Var(within=Reals, bounds=(0,None), initialize=0) m.x209 = Var(within=Reals, bounds=(0,None), initialize=0) m.x210 = Var(within=Reals, bounds=(0,None), initialize=0) m.x211 = Var(within=Reals, bounds=(0,None), initialize=0) m.x212 = Var(within=Reals, bounds=(0,None), initialize=0) m.x213 = Var(within=Reals, bounds=(0,None), initialize=0) m.x214 = Var(within=Reals, bounds=(0,None), initialize=0) m.x215 = Var(within=Reals, bounds=(0,None), initialize=0) m.x216 = Var(within=Reals, bounds=(0,None), initialize=0) m.x217 = Var(within=Reals, bounds=(0,None), initialize=0) m.x218 = Var(within=Reals, bounds=(0,None), initialize=0) m.x219 = Var(within=Reals, bounds=(0,None), initialize=0) m.x220 = Var(within=Reals, bounds=(0,None), initialize=0) m.x221 = Var(within=Reals, bounds=(0,None), initialize=0) m.x222 = Var(within=Reals, bounds=(0,None), initialize=0) m.x223 = Var(within=Reals, bounds=(0,None), initialize=0) m.x224 = Var(within=Reals, bounds=(0,None), initialize=0) m.x225 = Var(within=Reals, bounds=(0,None), initialize=0) m.x226 = Var(within=Reals, bounds=(0,None), initialize=0) m.x227 = Var(within=Reals, bounds=(0,None), initialize=0) m.x228 = Var(within=Reals, bounds=(0,None), initialize=0) m.x229 = Var(within=Reals, bounds=(0,None), initialize=0) m.x230 = Var(within=Reals, bounds=(0,None), initialize=0) m.x231 = Var(within=Reals, bounds=(0,None), initialize=0) m.x232 = Var(within=Reals, bounds=(0,None), initialize=0) m.x233 = Var(within=Reals, bounds=(0,None), initialize=0) m.x234 = Var(within=Reals, bounds=(0,None), initialize=0) m.x235 = Var(within=Reals, bounds=(0,None), initialize=0) m.x236 = Var(within=Reals, bounds=(0,None), initialize=0) m.x237 = Var(within=Reals, bounds=(0,None), initialize=0) m.x238 = Var(within=Reals, bounds=(0,None), initialize=0) m.x239 = Var(within=Reals, bounds=(0,None), initialize=0) m.x240 = Var(within=Reals, bounds=(0,None), initialize=0) m.x241 = Var(within=Reals, bounds=(0,None), initialize=0) m.x242 = Var(within=Reals, bounds=(0,None), initialize=0) m.b243 = Var(within=Binary, bounds=(0,1), initialize=0) m.b244 = Var(within=Binary, bounds=(0,1), initialize=0) m.b245 = Var(within=Binary, bounds=(0,1), initialize=0) m.b246 = Var(within=Binary, bounds=(0,1), initialize=0) m.b247 = Var(within=Binary, bounds=(0,1), initialize=0) m.b248 = Var(within=Binary, bounds=(0,1), initialize=0) m.b249 = Var(within=Binary, bounds=(0,1), initialize=0) m.b250 = Var(within=Binary, bounds=(0,1), initialize=0) m.b251 = Var(within=Binary, bounds=(0,1), initialize=0) m.b252 = Var(within=Binary, bounds=(0,1), initialize=0) m.b253 = Var(within=Binary, bounds=(0,1), initialize=0) m.b254 = Var(within=Binary, bounds=(0,1), initialize=0) m.b255 = Var(within=Binary, bounds=(0,1), initialize=0) m.b256 = Var(within=Binary, bounds=(0,1), initialize=0) m.b257 = Var(within=Binary, bounds=(0,1), initialize=0) m.b258 = Var(within=Binary, bounds=(0,1), initialize=0) m.b259 = Var(within=Binary, bounds=(0,1), initialize=0) m.b260 = Var(within=Binary, bounds=(0,1), initialize=0) m.b261 = Var(within=Binary, bounds=(0,1), initialize=0) m.b262 = Var(within=Binary, bounds=(0,1), initialize=0) m.b263 = Var(within=Binary, bounds=(0,1), initialize=0) m.b264 = Var(within=Binary, bounds=(0,1), initialize=0) m.b265 = Var(within=Binary, bounds=(0,1), initialize=0) m.b266 = Var(within=Binary, bounds=(0,1), initialize=0) m.b267 = Var(within=Binary, bounds=(0,1), initialize=0) m.b268 = Var(within=Binary, bounds=(0,1), initialize=0) m.b269 = Var(within=Binary, bounds=(0,1), initialize=0) m.b270 = Var(within=Binary, bounds=(0,1), initialize=0) m.b271 = Var(within=Binary, bounds=(0,1), initialize=0) m.b272 = Var(within=Binary, bounds=(0,1), initialize=0) m.b273 = Var(within=Binary, bounds=(0,1), initialize=0) m.b274 = Var(within=Binary, bounds=(0,1), initialize=0) m.b275 = Var(within=Binary, bounds=(0,1), initialize=0) m.b276 = Var(within=Binary, bounds=(0,1), initialize=0) m.b277 = Var(within=Binary, bounds=(0,1), initialize=0) m.b278 = Var(within=Binary, bounds=(0,1), initialize=0) m.b279 = Var(within=Binary, bounds=(0,1), initialize=0) m.b280 = Var(within=Binary, bounds=(0,1), initialize=0) m.b281 = Var(within=Binary, bounds=(0,1), initialize=0) m.b282 = Var(within=Binary, bounds=(0,1), initialize=0) m.b283 = Var(within=Binary, bounds=(0,1), initialize=0) m.b284 = Var(within=Binary, bounds=(0,1), initialize=0) m.b285 = Var(within=Binary, bounds=(0,1), initialize=0) m.b286 = Var(within=Binary, bounds=(0,1), initialize=0) m.b287 = Var(within=Binary, bounds=(0,1), initialize=0) m.b288 = Var(within=Binary, bounds=(0,1), initialize=0) m.b289 = Var(within=Binary, bounds=(0,1), initialize=0) m.b290 = Var(within=Binary, bounds=(0,1), initialize=0) m.b291 = Var(within=Binary, bounds=(0,1), initialize=0) m.b292 = Var(within=Binary, bounds=(0,1), initialize=0) m.b293 = Var(within=Binary, bounds=(0,1), initialize=0) m.b294 = Var(within=Binary, bounds=(0,1), initialize=0) m.b295 = Var(within=Binary, bounds=(0,1), initialize=0) m.b296 = Var(within=Binary, bounds=(0,1), initialize=0) m.b297 = Var(within=Binary, bounds=(0,1), initialize=0) m.b298 = Var(within=Binary, bounds=(0,1), initialize=0) m.b299 = Var(within=Binary, bounds=(0,1), initialize=0) m.b300 = Var(within=Binary, bounds=(0,1), initialize=0) m.b301 = Var(within=Binary, bounds=(0,1), initialize=0) m.b302 = Var(within=Binary, bounds=(0,1), initialize=0) m.obj = Objective(sense=maximize, expr= -m.x31 - m.x32 + 5 * m.x43 + 10 * m.x44 - 2 * m.x53 - m.x54 + 500 * m.x79 + 600 * m.x80 + 350 * m.x81 + 400 * m.x82 - 10 * m.x87 - 5 * m.x88 - 5 * m.x89 - 5 * m.x90 + 80 * m.x103 + 130 * m.x104 + 110 * m.x105 + 120 * m.x106 + 110 * m.x107 + 130 * m.x108 + 80 * m.x109 + 90 * m.x110 - 5 * m.b273 - 4 * m.b274 - 8 * m.b275 - 7 * m.b276 - 6 * m.b277 - 9 * m.b278 - 10 * m.b279 - 9 * m.b280 - 6 * m.b281 - 10 * m.b282 - 7 * m.b283 - 7 * m.b284 - 4 * m.b285 - 3 * m.b286 - 5 * m.b287 - 6 * m.b288 - 2 * m.b289 - 5 * m.b290 - 4 * m.b291 - 7 * m.b292 - 3 * m.b293 - 9 * m.b294 - 7 * m.b295 - 2 * m.b296 - 3 * m.b297 - m.b298 - 2 * m.b299 - 6 * m.b300 - 4 * m.b301 - 8 * m.b302) m.e1 = Constraint(expr= m.x31 - m.x33 - m.x35 == 0) m.e2 = Constraint(expr= m.x32 - m.x34 - m.x36 == 0) m.e3 = Constraint(expr= -m.x37 - m.x39 + m.x41 == 0) m.e4 = Constraint(expr= -m.x38 - m.x40 + m.x42 == 0) m.e5 = Constraint(expr= m.x41 - m.x43 - m.x45 == 0) m.e6 = Constraint(expr= m.x42 - m.x44 - m.x46 == 0) m.e7 = Constraint(expr= m.x45 - m.x47 - m.x49 - m.x51 == 0) m.e8 = Constraint(expr= m.x46 - m.x48 - m.x50 - m.x52 == 0) m.e9 = Constraint(expr= m.x55 - m.x61 - m.x63 == 0) m.e10 = Constraint(expr= m.x56 - m.x62 - m.x64 == 0) m.e11 = Constraint(expr= m.x59 - m.x65 - m.x67 - m.x69 == 0) m.e12 = Constraint(expr= m.x60 - m.x66 - m.x68 - m.x70 == 0) m.e13 = Constraint(expr= m.x75 - m.x83 - m.x85 == 0) m.e14 = Constraint(expr= m.x76 - m.x84 - m.x86 == 0) m.e15 = Constraint(expr= -m.x77 - m.x89 + m.x91 == 0) m.e16 = Constraint(expr= -m.x78 - m.x90 + m.x92 == 0) m.e17 = Constraint(expr= m.x79 - m.x93 - m.x95 == 0) m.e18 = Constraint(expr= m.x80 - m.x94 - m.x96 == 0) m.e19 = Constraint(expr= m.x81 - m.x97 - m.x99 - m.x101 == 0) m.e20 = Constraint(expr= m.x82 - m.x98 - m.x100 - m.x102 == 0) m.e21 = Constraint(expr= (m.x119 / (1e-06 + m.b243) - log(m.x111 / (1e-06 + m.b243) + 1)) * (1e-06 + m.b243) <= 0) m.e22 = Constraint(expr= (m.x120 / (1e-06 + m.b244) - log(m.x112 / (1e-06 + m.b244) + 1)) * (1e-06 + m.b244) <= 0) m.e23 = Constraint(expr= m.x113 == 0) m.e24 = Constraint(expr= m.x114 == 0) m.e25 = Constraint(expr= m.x121 == 0) m.e26 = Constraint(expr= m.x122 == 0) m.e27 = Constraint(expr= m.x33 - m.x111 - m.x113 == 0) m.e28 = Constraint(expr= m.x34 - m.x112 - m.x114 == 0) m.e29 = Constraint(expr= m.x37 - m.x119 - m.x121 == 0) m.e30 = Constraint(expr= m.x38 - m.x120 - m.x122 == 0) m.e31 = Constraint(expr= m.x111 - 40 * m.b243 <= 0) m.e32 = Constraint(expr= m.x112 - 40 * m.b244 <= 0) m.e33 = Constraint(expr= m.x113 + 40 * m.b243 <= 40) m.e34 = Constraint(expr= m.x114 + 40 * m.b244 <= 40) m.e35 = Constraint(expr= m.x119 - 3.71357206670431 * m.b243 <= 0) m.e36 = Constraint(expr= m.x120 - 3.71357206670431 * m.b244 <= 0) m.e37 = Constraint(expr= m.x121 + 3.71357206670431 * m.b243 <= 3.71357206670431) m.e38 = Constraint(expr= m.x122 + 3.71357206670431 * m.b244 <= 3.71357206670431) m.e39 = Constraint(expr= (m.x123 / (1e-06 + m.b245) - 1.2 * log(m.x115 / (1e-06 + m.b245) + 1)) * (1e-06 + m.b245) <= 0) m.e40 = Constraint(expr= (m.x124 / (1e-06 + m.b246) - 1.2 * log(m.x116 / (1e-06 + m.b246) + 1)) * (1e-06 + m.b246) <= 0) m.e41 = Constraint(expr= m.x117 == 0) m.e42 = Constraint(expr= m.x118 == 0) m.e43 = Constraint(expr= m.x125 == 0) m.e44 = Constraint(expr= m.x126 == 0) m.e45 = Constraint(expr= m.x35 - m.x115 - m.x117 == 0) m.e46 = Constraint(expr= m.x36 - m.x116 - m.x118 == 0) m.e47 = Constraint(expr= m.x39 - m.x123 - m.x125 == 0) m.e48 = Constraint(expr= m.x40 - m.x124 - m.x126 == 0) m.e49 = Constraint(expr= m.x115 - 40 * m.b245 <= 0) m.e50 = Constraint(expr= m.x116 - 40 * m.b246 <= 0) m.e51 = Constraint(expr= m.x117 + 40 * m.b245 <= 40) m.e52 = Constraint(expr= m.x118 + 40 * m.b246 <= 40) m.e53 = Constraint(expr= m.x123 - 4.45628648004517 * m.b245 <= 0) m.e54 = Constraint(expr= m.x124 - 4.45628648004517 * m.b246 <= 0) m.e55 = Constraint(expr= m.x125 + 4.45628648004517 * m.b245 <= 4.45628648004517) m.e56 = Constraint(expr= m.x126 + 4.45628648004517 * m.b246 <= 4.45628648004517) m.e57 = Constraint(expr= -0.75 * m.x127 + m.x143 == 0) m.e58 = Constraint(expr= -0.75 * m.x128 + m.x144 == 0) m.e59 = Constraint(expr= m.x129 == 0) m.e60 = Constraint(expr= m.x130 == 0) m.e61 = Constraint(expr= m.x145 == 0) m.e62 = Constraint(expr= m.x146 == 0) m.e63 = Constraint(expr= m.x47 - m.x127 - m.x129 == 0) m.e64 = Constraint(expr= m.x48 - m.x128 - m.x130 == 0) m.e65 = Constraint(expr= m.x55 - m.x143 - m.x145 == 0) m.e66 = Constraint(expr= m.x56 - m.x144 - m.x146 == 0) m.e67 = Constraint(expr= m.x127 - 4.45628648004517 * m.b247 <= 0) m.e68 = Constraint(expr= m.x128 - 4.45628648004517 * m.b248 <= 0) m.e69 = Constraint(expr= m.x129 + 4.45628648004517 * m.b247 <= 4.45628648004517) m.e70 = Constraint(expr= m.x130 + 4.45628648004517 * m.b248 <= 4.45628648004517) m.e71 = Constraint(expr= m.x143 - 3.34221486003388 * m.b247 <= 0) m.e72 = Constraint(expr= m.x144 - 3.34221486003388 * m.b248 <= 0) m.e73 = Constraint(expr= m.x145 + 3.34221486003388 * m.b247 <= 3.34221486003388) m.e74 = Constraint(expr= m.x146 + 3.34221486003388 * m.b248 <= 3.34221486003388) m.e75 = Constraint(expr= (m.x147 / (1e-06 + m.b249) - 1.5 * log(m.x131 / (1e-06 + m.b249) + 1)) * (1e-06 + m.b249) <= 0) m.e76 = Constraint(expr= (m.x148 / (1e-06 + m.b250) - 1.5 * log(m.x132 / (1e-06 + m.b250) + 1)) * (1e-06 + m.b250) <= 0) m.e77 = Constraint(expr= m.x133 == 0) m.e78 = Constraint(expr= m.x134 == 0) m.e79 = Constraint(expr= m.x151 == 0) m.e80 = Constraint(expr= m.x152 == 0) m.e81 = Constraint(expr= m.x49 - m.x131 - m.x133 == 0) m.e82 = Constraint(expr= m.x50 - m.x132 - m.x134 == 0) m.e83 = Constraint(expr= m.x57 - m.x147 - m.x151 == 0) m.e84 = Constraint(expr= m.x58 - m.x148 - m.x152 == 0) m.e85 = Constraint(expr= m.x131 - 4.45628648004517 * m.b249 <= 0) m.e86 = Constraint(expr= m.x132 - 4.45628648004517 * m.b250 <= 0) m.e87 = Constraint(expr= m.x133 + 4.45628648004517 * m.b249 <= 4.45628648004517) m.e88 = Constraint(expr= m.x134 + 4.45628648004517 * m.b250 <= 4.45628648004517) m.e89 = Constraint(expr= m.x147 - 2.54515263975353 * m.b249 <= 0) m.e90 = Constraint(expr= m.x148 - 2.54515263975353 * m.b250 <= 0) m.e91 = Constraint(expr= m.x151 + 2.54515263975353 * m.b249 <= 2.54515263975353) m.e92 = Constraint(expr= m.x152 + 2.54515263975353 * m.b250 <= 2.54515263975353) m.e93 = Constraint(expr= -m.x135 + m.x155 == 0) m.e94 = Constraint(expr= -m.x136 + m.x156 == 0) m.e95 = Constraint(expr= -0.5 * m.x139 + m.x155 == 0) m.e96 = Constraint(expr= -0.5 * m.x140 + m.x156 == 0) m.e97 = Constraint(expr= m.x137 == 0) m.e98 = Constraint(expr= m.x138 == 0) m.e99 = Constraint(expr= m.x141 == 0) m.e100 = Constraint(expr= m.x142 == 0) m.e101 = Constraint(expr= m.x157 == 0) m.e102 = Constraint(expr= m.x158 == 0) m.e103 = Constraint(expr= m.x51 - m.x135 - m.x137 == 0) m.e104 = Constraint(expr= m.x52 - m.x136 - m.x138 == 0) m.e105 = Constraint(expr= m.x53 - m.x139 - m.x141 == 0) m.e106 = Constraint(expr= m.x54 - m.x140 - m.x142 == 0) m.e107 = Constraint(expr= m.x59 - m.x155 - m.x157 == 0) m.e108 = Constraint(expr= m.x60 - m.x156 - m.x158 == 0) m.e109 = Constraint(expr= m.x135 - 4.45628648004517 * m.b251 <= 0) m.e110 = Constraint(expr= m.x136 - 4.45628648004517 * m.b252 <= 0) m.e111 = Constraint(expr= m.x137 + 4.45628648004517 * m.b251 <= 4.45628648004517) m.e112 = Constraint(expr= m.x138 + 4.45628648004517 * m.b252 <= 4.45628648004517) m.e113 = Constraint(expr= m.x139 - 30 * m.b251 <= 0) m.e114 = Constraint(expr= m.x140 - 30 * m.b252 <= 0) m.e115 = Constraint(expr= m.x141 + 30 * m.b251 <= 30) m.e116 = Constraint(expr= m.x142 + 30 * m.b252 <= 30) m.e117 = Constraint(expr= m.x155 - 15 * m.b251 <= 0) m.e118 = Constraint(expr= m.x156 - 15 * m.b252 <= 0) m.e119 = Constraint(expr= m.x157 + 15 * m.b251 <= 15) m.e120 = Constraint(expr= m.x158 + 15 * m.b252 <= 15) m.e121 = Constraint(expr= (m.x179 / (1e-06 + m.b253) - 1.25 * log(m.x159 / ( 1e-06 + m.b253) + 1)) * (1e-06 + m.b253) <= 0) m.e122 = Constraint(expr= (m.x180 / (1e-06 + m.b254) - 1.25 * log(m.x160 / ( 1e-06 + m.b254) + 1)) * (1e-06 + m.b254) <= 0) m.e123 = Constraint(expr= m.x161 == 0) m.e124 = Constraint(expr= m.x162 == 0) m.e125 = Constraint(expr= m.x183 == 0) m.e126 = Constraint(expr= m.x184 == 0) m.e127 = Constraint(expr= m.x61 - m.x159 - m.x161 == 0) m.e128 = Constraint(expr= m.x62 - m.x160 - m.x162 == 0) m.e129 = Constraint(expr= m.x71 - m.x179 - m.x183 == 0) m.e130 = Constraint(expr= m.x72 - m.x180 - m.x184 == 0) m.e131 = Constraint(expr= m.x159 - 3.34221486003388 * m.b253 <= 0) m.e132 = Constraint(expr= m.x160 - 3.34221486003388 * m.b254 <= 0) m.e133 = Constraint(expr= m.x161 + 3.34221486003388 * m.b253 <= 3.34221486003388) m.e134 = Constraint(expr= m.x162 + 3.34221486003388 * m.b254 <= 3.34221486003388) m.e135 = Constraint(expr= m.x179 - 1.83548069293539 * m.b253 <= 0) m.e136 = Constraint(expr= m.x180 - 1.83548069293539 * m.b254 <= 0) m.e137 = Constraint(expr= m.x183 + 1.83548069293539 * m.b253 <= 1.83548069293539) m.e138 = Constraint(expr= m.x184 + 1.83548069293539 * m.b254 <= 1.83548069293539) m.e139 = Constraint(expr= (m.x187 / (1e-06 + m.b255) - 0.9 * log(m.x163 / ( 1e-06 + m.b255) + 1)) * (1e-06 + m.b255) <= 0) m.e140 = Constraint(expr= (m.x188 / (1e-06 + m.b256) - 0.9 * log(m.x164 / ( 1e-06 + m.b256) + 1)) * (1e-06 + m.b256) <= 0) m.e141 = Constraint(expr= m.x165 == 0) m.e142 = Constraint(expr= m.x166 == 0) m.e143 = Constraint(expr= m.x191 == 0) m.e144 = Constraint(expr= m.x192 == 0) m.e145 = Constraint(expr= m.x63 - m.x163 - m.x165 == 0) m.e146 = Constraint(expr= m.x64 - m.x164 - m.x166 == 0) m.e147 = Constraint(expr= m.x73 - m.x187 - m.x191 == 0) m.e148 = Constraint(expr= m.x74 - m.x188 - m.x192 == 0) m.e149 = Constraint(expr= m.x163 - 3.34221486003388 * m.b255 <= 0) m.e150 = Constraint(expr= m.x164 - 3.34221486003388 * m.b256 <= 0) m.e151 = Constraint(expr= m.x165 + 3.34221486003388 * m.b255 <= 3.34221486003388) m.e152 = Constraint(expr= m.x166 + 3.34221486003388 * m.b256 <= 3.34221486003388) m.e153 = Constraint(expr= m.x187 - 1.32154609891348 * m.b255 <= 0) m.e154 = Constraint(expr= m.x188 - 1.32154609891348 * m.b256 <= 0) m.e155 = Constraint(expr= m.x191 + 1.32154609891348 * m.b255 <= 1.32154609891348) m.e156 = Constraint(expr= m.x192 + 1.32154609891348 * m.b256 <= 1.32154609891348) m.e157 = Constraint(expr= (m.x195 / (1e-06 + m.b257) - log(m.x149 / (1e-06 + m.b257) + 1)) * (1e-06 + m.b257) <= 0) m.e158 = Constraint(expr= (m.x196 / (1e-06 + m.b258) - log(m.x150 / (1e-06 + m.b258) + 1)) * (1e-06 + m.b258) <= 0) m.e159 = Constraint(expr= m.x153 == 0) m.e160 = Constraint(expr= m.x154 == 0) m.e161 = Constraint(expr= m.x197 == 0) m.e162 = Constraint(expr= m.x198 == 0) m.e163 = Constraint(expr= m.x57 - m.x149 - m.x153 == 0) m.e164 = Constraint(expr= m.x58 - m.x150 - m.x154 == 0) m.e165 = Constraint(expr= m.x75 - m.x195 - m.x197 == 0) m.e166 = Constraint(expr= m.x76 - m.x196 - m.x198 == 0) m.e167 = Constraint(expr= m.x149 - 2.54515263975353 * m.b257 <= 0) m.e168 = Constraint(expr= m.x150 - 2.54515263975353 * m.b258 <= 0) m.e169 = Constraint(expr= m.x153 + 2.54515263975353 * m.b257 <= 2.54515263975353) m.e170 = Constraint(expr= m.x154 + 2.54515263975353 * m.b258 <= 2.54515263975353) m.e171 = Constraint(expr= m.x195 - 1.26558121681553 * m.b257 <= 0) m.e172 = Constraint(expr= m.x196 - 1.26558121681553 * m.b258 <= 0) m.e173 = Constraint(expr= m.x197 + 1.26558121681553 * m.b257 <= 1.26558121681553) m.e174 = Constraint(expr= m.x198 + 1.26558121681553 * m.b258 <= 1.26558121681553) m.e175 = Constraint(expr= -0.9 * m.x167 + m.x199 == 0) m.e176 = Constraint(expr= -0.9 * m.x168 + m.x200 == 0) m.e177 = Constraint(expr= m.x169 == 0) m.e178 = Constraint(expr= m.x170 == 0) m.e179 = Constraint(expr= m.x201 == 0) m.e180 = Constraint(expr= m.x202 == 0) m.e181 = Constraint(expr= m.x65 - m.x167 - m.x169 == 0) m.e182 = Constraint(expr= m.x66 - m.x168 - m.x170 == 0) m.e183 = Constraint(expr= m.x77 - m.x199 - m.x201 == 0) m.e184 = Constraint(expr= m.x78 - m.x200 - m.x202 == 0) m.e185 = Constraint(expr= m.x167 - 15 * m.b259 <= 0) m.e186 = Constraint(expr= m.x168 - 15 * m.b260 <= 0) m.e187 = Constraint(expr= m.x169 + 15 * m.b259 <= 15) m.e188 = Constraint(expr= m.x170 + 15 * m.b260 <= 15) m.e189 = Constraint(expr= m.x199 - 13.5 * m.b259 <= 0) m.e190 = Constraint(expr= m.x200 - 13.5 * m.b260 <= 0) m.e191 = Constraint(expr= m.x201 + 13.5 * m.b259 <= 13.5) m.e192 = Constraint(expr= m.x202 + 13.5 * m.b260 <= 13.5) m.e193 = Constraint(expr= -0.6 * m.x171 + m.x203 == 0) m.e194 = Constraint(expr= -0.6 * m.x172 + m.x204 == 0) m.e195 = Constraint(expr= m.x173 == 0) m.e196 = Constraint(expr= m.x174 == 0) m.e197 = Constraint(expr= m.x205 == 0) m.e198 = Constraint(expr= m.x206 == 0) m.e199 = Constraint(expr= m.x67 - m.x171 - m.x173 == 0) m.e200 = Constraint(expr= m.x68 - m.x172 - m.x174 == 0) m.e201 = Constraint(expr= m.x79 - m.x203 - m.x205 == 0) m.e202 = Constraint(expr= m.x80 - m.x204 - m.x206 == 0) m.e203 = Constraint(expr= m.x171 - 15 * m.b261 <= 0) m.e204 = Constraint(expr= m.x172 - 15 * m.b262 <= 0) m.e205 = Constraint(expr= m.x173 + 15 * m.b261 <= 15) m.e206 = Constraint(expr= m.x174 + 15 * m.b262 <= 15) m.e207 = Constraint(expr= m.x203 - 9 * m.b261 <= 0) m.e208 = Constraint(expr= m.x204 - 9 * m.b262 <= 0) m.e209 = Constraint(expr= m.x205 + 9 * m.b261 <= 9) m.e210 = Constraint(expr= m.x206 + 9 * m.b262 <= 9) m.e211 = Constraint(expr= (m.x207 / (1e-06 + m.b263) - 1.1 * log(m.x175 / ( 1e-06 + m.b263) + 1)) * (1e-06 + m.b263) <= 0) m.e212 = Constraint(expr= (m.x208 / (1e-06 + m.b264) - 1.1 * log(m.x176 / ( 1e-06 + m.b264) + 1)) * (1e-06 + m.b264) <= 0) m.e213 = Constraint(expr= m.x177 == 0) m.e214 = Constraint(expr= m.x178 == 0) m.e215 = Constraint(expr= m.x209 == 0) m.e216 = Constraint(expr= m.x210 == 0) m.e217 = Constraint(expr= m.x69 - m.x175 - m.x177 == 0) m.e218 = Constraint(expr= m.x70 - m.x176 - m.x178 == 0) m.e219 = Constraint(expr= m.x81 - m.x207 - m.x209 == 0) m.e220 = Constraint(expr= m.x82 - m.x208 - m.x210 == 0) m.e221 = Constraint(expr= m.x175 - 15 * m.b263 <= 0) m.e222 = Constraint(expr= m.x176 - 15 * m.b264 <= 0) m.e223 = Constraint(expr= m.x177 + 15 * m.b263 <= 15) m.e224 = Constraint(expr= m.x178 + 15 * m.b264 <= 15) m.e225 = Constraint(expr= m.x207 - 3.04984759446376 * m.b263 <= 0) m.e226 = Constraint(expr= m.x208 - 3.04984759446376 * m.b264 <= 0) m.e227 = Constraint(expr= m.x209 + 3.04984759446376 * m.b263 <= 3.04984759446376) m.e228 = Constraint(expr= m.x210 + 3.04984759446376 * m.b264 <= 3.04984759446376) m.e229 = Constraint(expr= -0.9 * m.x181 + m.x227 == 0) m.e230 = Constraint(expr= -0.9 * m.x182 + m.x228 == 0) m.e231 = Constraint(expr= -m.x219 + m.x227 == 0) m.e232 = Constraint(expr= -m.x220 + m.x228 == 0) m.e233 = Constraint(expr= m.x185 == 0) m.e234 = Constraint(expr= m.x186 == 0) m.e235 = Constraint(expr= m.x221 == 0) m.e236 = Constraint(expr= m.x222 == 0) m.e237 = Constraint(expr= m.x229 == 0) m.e238 = Constraint(expr= m.x230 == 0) m.e239 = Constraint(expr= m.x71 - m.x181 - m.x185 == 0) m.e240 = Constraint(expr= m.x72 - m.x182 - m.x186 == 0) m.e241 = Constraint(expr= m.x87 - m.x219 - m.x221 == 0) m.e242 = Constraint(expr= m.x88 - m.x220 - m.x222 == 0) m.e243 = Constraint(expr= m.x103 - m.x227 - m.x229 == 0) m.e244 = Constraint(expr= m.x104 - m.x228 - m.x230 == 0) m.e245 = Constraint(expr= m.x181 - 1.83548069293539 * m.b265 <= 0) m.e246 = Constraint(expr= m.x182 - 1.83548069293539 * m.b266 <= 0) m.e247 = Constraint(expr= m.x185 + 1.83548069293539 * m.b265 <= 1.83548069293539) m.e248 = Constraint(expr= m.x186 + 1.83548069293539 * m.b266 <= 1.83548069293539) m.e249 = Constraint(expr= m.x219 - 20 * m.b265 <= 0) m.e250 = Constraint(expr= m.x220 - 20 * m.b266 <= 0) m.e251 = Constraint(expr= m.x221 + 20 * m.b265 <= 20) m.e252 = Constraint(expr= m.x222 + 20 * m.b266 <= 20) m.e253 = Constraint(expr= m.x227 - 20 * m.b265 <= 0) m.e254 = Constraint(expr= m.x228 - 20 * m.b266 <= 0) m.e255 = Constraint(expr= m.x229 + 20 * m.b265 <= 20) m.e256 = Constraint(expr= m.x230 + 20 * m.b266 <= 20) m.e257 = Constraint(expr= (m.x231 / (1e-06 + m.b267) - log(m.x189 / (1e-06 + m.b267) + 1)) * (1e-06 + m.b267) <= 0) m.e258 = Constraint(expr= (m.x232 / (1e-06 + m.b268) - log(m.x190 / (1e-06 + m.b268) + 1)) * (1e-06 + m.b268) <= 0) m.e259 = Constraint(expr= m.x193 == 0) m.e260 = Constraint(expr= m.x194 == 0) m.e261 = Constraint(expr= m.x233 == 0) m.e262 = Constraint(expr= m.x234 == 0) m.e263 = Constraint(expr= m.x73 - m.x189 - m.x193 == 0) m.e264 = Constraint(expr= m.x74 - m.x190 - m.x194 == 0) m.e265 = Constraint(expr= m.x105 - m.x231 - m.x233 == 0) m.e266 = Constraint(expr= m.x106 - m.x232 - m.x234 == 0) m.e267 = Constraint(expr= m.x189 - 1.32154609891348 * m.b267 <= 0) m.e268 = Constraint(expr= m.x190 - 1.32154609891348 * m.b268 <= 0) m.e269 = Constraint(expr= m.x193 + 1.32154609891348 * m.b267 <= 1.32154609891348) m.e270 = Constraint(expr= m.x194 + 1.32154609891348 * m.b268 <= 1.32154609891348) m.e271 = Constraint(expr= m.x231 - 0.842233385663186 * m.b267 <= 0) m.e272 = Constraint(expr= m.x232 - 0.842233385663186 * m.b268 <= 0) m.e273 = Constraint(expr= m.x233 + 0.842233385663186 * m.b267 <= 0.842233385663186) m.e274 = Constraint(expr= m.x234 + 0.842233385663186 * m.b268 <= 0.842233385663186) m.e275 = Constraint(expr= (m.x235 / (1e-06 + m.b269) - 0.7 * log(m.x211 / ( 1e-06 + m.b269) + 1)) * (1e-06 + m.b269) <= 0) m.e276 = Constraint(expr= (m.x236 / (1e-06 + m.b270) - 0.7 * log(m.x212 / ( 1e-06 + m.b270) + 1)) * (1e-06 + m.b270) <= 0) m.e277 = Constraint(expr= m.x213 == 0) m.e278 = Constraint(expr= m.x214 == 0) m.e279 = Constraint(expr= m.x237 == 0) m.e280 = Constraint(expr= m.x238 == 0) m.e281 = Constraint(expr= m.x83 - m.x211 - m.x213 == 0) m.e282 = Constraint(expr= m.x84 - m.x212 - m.x214 == 0) m.e283 = Constraint(expr= m.x107 - m.x235 - m.x237 == 0) m.e284 = Constraint(expr= m.x108 - m.x236 - m.x238 == 0) m.e285 = Constraint(expr= m.x211 - 1.26558121681553 * m.b269 <= 0) m.e286 = Constraint(expr= m.x212 - 1.26558121681553 * m.b270 <= 0) m.e287 = Constraint(expr= m.x213 + 1.26558121681553 * m.b269 <= 1.26558121681553) m.e288 = Constraint(expr= m.x214 + 1.26558121681553 * m.b270 <= 1.26558121681553) m.e289 = Constraint(expr= m.x235 - 0.572481933717686 * m.b269 <= 0) m.e290 = Constraint(expr= m.x236 - 0.572481933717686 * m.b270 <= 0) m.e291 = Constraint(expr= m.x237 + 0.572481933717686 * m.b269 <= 0.572481933717686) m.e292 = Constraint(expr= m.x238 + 0.572481933717686 * m.b270 <= 0.572481933717686) m.e293 = Constraint(expr= (m.x239 / (1e-06 + m.b271) - 0.65 * log(m.x215 / ( 1e-06 + m.b271) + 1)) * (1e-06 + m.b271) <= 0) m.e294 = Constraint(expr= (m.x240 / (1e-06 + m.b272) - 0.65 * log(m.x216 / ( 1e-06 + m.b272) + 1)) * (1e-06 + m.b272) <= 0) m.e295 = Constraint(expr= (m.x239 / (1e-06 + m.b271) - 0.65 * log(m.x223 / ( 1e-06 + m.b271) + 1)) * (1e-06 + m.b271) <= 0) m.e296 = Constraint(expr= (m.x240 / (1e-06 + m.b272) - 0.65 * log(m.x224 / ( 1e-06 + m.b272) + 1)) * (1e-06 + m.b272) <= 0) m.e297 = Constraint(expr= m.x217 == 0) m.e298 = Constraint(expr= m.x218 == 0) m.e299 = Constraint(expr= m.x225 == 0) m.e300 = Constraint(expr= m.x226 == 0) m.e301 = Constraint(expr= m.x241 == 0) m.e302 = Constraint(expr= m.x242 == 0) m.e303 = Constraint(expr= m.x85 - m.x215 - m.x217 == 0) m.e304 = Constraint(expr= m.x86 - m.x216 - m.x218 == 0) m.e305 = Constraint(expr= m.x91 - m.x223 - m.x225 == 0) m.e306 = Constraint(expr= m.x92 - m.x224 - m.x226 == 0) m.e307 = Constraint(expr= m.x109 - m.x239 - m.x241 == 0) m.e308 = Constraint(expr= m.x110 - m.x240 - m.x242 == 0) m.e309 = Constraint(expr= m.x215 - 1.26558121681553 * m.b271 <= 0) m.e310 = Constraint(expr= m.x216 - 1.26558121681553 * m.b272 <= 0) m.e311 = Constraint(expr= m.x217 + 1.26558121681553 * m.b271 <= 1.26558121681553) m.e312 = Constraint(expr= m.x218 + 1.26558121681553 * m.b272 <= 1.26558121681553) m.e313 = Constraint(expr= m.x223 - 33.5 * m.b271 <= 0) m.e314 = Constraint(expr= m.x224 - 33.5 * m.b272 <= 0) m.e315 = Constraint(expr= m.x225 + 33.5 * m.b271 <= 33.5) m.e316 = Constraint(expr= m.x226 + 33.5 * m.b272 <= 33.5) m.e317 = Constraint(expr= m.x239 - 2.30162356062425 * m.b271 <= 0) m.e318 = Constraint(expr= m.x240 - 2.30162356062425 * m.b272 <= 0) m.e319 = Constraint(expr= m.x241 + 2.30162356062425 * m.b271 <= 2.30162356062425) m.e320 = Constraint(expr= m.x242 + 2.30162356062425 * m.b272 <= 2.30162356062425) m.e321 = Constraint(expr= m.x1 + 5 * m.b273 == 0) m.e322 = Constraint(expr= m.x2 + 4 * m.b274 == 0) m.e323 = Constraint(expr= m.x3 + 8 * m.b275 == 0) m.e324 = Constraint(expr= m.x4 + 7 * m.b276 == 0) m.e325 = Constraint(expr= m.x5 + 6 * m.b277 == 0) m.e326 = Constraint(expr= m.x6 + 9 * m.b278 == 0) m.e327 = Constraint(expr= m.x7 + 10 * m.b279 == 0) m.e328 = Constraint(expr= m.x8 + 9 * m.b280 == 0) m.e329 = Constraint(expr= m.x9 + 6 * m.b281 == 0) m.e330 = Constraint(expr= m.x10 + 10 * m.b282 == 0) m.e331 = Constraint(expr= m.x11 + 7 * m.b283 == 0) m.e332 = Constraint(expr= m.x12 + 7 * m.b284 == 0) m.e333 = Constraint(expr= m.x13 + 4 * m.b285 == 0) m.e334 = Constraint(expr= m.x14 + 3 * m.b286 == 0) m.e335 = Constraint(expr= m.x15 + 5 * m.b287 == 0) m.e336 = Constraint(expr= m.x16 + 6 * m.b288 == 0) m.e337 = Constraint(expr= m.x17 + 2 * m.b289 == 0) m.e338 = Constraint(expr= m.x18 + 5 * m.b290 == 0) m.e339 = Constraint(expr= m.x19 + 4 * m.b291 == 0) m.e340 = Constraint(expr= m.x20 + 7 * m.b292 == 0) m.e341 = Constraint(expr= m.x21 + 3 * m.b293 == 0) m.e342 = Constraint(expr= m.x22 + 9 * m.b294 == 0) m.e343 = Constraint(expr= m.x23 + 7 * m.b295 == 0) m.e344 = Constraint(expr= m.x24 + 2 * m.b296 == 0) m.e345 = Constraint(expr= m.x25 + 3 * m.b297 == 0) m.e346 = Constraint(expr= m.x26 + m.b298 == 0) m.e347 = Constraint(expr= m.x27 + 2 * m.b299 == 0) m.e348 = Constraint(expr= m.x28 + 6 * m.b300 == 0) m.e349 = Constraint(expr= m.x29 + 4 * m.b301 == 0) m.e350 = Constraint(expr= m.x30 + 8 * m.b302 == 0) m.e351 = Constraint(expr= m.b243 - m.b244 <= 0) m.e352 = Constraint(expr= m.b245 - m.b246 <= 0) m.e353 = Constraint(expr= m.b247 - m.b248 <= 0) m.e354 = Constraint(expr= m.b249 - m.b250 <= 0) m.e355 = Constraint(expr= m.b251 - m.b252 <= 0) m.e356 = Constraint(expr= m.b253 - m.b254 <= 0) m.e357 = Constraint(expr= m.b255 - m.b256 <= 0) m.e358 = Constraint(expr= m.b257 - m.b258 <= 0) m.e359 = Constraint(expr= m.b259 - m.b260 <= 0) m.e360 = Constraint(expr= m.b261 - m.b262 <= 0) m.e361 = Constraint(expr= m.b263 - m.b264 <= 0) m.e362 = Constraint(expr= m.b265 - m.b266 <= 0) m.e363 = Constraint(expr= m.b267 - m.b268 <= 0) m.e364 = Constraint(expr= m.b269 - m.b270 <= 0) m.e365 = Constraint(expr= m.b271 - m.b272 <= 0) m.e366 = Constraint(expr= m.b273 + m.b274 <= 1) m.e367 = Constraint(expr= m.b273 + m.b274 <= 1) m.e368 = Constraint(expr= m.b275 + m.b276 <= 1) m.e369 = Constraint(expr= m.b275 + m.b276 <= 1) m.e370 = Constraint(expr= m.b277 + m.b278 <= 1) m.e371 = Constraint(expr= m.b277 + m.b278 <= 1) m.e372 = Constraint(expr= m.b279 + m.b280 <= 1) m.e373 = Constraint(expr= m.b279 + m.b280 <= 1) m.e374 = Constraint(expr= m.b281 + m.b282 <= 1) m.e375 = Constraint(expr= m.b281 + m.b282 <= 1) m.e376 = Constraint(expr= m.b283 + m.b284 <= 1) m.e377 = Constraint(expr= m.b283 + m.b284 <= 1) m.e378 = Constraint(expr= m.b285 + m.b286 <= 1) m.e379 = Constraint(expr= m.b285 + m.b286 <= 1) m.e380 = Constraint(expr= m.b287 + m.b288 <= 1) m.e381 = Constraint(expr= m.b287 + m.b288 <= 1) m.e382 = Constraint(expr= m.b289 + m.b290 <= 1) m.e383 = Constraint(expr= m.b289 + m.b290 <= 1) m.e384 = Constraint(expr= m.b291 + m.b292 <= 1) m.e385 = Constraint(expr= m.b291 + m.b292 <= 1) m.e386 = Constraint(expr= m.b293 + m.b294 <= 1) m.e387 = Constraint(expr= m.b293 + m.b294 <= 1) m.e388 = Constraint(expr= m.b295 + m.b296 <= 1) m.e389 = Constraint(expr= m.b295 + m.b296 <= 1) m.e390 = Constraint(expr= m.b297 + m.b298 <= 1) m.e391 = Constraint(expr= m.b297 + m.b298 <= 1) m.e392 = Constraint(expr= m.b299 + m.b300 <= 1) m.e393 = Constraint(expr= m.b299 + m.b300 <= 1) m.e394 = Constraint(expr= m.b301 + m.b302 <= 1) m.e395 = Constraint(expr= m.b301 + m.b302 <= 1) m.e396 = Constraint(expr= m.b243 - m.b273 <= 0) m.e397 = Constraint(expr= -m.b243 + m.b244 - m.b274 <= 0) m.e398 = Constraint(expr= m.b245 - m.b275 <= 0) m.e399 = Constraint(expr= -m.b245 + m.b246 - m.b276 <= 0) m.e400 = Constraint(expr= m.b247 - m.b277 <= 0) m.e401 = Constraint(expr= -m.b247 + m.b248 - m.b278 <= 0) m.e402 = Constraint(expr= m.b249 - m.b279 <= 0) m.e403 = Constraint(expr= -m.b249 + m.b250 - m.b280 <= 0) m.e404 = Constraint(expr= m.b251 - m.b281 <= 0) m.e405 = Constraint(expr= -m.b251 + m.b252 - m.b282 <= 0) m.e406 = Constraint(expr= m.b253 - m.b283 <= 0) m.e407 = Constraint(expr= -m.b253 + m.b254 - m.b284 <= 0) m.e408 = Constraint(expr= m.b255 - m.b285 <= 0) m.e409 = Constraint(expr= -m.b255 + m.b256 - m.b286 <= 0) m.e410 = Constraint(expr= m.b257 - m.b287 <= 0) m.e411 = Constraint(expr= -m.b257 + m.b258 - m.b288 <= 0) m.e412 = Constraint(expr= m.b259 - m.b289 <= 0) m.e413 = Constraint(expr= -m.b259 + m.b260 - m.b290 <= 0) m.e414 = Constraint(expr= m.b261 - m.b291 <= 0) m.e415 = Constraint(expr= -m.b261 + m.b262 - m.b292 <= 0) m.e416 = Constraint(expr= m.b263 - m.b293 <= 0) m.e417 = Constraint(expr= -m.b263 + m.b264 - m.b294 <= 0) m.e418 = Constraint(expr= m.b265 - m.b295 <= 0) m.e419 = Constraint(expr= -m.b265 + m.b266 - m.b296 <= 0) m.e420 = Constraint(expr= m.b267 - m.b297 <= 0) m.e421 = Constraint(expr= -m.b267 + m.b268 - m.b298 <= 0) m.e422 = Constraint(expr= m.b269 - m.b299 <= 0) m.e423 = Constraint(expr= -m.b269 + m.b270 - m.b300 <= 0) m.e424 = Constraint(expr= m.b271 - m.b301 <= 0) m.e425 = Constraint(expr= -m.b271 + m.b272 - m.b302 <= 0) m.e426 = Constraint(expr= m.b243 + m.b245 == 1) m.e427 = Constraint(expr= m.b244 + m.b246 == 1) m.e428 = Constraint(expr= -m.b247 + m.b253 + m.b255 >= 0) m.e429 = Constraint(expr= -m.b248 + m.b254 + m.b256 >= 0) m.e430 = Constraint(expr= -m.b253 + m.b265 >= 0) m.e431 = Constraint(expr= -m.b254 + m.b266 >= 0) m.e432 = Constraint(expr= -m.b255 + m.b267 >= 0) m.e433 = Constraint(expr= -m.b256 + m.b268 >= 0) m.e434 = Constraint(expr= -m.b249 + m.b257 >= 0) m.e435 = Constraint(expr= -m.b250 + m.b258 >= 0) m.e436 = Constraint(expr= -m.b257 + m.b269 + m.b271 >= 0) m.e437 = Constraint(expr= -m.b258 + m.b270 + m.b272 >= 0) m.e438 = Constraint(expr= -m.b251 + m.b259 + m.b261 + m.b263 >= 0) m.e439 = Constraint(expr= -m.b252 + m.b260 + m.b262 + m.b264 >= 0) m.e440 = Constraint(expr= -m.b259 + m.b271 >= 0) m.e441 = Constraint(expr= -m.b260 + m.b272 >= 0) m.e442 = Constraint(expr= m.b243 + m.b245 - m.b247 >= 0) m.e443 = Constraint(expr= m.b244 + m.b246 - m.b248 >= 0) m.e444 = Constraint(expr= m.b243 + m.b245 - m.b249 >= 0) m.e445 = Constraint(expr= m.b244 + m.b246 - m.b250 >= 0) m.e446 = Constraint(expr= m.b243 + m.b245 - m.b251 >= 0) m.e447 = Constraint(expr= m.b244 + m.b246 - m.b252 >= 0) m.e448 = Constraint(expr= m.b247 - m.b253 >= 0) m.e449 = Constraint(expr= m.b248 - m.b254 >= 0) m.e450 = Constraint(expr= m.b247 - m.b255 >= 0) m.e451 = Constraint(expr= m.b248 - m.b256 >= 0) m.e452 = Constraint(expr= m.b249 - m.b257 >= 0) m.e453 = Constraint(expr= m.b250 - m.b258 >= 0) m.e454 = Constraint(expr= m.b251 - m.b259 >= 0) m.e455 = Constraint(expr= m.b252 - m.b260 >= 0) m.e456 = Constraint(expr= m.b251 - m.b261 >= 0) m.e457 = Constraint(expr= m.b252 - m.b262 >= 0) m.e458 = Constraint(expr= m.b251 - m.b263 >= 0) m.e459 = Constraint(expr= m.b252 - m.b264 >= 0) m.e460 = Constraint(expr= m.b253 - m.b265 >= 0) m.e461 = Constraint(expr= m.b254 - m.b266 >= 0) m.e462 = Constraint(expr= m.b255 - m.b267 >= 0) m.e463 = Constraint(expr= m.b256 - m.b268 >= 0) m.e464 = Constraint(expr= m.b257 - m.b269 >= 0) m.e465 = Constraint(expr= m.b258 - m.b270 >= 0) m.e466 = Constraint(expr= m.b257 - m.b271 >= 0) m.e467 = Constraint(expr= m.b258 - m.b272 >= 0)