<html><head><meta name="color-scheme" content="light dark"></head><body><pre style="word-wrap: break-word; white-space: pre-wrap;"># MINLP written by GAMS Convert at 03/03/22 14:44:51
#
# Equation counts
#     Total        E        G        L        N        X        C        B
#       708      256       40      412        0        0        0        0
#
# Variable counts
#                  x        b        i      s1s      s2s       sc       si
#     Total     cont   binary  integer     sos1     sos2    scont     sint
#       388      308       80        0        0        0        0        0
# FX      0
#
# Nonzero counts
#     Total    const       NL
#      1548     1476       72
#
# 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=(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=(0,40), initialize=0)
m.x42 = Var(within=Reals, bounds=(0,40), initialize=0)
m.x43 = Var(within=Reals, bounds=(0,40), initialize=0)
m.x44 = Var(within=Reals, bounds=(0,40), 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,None), initialize=0)
m.x54 = Var(within=Reals, bounds=(0,None), 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,30), initialize=0)
m.x86 = Var(within=Reals, bounds=(0,30), initialize=0)
m.x87 = Var(within=Reals, bounds=(0,30), initialize=0)
m.x88 = Var(within=Reals, bounds=(0,30), initialize=0)
m.x89 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x90 = Var(within=Reals, bounds=(0,None), 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.x243 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x244 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x245 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x246 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x247 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x248 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x249 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x250 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x251 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x252 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x253 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x254 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x255 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x256 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x257 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x258 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x259 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x260 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x261 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x262 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x263 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x264 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x265 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x266 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x267 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x268 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x269 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x270 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x271 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x272 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x273 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x274 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x275 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x276 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x277 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x278 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x279 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x280 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x281 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x282 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x283 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x284 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x285 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x286 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x287 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x288 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x289 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x290 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x291 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x292 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x293 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x294 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x295 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x296 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x297 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x298 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x299 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x300 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x301 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x302 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x303 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x304 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x305 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x306 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x307 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x308 = Var(within=Reals, bounds=(0,None), initialize=0)
m.b309 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b310 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b311 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b312 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b313 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b314 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b315 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b316 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b317 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b318 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b319 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b320 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b321 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b322 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b323 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b324 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b325 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b326 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b327 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b328 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b329 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b330 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b331 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b332 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b333 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b334 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b335 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b336 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b337 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b338 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b339 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b340 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b341 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b342 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b343 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b344 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b345 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b346 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b347 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b348 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b349 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b350 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b351 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b352 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b353 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b354 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b355 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b356 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b357 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b358 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b359 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b360 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b361 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b362 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b363 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b364 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b365 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b366 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b367 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b368 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b369 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b370 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b371 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b372 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b373 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b374 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b375 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b376 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b377 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b378 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b379 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b380 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b381 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b382 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b383 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b384 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b385 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b386 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b387 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b388 = Var(within=Binary, bounds=(0,1), initialize=0)

m.obj = Objective(sense=maximize, expr= -m.x41 - m.x42 - m.x43 - m.x44 + 5 *
    m.x65 + 10 * m.x66 + 5 * m.x67 + 10 * m.x68 - 2 * m.x85 - m.x86 - 2 * m.x87
    - m.x88 + 80 * m.x117 + 90 * m.x118 + 120 * m.x119 + 100 * m.x120 + 285 *
    m.x121 + 390 * m.x122 + 350 * m.x123 + 300 * m.x124 + 290 * m.x125 + 405 *
    m.x126 + 190 * m.x127 + 340 * m.x128 + 280 * m.x129 + 400 * m.x130 + 430 *
    m.x131 + 260 * m.x132 + 290 * m.x133 + 300 * m.x134 + 240 * m.x135 + 310 *
    m.x136 + 350 * m.x137 + 250 * m.x138 + 300 * m.x139 + 400 * m.x140 - 5 *
    m.b349 - 4 * m.b350 - 6 * m.b351 - 3 * m.b352 - 8 * m.b353 - 7 * m.b354 - 6
    * m.b355 - 5 * m.b356 - 6 * m.b357 - 9 * m.b358 - 4 * m.b359 - 3 * m.b360
    - 10 * m.b361 - 9 * m.b362 - 5 * m.b363 - 6 * m.b364 - 6 * m.b365 - 10 *
    m.b366 - 6 * m.b367 - 9 * m.b368 - 7 * m.b369 - 7 * m.b370 - 4 * m.b371 - 2
    * m.b372 - 4 * m.b373 - 3 * m.b374 - 2 * m.b375 - 8 * m.b376 - 5 * m.b377
    - 6 * m.b378 - 7 * m.b379 - 4 * m.b380 - 2 * m.b381 - 5 * m.b382 - 2 *
    m.b383 - 6 * m.b384 - 4 * m.b385 - 7 * m.b386 - 4 * m.b387 - 7 * m.b388)

m.e1 = Constraint(expr= m.x41 - m.x45 - m.x49 == 0)
m.e2 = Constraint(expr= m.x42 - m.x46 - m.x50 == 0)
m.e3 = Constraint(expr= m.x43 - m.x47 - m.x51 == 0)
m.e4 = Constraint(expr= m.x44 - m.x48 - m.x52 == 0)
m.e5 = Constraint(expr= -m.x53 - m.x57 + m.x61 == 0)
m.e6 = Constraint(expr= -m.x54 - m.x58 + m.x62 == 0)
m.e7 = Constraint(expr= -m.x55 - m.x59 + m.x63 == 0)
m.e8 = Constraint(expr= -m.x56 - m.x60 + m.x64 == 0)
m.e9 = Constraint(expr= m.x61 - m.x65 - m.x69 == 0)
m.e10 = Constraint(expr= m.x62 - m.x66 - m.x70 == 0)
m.e11 = Constraint(expr= m.x63 - m.x67 - m.x71 == 0)
m.e12 = Constraint(expr= m.x64 - m.x68 - m.x72 == 0)
m.e13 = Constraint(expr= m.x69 - m.x73 - m.x77 - m.x81 == 0)
m.e14 = Constraint(expr= m.x70 - m.x74 - m.x78 - m.x82 == 0)
m.e15 = Constraint(expr= m.x71 - m.x75 - m.x79 - m.x83 == 0)
m.e16 = Constraint(expr= m.x72 - m.x76 - m.x80 - m.x84 == 0)
m.e17 = Constraint(expr= m.x89 - m.x101 - m.x105 == 0)
m.e18 = Constraint(expr= m.x90 - m.x102 - m.x106 == 0)
m.e19 = Constraint(expr= m.x91 - m.x103 - m.x107 == 0)
m.e20 = Constraint(expr= m.x92 - m.x104 - m.x108 == 0)
m.e21 = Constraint(expr= m.x97 - m.x109 - m.x113 - m.x117 == 0)
m.e22 = Constraint(expr= m.x98 - m.x110 - m.x114 - m.x118 == 0)
m.e23 = Constraint(expr= m.x99 - m.x111 - m.x115 - m.x119 == 0)
m.e24 = Constraint(expr= m.x100 - m.x112 - m.x116 - m.x120 == 0)
m.e25 = Constraint(expr= (m.x157 / (1e-06 + m.b309) - log(m.x141 / (1e-06 +
    m.b309) + 1)) * (1e-06 + m.b309) &lt;= 0)
m.e26 = Constraint(expr= (m.x158 / (1e-06 + m.b310) - log(m.x142 / (1e-06 +
    m.b310) + 1)) * (1e-06 + m.b310) &lt;= 0)
m.e27 = Constraint(expr= (m.x159 / (1e-06 + m.b311) - log(m.x143 / (1e-06 +
    m.b311) + 1)) * (1e-06 + m.b311) &lt;= 0)
m.e28 = Constraint(expr= (m.x160 / (1e-06 + m.b312) - log(m.x144 / (1e-06 +
    m.b312) + 1)) * (1e-06 + m.b312) &lt;= 0)
m.e29 = Constraint(expr= m.x145 == 0)
m.e30 = Constraint(expr= m.x146 == 0)
m.e31 = Constraint(expr= m.x147 == 0)
m.e32 = Constraint(expr= m.x148 == 0)
m.e33 = Constraint(expr= m.x161 == 0)
m.e34 = Constraint(expr= m.x162 == 0)
m.e35 = Constraint(expr= m.x163 == 0)
m.e36 = Constraint(expr= m.x164 == 0)
m.e37 = Constraint(expr= m.x45 - m.x141 - m.x145 == 0)
m.e38 = Constraint(expr= m.x46 - m.x142 - m.x146 == 0)
m.e39 = Constraint(expr= m.x47 - m.x143 - m.x147 == 0)
m.e40 = Constraint(expr= m.x48 - m.x144 - m.x148 == 0)
m.e41 = Constraint(expr= m.x53 - m.x157 - m.x161 == 0)
m.e42 = Constraint(expr= m.x54 - m.x158 - m.x162 == 0)
m.e43 = Constraint(expr= m.x55 - m.x159 - m.x163 == 0)
m.e44 = Constraint(expr= m.x56 - m.x160 - m.x164 == 0)
m.e45 = Constraint(expr= m.x141 - 40 * m.b309 &lt;= 0)
m.e46 = Constraint(expr= m.x142 - 40 * m.b310 &lt;= 0)
m.e47 = Constraint(expr= m.x143 - 40 * m.b311 &lt;= 0)
m.e48 = Constraint(expr= m.x144 - 40 * m.b312 &lt;= 0)
m.e49 = Constraint(expr= m.x145 + 40 * m.b309 &lt;= 40)
m.e50 = Constraint(expr= m.x146 + 40 * m.b310 &lt;= 40)
m.e51 = Constraint(expr= m.x147 + 40 * m.b311 &lt;= 40)
m.e52 = Constraint(expr= m.x148 + 40 * m.b312 &lt;= 40)
m.e53 = Constraint(expr= m.x157 - 3.71357206670431 * m.b309 &lt;= 0)
m.e54 = Constraint(expr= m.x158 - 3.71357206670431 * m.b310 &lt;= 0)
m.e55 = Constraint(expr= m.x159 - 3.71357206670431 * m.b311 &lt;= 0)
m.e56 = Constraint(expr= m.x160 - 3.71357206670431 * m.b312 &lt;= 0)
m.e57 = Constraint(expr= m.x161 + 3.71357206670431 * m.b309
    &lt;= 3.71357206670431)
m.e58 = Constraint(expr= m.x162 + 3.71357206670431 * m.b310
    &lt;= 3.71357206670431)
m.e59 = Constraint(expr= m.x163 + 3.71357206670431 * m.b311
    &lt;= 3.71357206670431)
m.e60 = Constraint(expr= m.x164 + 3.71357206670431 * m.b312
    &lt;= 3.71357206670431)
m.e61 = Constraint(expr= (m.x165 / (1e-06 + m.b313) - 1.2 * log(m.x149 / (1e-06
    + m.b313) + 1)) * (1e-06 + m.b313) &lt;= 0)
m.e62 = Constraint(expr= (m.x166 / (1e-06 + m.b314) - 1.2 * log(m.x150 / (1e-06
    + m.b314) + 1)) * (1e-06 + m.b314) &lt;= 0)
m.e63 = Constraint(expr= (m.x167 / (1e-06 + m.b315) - 1.2 * log(m.x151 / (1e-06
    + m.b315) + 1)) * (1e-06 + m.b315) &lt;= 0)
m.e64 = Constraint(expr= (m.x168 / (1e-06 + m.b316) - 1.2 * log(m.x152 / (1e-06
    + m.b316) + 1)) * (1e-06 + m.b316) &lt;= 0)
m.e65 = Constraint(expr= m.x153 == 0)
m.e66 = Constraint(expr= m.x154 == 0)
m.e67 = Constraint(expr= m.x155 == 0)
m.e68 = Constraint(expr= m.x156 == 0)
m.e69 = Constraint(expr= m.x169 == 0)
m.e70 = Constraint(expr= m.x170 == 0)
m.e71 = Constraint(expr= m.x171 == 0)
m.e72 = Constraint(expr= m.x172 == 0)
m.e73 = Constraint(expr= m.x49 - m.x149 - m.x153 == 0)
m.e74 = Constraint(expr= m.x50 - m.x150 - m.x154 == 0)
m.e75 = Constraint(expr= m.x51 - m.x151 - m.x155 == 0)
m.e76 = Constraint(expr= m.x52 - m.x152 - m.x156 == 0)
m.e77 = Constraint(expr= m.x57 - m.x165 - m.x169 == 0)
m.e78 = Constraint(expr= m.x58 - m.x166 - m.x170 == 0)
m.e79 = Constraint(expr= m.x59 - m.x167 - m.x171 == 0)
m.e80 = Constraint(expr= m.x60 - m.x168 - m.x172 == 0)
m.e81 = Constraint(expr= m.x149 - 40 * m.b313 &lt;= 0)
m.e82 = Constraint(expr= m.x150 - 40 * m.b314 &lt;= 0)
m.e83 = Constraint(expr= m.x151 - 40 * m.b315 &lt;= 0)
m.e84 = Constraint(expr= m.x152 - 40 * m.b316 &lt;= 0)
m.e85 = Constraint(expr= m.x153 + 40 * m.b313 &lt;= 40)
m.e86 = Constraint(expr= m.x154 + 40 * m.b314 &lt;= 40)
m.e87 = Constraint(expr= m.x155 + 40 * m.b315 &lt;= 40)
m.e88 = Constraint(expr= m.x156 + 40 * m.b316 &lt;= 40)
m.e89 = Constraint(expr= m.x165 - 4.45628648004517 * m.b313 &lt;= 0)
m.e90 = Constraint(expr= m.x166 - 4.45628648004517 * m.b314 &lt;= 0)
m.e91 = Constraint(expr= m.x167 - 4.45628648004517 * m.b315 &lt;= 0)
m.e92 = Constraint(expr= m.x168 - 4.45628648004517 * m.b316 &lt;= 0)
m.e93 = Constraint(expr= m.x169 + 4.45628648004517 * m.b313
    &lt;= 4.45628648004517)
m.e94 = Constraint(expr= m.x170 + 4.45628648004517 * m.b314
    &lt;= 4.45628648004517)
m.e95 = Constraint(expr= m.x171 + 4.45628648004517 * m.b315
    &lt;= 4.45628648004517)
m.e96 = Constraint(expr= m.x172 + 4.45628648004517 * m.b316
    &lt;= 4.45628648004517)
m.e97 = Constraint(expr= -0.75 * m.x173 + m.x205 == 0)
m.e98 = Constraint(expr= -0.75 * m.x174 + m.x206 == 0)
m.e99 = Constraint(expr= -0.75 * m.x175 + m.x207 == 0)
m.e100 = Constraint(expr= -0.75 * m.x176 + m.x208 == 0)
m.e101 = Constraint(expr= m.x177 == 0)
m.e102 = Constraint(expr= m.x178 == 0)
m.e103 = Constraint(expr= m.x179 == 0)
m.e104 = Constraint(expr= m.x180 == 0)
m.e105 = Constraint(expr= m.x209 == 0)
m.e106 = Constraint(expr= m.x210 == 0)
m.e107 = Constraint(expr= m.x211 == 0)
m.e108 = Constraint(expr= m.x212 == 0)
m.e109 = Constraint(expr= m.x73 - m.x173 - m.x177 == 0)
m.e110 = Constraint(expr= m.x74 - m.x174 - m.x178 == 0)
m.e111 = Constraint(expr= m.x75 - m.x175 - m.x179 == 0)
m.e112 = Constraint(expr= m.x76 - m.x176 - m.x180 == 0)
m.e113 = Constraint(expr= m.x89 - m.x205 - m.x209 == 0)
m.e114 = Constraint(expr= m.x90 - m.x206 - m.x210 == 0)
m.e115 = Constraint(expr= m.x91 - m.x207 - m.x211 == 0)
m.e116 = Constraint(expr= m.x92 - m.x208 - m.x212 == 0)
m.e117 = Constraint(expr= m.x173 - 4.45628648004517 * m.b317 &lt;= 0)
m.e118 = Constraint(expr= m.x174 - 4.45628648004517 * m.b318 &lt;= 0)
m.e119 = Constraint(expr= m.x175 - 4.45628648004517 * m.b319 &lt;= 0)
m.e120 = Constraint(expr= m.x176 - 4.45628648004517 * m.b320 &lt;= 0)
m.e121 = Constraint(expr= m.x177 + 4.45628648004517 * m.b317
    &lt;= 4.45628648004517)
m.e122 = Constraint(expr= m.x178 + 4.45628648004517 * m.b318
    &lt;= 4.45628648004517)
m.e123 = Constraint(expr= m.x179 + 4.45628648004517 * m.b319
    &lt;= 4.45628648004517)
m.e124 = Constraint(expr= m.x180 + 4.45628648004517 * m.b320
    &lt;= 4.45628648004517)
m.e125 = Constraint(expr= m.x205 - 3.34221486003388 * m.b317 &lt;= 0)
m.e126 = Constraint(expr= m.x206 - 3.34221486003388 * m.b318 &lt;= 0)
m.e127 = Constraint(expr= m.x207 - 3.34221486003388 * m.b319 &lt;= 0)
m.e128 = Constraint(expr= m.x208 - 3.34221486003388 * m.b320 &lt;= 0)
m.e129 = Constraint(expr= m.x209 + 3.34221486003388 * m.b317
    &lt;= 3.34221486003388)
m.e130 = Constraint(expr= m.x210 + 3.34221486003388 * m.b318
    &lt;= 3.34221486003388)
m.e131 = Constraint(expr= m.x211 + 3.34221486003388 * m.b319
    &lt;= 3.34221486003388)
m.e132 = Constraint(expr= m.x212 + 3.34221486003388 * m.b320
    &lt;= 3.34221486003388)
m.e133 = Constraint(expr= (m.x213 / (1e-06 + m.b321) - 1.5 * log(m.x181 / (
    1e-06 + m.b321) + 1)) * (1e-06 + m.b321) &lt;= 0)
m.e134 = Constraint(expr= (m.x214 / (1e-06 + m.b322) - 1.5 * log(m.x182 / (
    1e-06 + m.b322) + 1)) * (1e-06 + m.b322) &lt;= 0)
m.e135 = Constraint(expr= (m.x215 / (1e-06 + m.b323) - 1.5 * log(m.x183 / (
    1e-06 + m.b323) + 1)) * (1e-06 + m.b323) &lt;= 0)
m.e136 = Constraint(expr= (m.x216 / (1e-06 + m.b324) - 1.5 * log(m.x184 / (
    1e-06 + m.b324) + 1)) * (1e-06 + m.b324) &lt;= 0)
m.e137 = Constraint(expr= m.x185 == 0)
m.e138 = Constraint(expr= m.x186 == 0)
m.e139 = Constraint(expr= m.x187 == 0)
m.e140 = Constraint(expr= m.x188 == 0)
m.e141 = Constraint(expr= m.x221 == 0)
m.e142 = Constraint(expr= m.x222 == 0)
m.e143 = Constraint(expr= m.x223 == 0)
m.e144 = Constraint(expr= m.x224 == 0)
m.e145 = Constraint(expr= m.x77 - m.x181 - m.x185 == 0)
m.e146 = Constraint(expr= m.x78 - m.x182 - m.x186 == 0)
m.e147 = Constraint(expr= m.x79 - m.x183 - m.x187 == 0)
m.e148 = Constraint(expr= m.x80 - m.x184 - m.x188 == 0)
m.e149 = Constraint(expr= m.x93 - m.x213 - m.x221 == 0)
m.e150 = Constraint(expr= m.x94 - m.x214 - m.x222 == 0)
m.e151 = Constraint(expr= m.x95 - m.x215 - m.x223 == 0)
m.e152 = Constraint(expr= m.x96 - m.x216 - m.x224 == 0)
m.e153 = Constraint(expr= m.x181 - 4.45628648004517 * m.b321 &lt;= 0)
m.e154 = Constraint(expr= m.x182 - 4.45628648004517 * m.b322 &lt;= 0)
m.e155 = Constraint(expr= m.x183 - 4.45628648004517 * m.b323 &lt;= 0)
m.e156 = Constraint(expr= m.x184 - 4.45628648004517 * m.b324 &lt;= 0)
m.e157 = Constraint(expr= m.x185 + 4.45628648004517 * m.b321
    &lt;= 4.45628648004517)
m.e158 = Constraint(expr= m.x186 + 4.45628648004517 * m.b322
    &lt;= 4.45628648004517)
m.e159 = Constraint(expr= m.x187 + 4.45628648004517 * m.b323
    &lt;= 4.45628648004517)
m.e160 = Constraint(expr= m.x188 + 4.45628648004517 * m.b324
    &lt;= 4.45628648004517)
m.e161 = Constraint(expr= m.x213 - 2.54515263975353 * m.b321 &lt;= 0)
m.e162 = Constraint(expr= m.x214 - 2.54515263975353 * m.b322 &lt;= 0)
m.e163 = Constraint(expr= m.x215 - 2.54515263975353 * m.b323 &lt;= 0)
m.e164 = Constraint(expr= m.x216 - 2.54515263975353 * m.b324 &lt;= 0)
m.e165 = Constraint(expr= m.x221 + 2.54515263975353 * m.b321
    &lt;= 2.54515263975353)
m.e166 = Constraint(expr= m.x222 + 2.54515263975353 * m.b322
    &lt;= 2.54515263975353)
m.e167 = Constraint(expr= m.x223 + 2.54515263975353 * m.b323
    &lt;= 2.54515263975353)
m.e168 = Constraint(expr= m.x224 + 2.54515263975353 * m.b324
    &lt;= 2.54515263975353)
m.e169 = Constraint(expr= -m.x189 + m.x229 == 0)
m.e170 = Constraint(expr= -m.x190 + m.x230 == 0)
m.e171 = Constraint(expr= -m.x191 + m.x231 == 0)
m.e172 = Constraint(expr= -m.x192 + m.x232 == 0)
m.e173 = Constraint(expr= -0.5 * m.x197 + m.x229 == 0)
m.e174 = Constraint(expr= -0.5 * m.x198 + m.x230 == 0)
m.e175 = Constraint(expr= -0.5 * m.x199 + m.x231 == 0)
m.e176 = Constraint(expr= -0.5 * m.x200 + m.x232 == 0)
m.e177 = Constraint(expr= m.x193 == 0)
m.e178 = Constraint(expr= m.x194 == 0)
m.e179 = Constraint(expr= m.x195 == 0)
m.e180 = Constraint(expr= m.x196 == 0)
m.e181 = Constraint(expr= m.x201 == 0)
m.e182 = Constraint(expr= m.x202 == 0)
m.e183 = Constraint(expr= m.x203 == 0)
m.e184 = Constraint(expr= m.x204 == 0)
m.e185 = Constraint(expr= m.x233 == 0)
m.e186 = Constraint(expr= m.x234 == 0)
m.e187 = Constraint(expr= m.x235 == 0)
m.e188 = Constraint(expr= m.x236 == 0)
m.e189 = Constraint(expr= m.x81 - m.x189 - m.x193 == 0)
m.e190 = Constraint(expr= m.x82 - m.x190 - m.x194 == 0)
m.e191 = Constraint(expr= m.x83 - m.x191 - m.x195 == 0)
m.e192 = Constraint(expr= m.x84 - m.x192 - m.x196 == 0)
m.e193 = Constraint(expr= m.x85 - m.x197 - m.x201 == 0)
m.e194 = Constraint(expr= m.x86 - m.x198 - m.x202 == 0)
m.e195 = Constraint(expr= m.x87 - m.x199 - m.x203 == 0)
m.e196 = Constraint(expr= m.x88 - m.x200 - m.x204 == 0)
m.e197 = Constraint(expr= m.x97 - m.x229 - m.x233 == 0)
m.e198 = Constraint(expr= m.x98 - m.x230 - m.x234 == 0)
m.e199 = Constraint(expr= m.x99 - m.x231 - m.x235 == 0)
m.e200 = Constraint(expr= m.x100 - m.x232 - m.x236 == 0)
m.e201 = Constraint(expr= m.x189 - 4.45628648004517 * m.b325 &lt;= 0)
m.e202 = Constraint(expr= m.x190 - 4.45628648004517 * m.b326 &lt;= 0)
m.e203 = Constraint(expr= m.x191 - 4.45628648004517 * m.b327 &lt;= 0)
m.e204 = Constraint(expr= m.x192 - 4.45628648004517 * m.b328 &lt;= 0)
m.e205 = Constraint(expr= m.x193 + 4.45628648004517 * m.b325
    &lt;= 4.45628648004517)
m.e206 = Constraint(expr= m.x194 + 4.45628648004517 * m.b326
    &lt;= 4.45628648004517)
m.e207 = Constraint(expr= m.x195 + 4.45628648004517 * m.b327
    &lt;= 4.45628648004517)
m.e208 = Constraint(expr= m.x196 + 4.45628648004517 * m.b328
    &lt;= 4.45628648004517)
m.e209 = Constraint(expr= m.x197 - 30 * m.b325 &lt;= 0)
m.e210 = Constraint(expr= m.x198 - 30 * m.b326 &lt;= 0)
m.e211 = Constraint(expr= m.x199 - 30 * m.b327 &lt;= 0)
m.e212 = Constraint(expr= m.x200 - 30 * m.b328 &lt;= 0)
m.e213 = Constraint(expr= m.x201 + 30 * m.b325 &lt;= 30)
m.e214 = Constraint(expr= m.x202 + 30 * m.b326 &lt;= 30)
m.e215 = Constraint(expr= m.x203 + 30 * m.b327 &lt;= 30)
m.e216 = Constraint(expr= m.x204 + 30 * m.b328 &lt;= 30)
m.e217 = Constraint(expr= m.x229 - 15 * m.b325 &lt;= 0)
m.e218 = Constraint(expr= m.x230 - 15 * m.b326 &lt;= 0)
m.e219 = Constraint(expr= m.x231 - 15 * m.b327 &lt;= 0)
m.e220 = Constraint(expr= m.x232 - 15 * m.b328 &lt;= 0)
m.e221 = Constraint(expr= m.x233 + 15 * m.b325 &lt;= 15)
m.e222 = Constraint(expr= m.x234 + 15 * m.b326 &lt;= 15)
m.e223 = Constraint(expr= m.x235 + 15 * m.b327 &lt;= 15)
m.e224 = Constraint(expr= m.x236 + 15 * m.b328 &lt;= 15)
m.e225 = Constraint(expr= (m.x269 / (1e-06 + m.b329) - 1.25 * log(m.x237 / (
    1e-06 + m.b329) + 1)) * (1e-06 + m.b329) &lt;= 0)
m.e226 = Constraint(expr= (m.x270 / (1e-06 + m.b330) - 1.25 * log(m.x238 / (
    1e-06 + m.b330) + 1)) * (1e-06 + m.b330) &lt;= 0)
m.e227 = Constraint(expr= (m.x271 / (1e-06 + m.b331) - 1.25 * log(m.x239 / (
    1e-06 + m.b331) + 1)) * (1e-06 + m.b331) &lt;= 0)
m.e228 = Constraint(expr= (m.x272 / (1e-06 + m.b332) - 1.25 * log(m.x240 / (
    1e-06 + m.b332) + 1)) * (1e-06 + m.b332) &lt;= 0)
m.e229 = Constraint(expr= m.x241 == 0)
m.e230 = Constraint(expr= m.x242 == 0)
m.e231 = Constraint(expr= m.x243 == 0)
m.e232 = Constraint(expr= m.x244 == 0)
m.e233 = Constraint(expr= m.x273 == 0)
m.e234 = Constraint(expr= m.x274 == 0)
m.e235 = Constraint(expr= m.x275 == 0)
m.e236 = Constraint(expr= m.x276 == 0)
m.e237 = Constraint(expr= m.x101 - m.x237 - m.x241 == 0)
m.e238 = Constraint(expr= m.x102 - m.x238 - m.x242 == 0)
m.e239 = Constraint(expr= m.x103 - m.x239 - m.x243 == 0)
m.e240 = Constraint(expr= m.x104 - m.x240 - m.x244 == 0)
m.e241 = Constraint(expr= m.x121 - m.x269 - m.x273 == 0)
m.e242 = Constraint(expr= m.x122 - m.x270 - m.x274 == 0)
m.e243 = Constraint(expr= m.x123 - m.x271 - m.x275 == 0)
m.e244 = Constraint(expr= m.x124 - m.x272 - m.x276 == 0)
m.e245 = Constraint(expr= m.x237 - 3.34221486003388 * m.b329 &lt;= 0)
m.e246 = Constraint(expr= m.x238 - 3.34221486003388 * m.b330 &lt;= 0)
m.e247 = Constraint(expr= m.x239 - 3.34221486003388 * m.b331 &lt;= 0)
m.e248 = Constraint(expr= m.x240 - 3.34221486003388 * m.b332 &lt;= 0)
m.e249 = Constraint(expr= m.x241 + 3.34221486003388 * m.b329
    &lt;= 3.34221486003388)
m.e250 = Constraint(expr= m.x242 + 3.34221486003388 * m.b330
    &lt;= 3.34221486003388)
m.e251 = Constraint(expr= m.x243 + 3.34221486003388 * m.b331
    &lt;= 3.34221486003388)
m.e252 = Constraint(expr= m.x244 + 3.34221486003388 * m.b332
    &lt;= 3.34221486003388)
m.e253 = Constraint(expr= m.x269 - 1.83548069293539 * m.b329 &lt;= 0)
m.e254 = Constraint(expr= m.x270 - 1.83548069293539 * m.b330 &lt;= 0)
m.e255 = Constraint(expr= m.x271 - 1.83548069293539 * m.b331 &lt;= 0)
m.e256 = Constraint(expr= m.x272 - 1.83548069293539 * m.b332 &lt;= 0)
m.e257 = Constraint(expr= m.x273 + 1.83548069293539 * m.b329
    &lt;= 1.83548069293539)
m.e258 = Constraint(expr= m.x274 + 1.83548069293539 * m.b330
    &lt;= 1.83548069293539)
m.e259 = Constraint(expr= m.x275 + 1.83548069293539 * m.b331
    &lt;= 1.83548069293539)
m.e260 = Constraint(expr= m.x276 + 1.83548069293539 * m.b332
    &lt;= 1.83548069293539)
m.e261 = Constraint(expr= (m.x277 / (1e-06 + m.b333) - 0.9 * log(m.x245 / (
    1e-06 + m.b333) + 1)) * (1e-06 + m.b333) &lt;= 0)
m.e262 = Constraint(expr= (m.x278 / (1e-06 + m.b334) - 0.9 * log(m.x246 / (
    1e-06 + m.b334) + 1)) * (1e-06 + m.b334) &lt;= 0)
m.e263 = Constraint(expr= (m.x279 / (1e-06 + m.b335) - 0.9 * log(m.x247 / (
    1e-06 + m.b335) + 1)) * (1e-06 + m.b335) &lt;= 0)
m.e264 = Constraint(expr= (m.x280 / (1e-06 + m.b336) - 0.9 * log(m.x248 / (
    1e-06 + m.b336) + 1)) * (1e-06 + m.b336) &lt;= 0)
m.e265 = Constraint(expr= m.x249 == 0)
m.e266 = Constraint(expr= m.x250 == 0)
m.e267 = Constraint(expr= m.x251 == 0)
m.e268 = Constraint(expr= m.x252 == 0)
m.e269 = Constraint(expr= m.x281 == 0)
m.e270 = Constraint(expr= m.x282 == 0)
m.e271 = Constraint(expr= m.x283 == 0)
m.e272 = Constraint(expr= m.x284 == 0)
m.e273 = Constraint(expr= m.x105 - m.x245 - m.x249 == 0)
m.e274 = Constraint(expr= m.x106 - m.x246 - m.x250 == 0)
m.e275 = Constraint(expr= m.x107 - m.x247 - m.x251 == 0)
m.e276 = Constraint(expr= m.x108 - m.x248 - m.x252 == 0)
m.e277 = Constraint(expr= m.x125 - m.x277 - m.x281 == 0)
m.e278 = Constraint(expr= m.x126 - m.x278 - m.x282 == 0)
m.e279 = Constraint(expr= m.x127 - m.x279 - m.x283 == 0)
m.e280 = Constraint(expr= m.x128 - m.x280 - m.x284 == 0)
m.e281 = Constraint(expr= m.x245 - 3.34221486003388 * m.b333 &lt;= 0)
m.e282 = Constraint(expr= m.x246 - 3.34221486003388 * m.b334 &lt;= 0)
m.e283 = Constraint(expr= m.x247 - 3.34221486003388 * m.b335 &lt;= 0)
m.e284 = Constraint(expr= m.x248 - 3.34221486003388 * m.b336 &lt;= 0)
m.e285 = Constraint(expr= m.x249 + 3.34221486003388 * m.b333
    &lt;= 3.34221486003388)
m.e286 = Constraint(expr= m.x250 + 3.34221486003388 * m.b334
    &lt;= 3.34221486003388)
m.e287 = Constraint(expr= m.x251 + 3.34221486003388 * m.b335
    &lt;= 3.34221486003388)
m.e288 = Constraint(expr= m.x252 + 3.34221486003388 * m.b336
    &lt;= 3.34221486003388)
m.e289 = Constraint(expr= m.x277 - 1.32154609891348 * m.b333 &lt;= 0)
m.e290 = Constraint(expr= m.x278 - 1.32154609891348 * m.b334 &lt;= 0)
m.e291 = Constraint(expr= m.x279 - 1.32154609891348 * m.b335 &lt;= 0)
m.e292 = Constraint(expr= m.x280 - 1.32154609891348 * m.b336 &lt;= 0)
m.e293 = Constraint(expr= m.x281 + 1.32154609891348 * m.b333
    &lt;= 1.32154609891348)
m.e294 = Constraint(expr= m.x282 + 1.32154609891348 * m.b334
    &lt;= 1.32154609891348)
m.e295 = Constraint(expr= m.x283 + 1.32154609891348 * m.b335
    &lt;= 1.32154609891348)
m.e296 = Constraint(expr= m.x284 + 1.32154609891348 * m.b336
    &lt;= 1.32154609891348)
m.e297 = Constraint(expr= (m.x285 / (1e-06 + m.b337) - log(m.x217 / (1e-06 +
    m.b337) + 1)) * (1e-06 + m.b337) &lt;= 0)
m.e298 = Constraint(expr= (m.x286 / (1e-06 + m.b338) - log(m.x218 / (1e-06 +
    m.b338) + 1)) * (1e-06 + m.b338) &lt;= 0)
m.e299 = Constraint(expr= (m.x287 / (1e-06 + m.b339) - log(m.x219 / (1e-06 +
    m.b339) + 1)) * (1e-06 + m.b339) &lt;= 0)
m.e300 = Constraint(expr= (m.x288 / (1e-06 + m.b340) - log(m.x220 / (1e-06 +
    m.b340) + 1)) * (1e-06 + m.b340) &lt;= 0)
m.e301 = Constraint(expr= m.x225 == 0)
m.e302 = Constraint(expr= m.x226 == 0)
m.e303 = Constraint(expr= m.x227 == 0)
m.e304 = Constraint(expr= m.x228 == 0)
m.e305 = Constraint(expr= m.x289 == 0)
m.e306 = Constraint(expr= m.x290 == 0)
m.e307 = Constraint(expr= m.x291 == 0)
m.e308 = Constraint(expr= m.x292 == 0)
m.e309 = Constraint(expr= m.x93 - m.x217 - m.x225 == 0)
m.e310 = Constraint(expr= m.x94 - m.x218 - m.x226 == 0)
m.e311 = Constraint(expr= m.x95 - m.x219 - m.x227 == 0)
m.e312 = Constraint(expr= m.x96 - m.x220 - m.x228 == 0)
m.e313 = Constraint(expr= m.x129 - m.x285 - m.x289 == 0)
m.e314 = Constraint(expr= m.x130 - m.x286 - m.x290 == 0)
m.e315 = Constraint(expr= m.x131 - m.x287 - m.x291 == 0)
m.e316 = Constraint(expr= m.x132 - m.x288 - m.x292 == 0)
m.e317 = Constraint(expr= m.x217 - 2.54515263975353 * m.b337 &lt;= 0)
m.e318 = Constraint(expr= m.x218 - 2.54515263975353 * m.b338 &lt;= 0)
m.e319 = Constraint(expr= m.x219 - 2.54515263975353 * m.b339 &lt;= 0)
m.e320 = Constraint(expr= m.x220 - 2.54515263975353 * m.b340 &lt;= 0)
m.e321 = Constraint(expr= m.x225 + 2.54515263975353 * m.b337
    &lt;= 2.54515263975353)
m.e322 = Constraint(expr= m.x226 + 2.54515263975353 * m.b338
    &lt;= 2.54515263975353)
m.e323 = Constraint(expr= m.x227 + 2.54515263975353 * m.b339
    &lt;= 2.54515263975353)
m.e324 = Constraint(expr= m.x228 + 2.54515263975353 * m.b340
    &lt;= 2.54515263975353)
m.e325 = Constraint(expr= m.x285 - 1.26558121681553 * m.b337 &lt;= 0)
m.e326 = Constraint(expr= m.x286 - 1.26558121681553 * m.b338 &lt;= 0)
m.e327 = Constraint(expr= m.x287 - 1.26558121681553 * m.b339 &lt;= 0)
m.e328 = Constraint(expr= m.x288 - 1.26558121681553 * m.b340 &lt;= 0)
m.e329 = Constraint(expr= m.x289 + 1.26558121681553 * m.b337
    &lt;= 1.26558121681553)
m.e330 = Constraint(expr= m.x290 + 1.26558121681553 * m.b338
    &lt;= 1.26558121681553)
m.e331 = Constraint(expr= m.x291 + 1.26558121681553 * m.b339
    &lt;= 1.26558121681553)
m.e332 = Constraint(expr= m.x292 + 1.26558121681553 * m.b340
    &lt;= 1.26558121681553)
m.e333 = Constraint(expr= -0.9 * m.x253 + m.x293 == 0)
m.e334 = Constraint(expr= -0.9 * m.x254 + m.x294 == 0)
m.e335 = Constraint(expr= -0.9 * m.x255 + m.x295 == 0)
m.e336 = Constraint(expr= -0.9 * m.x256 + m.x296 == 0)
m.e337 = Constraint(expr= m.x257 == 0)
m.e338 = Constraint(expr= m.x258 == 0)
m.e339 = Constraint(expr= m.x259 == 0)
m.e340 = Constraint(expr= m.x260 == 0)
m.e341 = Constraint(expr= m.x297 == 0)
m.e342 = Constraint(expr= m.x298 == 0)
m.e343 = Constraint(expr= m.x299 == 0)
m.e344 = Constraint(expr= m.x300 == 0)
m.e345 = Constraint(expr= m.x109 - m.x253 - m.x257 == 0)
m.e346 = Constraint(expr= m.x110 - m.x254 - m.x258 == 0)
m.e347 = Constraint(expr= m.x111 - m.x255 - m.x259 == 0)
m.e348 = Constraint(expr= m.x112 - m.x256 - m.x260 == 0)
m.e349 = Constraint(expr= m.x133 - m.x293 - m.x297 == 0)
m.e350 = Constraint(expr= m.x134 - m.x294 - m.x298 == 0)
m.e351 = Constraint(expr= m.x135 - m.x295 - m.x299 == 0)
m.e352 = Constraint(expr= m.x136 - m.x296 - m.x300 == 0)
m.e353 = Constraint(expr= m.x253 - 15 * m.b341 &lt;= 0)
m.e354 = Constraint(expr= m.x254 - 15 * m.b342 &lt;= 0)
m.e355 = Constraint(expr= m.x255 - 15 * m.b343 &lt;= 0)
m.e356 = Constraint(expr= m.x256 - 15 * m.b344 &lt;= 0)
m.e357 = Constraint(expr= m.x257 + 15 * m.b341 &lt;= 15)
m.e358 = Constraint(expr= m.x258 + 15 * m.b342 &lt;= 15)
m.e359 = Constraint(expr= m.x259 + 15 * m.b343 &lt;= 15)
m.e360 = Constraint(expr= m.x260 + 15 * m.b344 &lt;= 15)
m.e361 = Constraint(expr= m.x293 - 13.5 * m.b341 &lt;= 0)
m.e362 = Constraint(expr= m.x294 - 13.5 * m.b342 &lt;= 0)
m.e363 = Constraint(expr= m.x295 - 13.5 * m.b343 &lt;= 0)
m.e364 = Constraint(expr= m.x296 - 13.5 * m.b344 &lt;= 0)
m.e365 = Constraint(expr= m.x297 + 13.5 * m.b341 &lt;= 13.5)
m.e366 = Constraint(expr= m.x298 + 13.5 * m.b342 &lt;= 13.5)
m.e367 = Constraint(expr= m.x299 + 13.5 * m.b343 &lt;= 13.5)
m.e368 = Constraint(expr= m.x300 + 13.5 * m.b344 &lt;= 13.5)
m.e369 = Constraint(expr= -0.6 * m.x261 + m.x301 == 0)
m.e370 = Constraint(expr= -0.6 * m.x262 + m.x302 == 0)
m.e371 = Constraint(expr= -0.6 * m.x263 + m.x303 == 0)
m.e372 = Constraint(expr= -0.6 * m.x264 + m.x304 == 0)
m.e373 = Constraint(expr= m.x265 == 0)
m.e374 = Constraint(expr= m.x266 == 0)
m.e375 = Constraint(expr= m.x267 == 0)
m.e376 = Constraint(expr= m.x268 == 0)
m.e377 = Constraint(expr= m.x305 == 0)
m.e378 = Constraint(expr= m.x306 == 0)
m.e379 = Constraint(expr= m.x307 == 0)
m.e380 = Constraint(expr= m.x308 == 0)
m.e381 = Constraint(expr= m.x113 - m.x261 - m.x265 == 0)
m.e382 = Constraint(expr= m.x114 - m.x262 - m.x266 == 0)
m.e383 = Constraint(expr= m.x115 - m.x263 - m.x267 == 0)
m.e384 = Constraint(expr= m.x116 - m.x264 - m.x268 == 0)
m.e385 = Constraint(expr= m.x137 - m.x301 - m.x305 == 0)
m.e386 = Constraint(expr= m.x138 - m.x302 - m.x306 == 0)
m.e387 = Constraint(expr= m.x139 - m.x303 - m.x307 == 0)
m.e388 = Constraint(expr= m.x140 - m.x304 - m.x308 == 0)
m.e389 = Constraint(expr= m.x261 - 15 * m.b345 &lt;= 0)
m.e390 = Constraint(expr= m.x262 - 15 * m.b346 &lt;= 0)
m.e391 = Constraint(expr= m.x263 - 15 * m.b347 &lt;= 0)
m.e392 = Constraint(expr= m.x264 - 15 * m.b348 &lt;= 0)
m.e393 = Constraint(expr= m.x265 + 15 * m.b345 &lt;= 15)
m.e394 = Constraint(expr= m.x266 + 15 * m.b346 &lt;= 15)
m.e395 = Constraint(expr= m.x267 + 15 * m.b347 &lt;= 15)
m.e396 = Constraint(expr= m.x268 + 15 * m.b348 &lt;= 15)
m.e397 = Constraint(expr= m.x301 - 9 * m.b345 &lt;= 0)
m.e398 = Constraint(expr= m.x302 - 9 * m.b346 &lt;= 0)
m.e399 = Constraint(expr= m.x303 - 9 * m.b347 &lt;= 0)
m.e400 = Constraint(expr= m.x304 - 9 * m.b348 &lt;= 0)
m.e401 = Constraint(expr= m.x305 + 9 * m.b345 &lt;= 9)
m.e402 = Constraint(expr= m.x306 + 9 * m.b346 &lt;= 9)
m.e403 = Constraint(expr= m.x307 + 9 * m.b347 &lt;= 9)
m.e404 = Constraint(expr= m.x308 + 9 * m.b348 &lt;= 9)
m.e405 = Constraint(expr= m.x1 + 5 * m.b349 == 0)
m.e406 = Constraint(expr= m.x2 + 4 * m.b350 == 0)
m.e407 = Constraint(expr= m.x3 + 6 * m.b351 == 0)
m.e408 = Constraint(expr= m.x4 + 3 * m.b352 == 0)
m.e409 = Constraint(expr= m.x5 + 8 * m.b353 == 0)
m.e410 = Constraint(expr= m.x6 + 7 * m.b354 == 0)
m.e411 = Constraint(expr= m.x7 + 6 * m.b355 == 0)
m.e412 = Constraint(expr= m.x8 + 5 * m.b356 == 0)
m.e413 = Constraint(expr= m.x9 + 6 * m.b357 == 0)
m.e414 = Constraint(expr= m.x10 + 9 * m.b358 == 0)
m.e415 = Constraint(expr= m.x11 + 4 * m.b359 == 0)
m.e416 = Constraint(expr= m.x12 + 3 * m.b360 == 0)
m.e417 = Constraint(expr= m.x13 + 10 * m.b361 == 0)
m.e418 = Constraint(expr= m.x14 + 9 * m.b362 == 0)
m.e419 = Constraint(expr= m.x15 + 5 * m.b363 == 0)
m.e420 = Constraint(expr= m.x16 + 6 * m.b364 == 0)
m.e421 = Constraint(expr= m.x17 + 6 * m.b365 == 0)
m.e422 = Constraint(expr= m.x18 + 10 * m.b366 == 0)
m.e423 = Constraint(expr= m.x19 + 6 * m.b367 == 0)
m.e424 = Constraint(expr= m.x20 + 9 * m.b368 == 0)
m.e425 = Constraint(expr= m.x21 + 7 * m.b369 == 0)
m.e426 = Constraint(expr= m.x22 + 7 * m.b370 == 0)
m.e427 = Constraint(expr= m.x23 + 4 * m.b371 == 0)
m.e428 = Constraint(expr= m.x24 + 2 * m.b372 == 0)
m.e429 = Constraint(expr= m.x25 + 4 * m.b373 == 0)
m.e430 = Constraint(expr= m.x26 + 3 * m.b374 == 0)
m.e431 = Constraint(expr= m.x27 + 2 * m.b375 == 0)
m.e432 = Constraint(expr= m.x28 + 8 * m.b376 == 0)
m.e433 = Constraint(expr= m.x29 + 5 * m.b377 == 0)
m.e434 = Constraint(expr= m.x30 + 6 * m.b378 == 0)
m.e435 = Constraint(expr= m.x31 + 7 * m.b379 == 0)
m.e436 = Constraint(expr= m.x32 + 4 * m.b380 == 0)
m.e437 = Constraint(expr= m.x33 + 2 * m.b381 == 0)
m.e438 = Constraint(expr= m.x34 + 5 * m.b382 == 0)
m.e439 = Constraint(expr= m.x35 + 2 * m.b383 == 0)
m.e440 = Constraint(expr= m.x36 + 6 * m.b384 == 0)
m.e441 = Constraint(expr= m.x37 + 4 * m.b385 == 0)
m.e442 = Constraint(expr= m.x38 + 7 * m.b386 == 0)
m.e443 = Constraint(expr= m.x39 + 4 * m.b387 == 0)
m.e444 = Constraint(expr= m.x40 + 7 * m.b388 == 0)
m.e445 = Constraint(expr= m.b309 - m.b310 &lt;= 0)
m.e446 = Constraint(expr= m.b309 - m.b311 &lt;= 0)
m.e447 = Constraint(expr= m.b309 - m.b312 &lt;= 0)
m.e448 = Constraint(expr= m.b310 - m.b311 &lt;= 0)
m.e449 = Constraint(expr= m.b310 - m.b312 &lt;= 0)
m.e450 = Constraint(expr= m.b311 - m.b312 &lt;= 0)
m.e451 = Constraint(expr= m.b313 - m.b314 &lt;= 0)
m.e452 = Constraint(expr= m.b313 - m.b315 &lt;= 0)
m.e453 = Constraint(expr= m.b313 - m.b316 &lt;= 0)
m.e454 = Constraint(expr= m.b314 - m.b315 &lt;= 0)
m.e455 = Constraint(expr= m.b314 - m.b316 &lt;= 0)
m.e456 = Constraint(expr= m.b315 - m.b316 &lt;= 0)
m.e457 = Constraint(expr= m.b317 - m.b318 &lt;= 0)
m.e458 = Constraint(expr= m.b317 - m.b319 &lt;= 0)
m.e459 = Constraint(expr= m.b317 - m.b320 &lt;= 0)
m.e460 = Constraint(expr= m.b318 - m.b319 &lt;= 0)
m.e461 = Constraint(expr= m.b318 - m.b320 &lt;= 0)
m.e462 = Constraint(expr= m.b319 - m.b320 &lt;= 0)
m.e463 = Constraint(expr= m.b321 - m.b322 &lt;= 0)
m.e464 = Constraint(expr= m.b321 - m.b323 &lt;= 0)
m.e465 = Constraint(expr= m.b321 - m.b324 &lt;= 0)
m.e466 = Constraint(expr= m.b322 - m.b323 &lt;= 0)
m.e467 = Constraint(expr= m.b322 - m.b324 &lt;= 0)
m.e468 = Constraint(expr= m.b323 - m.b324 &lt;= 0)
m.e469 = Constraint(expr= m.b325 - m.b326 &lt;= 0)
m.e470 = Constraint(expr= m.b325 - m.b327 &lt;= 0)
m.e471 = Constraint(expr= m.b325 - m.b328 &lt;= 0)
m.e472 = Constraint(expr= m.b326 - m.b327 &lt;= 0)
m.e473 = Constraint(expr= m.b326 - m.b328 &lt;= 0)
m.e474 = Constraint(expr= m.b327 - m.b328 &lt;= 0)
m.e475 = Constraint(expr= m.b329 - m.b330 &lt;= 0)
m.e476 = Constraint(expr= m.b329 - m.b331 &lt;= 0)
m.e477 = Constraint(expr= m.b329 - m.b332 &lt;= 0)
m.e478 = Constraint(expr= m.b330 - m.b331 &lt;= 0)
m.e479 = Constraint(expr= m.b330 - m.b332 &lt;= 0)
m.e480 = Constraint(expr= m.b331 - m.b332 &lt;= 0)
m.e481 = Constraint(expr= m.b333 - m.b334 &lt;= 0)
m.e482 = Constraint(expr= m.b333 - m.b335 &lt;= 0)
m.e483 = Constraint(expr= m.b333 - m.b336 &lt;= 0)
m.e484 = Constraint(expr= m.b334 - m.b335 &lt;= 0)
m.e485 = Constraint(expr= m.b334 - m.b336 &lt;= 0)
m.e486 = Constraint(expr= m.b335 - m.b336 &lt;= 0)
m.e487 = Constraint(expr= m.b337 - m.b338 &lt;= 0)
m.e488 = Constraint(expr= m.b337 - m.b339 &lt;= 0)
m.e489 = Constraint(expr= m.b337 - m.b340 &lt;= 0)
m.e490 = Constraint(expr= m.b338 - m.b339 &lt;= 0)
m.e491 = Constraint(expr= m.b338 - m.b340 &lt;= 0)
m.e492 = Constraint(expr= m.b339 - m.b340 &lt;= 0)
m.e493 = Constraint(expr= m.b341 - m.b342 &lt;= 0)
m.e494 = Constraint(expr= m.b341 - m.b343 &lt;= 0)
m.e495 = Constraint(expr= m.b341 - m.b344 &lt;= 0)
m.e496 = Constraint(expr= m.b342 - m.b343 &lt;= 0)
m.e497 = Constraint(expr= m.b342 - m.b344 &lt;= 0)
m.e498 = Constraint(expr= m.b343 - m.b344 &lt;= 0)
m.e499 = Constraint(expr= m.b345 - m.b346 &lt;= 0)
m.e500 = Constraint(expr= m.b345 - m.b347 &lt;= 0)
m.e501 = Constraint(expr= m.b345 - m.b348 &lt;= 0)
m.e502 = Constraint(expr= m.b346 - m.b347 &lt;= 0)
m.e503 = Constraint(expr= m.b346 - m.b348 &lt;= 0)
m.e504 = Constraint(expr= m.b347 - m.b348 &lt;= 0)
m.e505 = Constraint(expr= m.b349 + m.b350 &lt;= 1)
m.e506 = Constraint(expr= m.b349 + m.b351 &lt;= 1)
m.e507 = Constraint(expr= m.b349 + m.b352 &lt;= 1)
m.e508 = Constraint(expr= m.b349 + m.b350 &lt;= 1)
m.e509 = Constraint(expr= m.b350 + m.b351 &lt;= 1)
m.e510 = Constraint(expr= m.b350 + m.b352 &lt;= 1)
m.e511 = Constraint(expr= m.b349 + m.b351 &lt;= 1)
m.e512 = Constraint(expr= m.b350 + m.b351 &lt;= 1)
m.e513 = Constraint(expr= m.b351 + m.b352 &lt;= 1)
m.e514 = Constraint(expr= m.b349 + m.b352 &lt;= 1)
m.e515 = Constraint(expr= m.b350 + m.b352 &lt;= 1)
m.e516 = Constraint(expr= m.b351 + m.b352 &lt;= 1)
m.e517 = Constraint(expr= m.b353 + m.b354 &lt;= 1)
m.e518 = Constraint(expr= m.b353 + m.b355 &lt;= 1)
m.e519 = Constraint(expr= m.b353 + m.b356 &lt;= 1)
m.e520 = Constraint(expr= m.b353 + m.b354 &lt;= 1)
m.e521 = Constraint(expr= m.b354 + m.b355 &lt;= 1)
m.e522 = Constraint(expr= m.b354 + m.b356 &lt;= 1)
m.e523 = Constraint(expr= m.b353 + m.b355 &lt;= 1)
m.e524 = Constraint(expr= m.b354 + m.b355 &lt;= 1)
m.e525 = Constraint(expr= m.b355 + m.b356 &lt;= 1)
m.e526 = Constraint(expr= m.b353 + m.b356 &lt;= 1)
m.e527 = Constraint(expr= m.b354 + m.b356 &lt;= 1)
m.e528 = Constraint(expr= m.b355 + m.b356 &lt;= 1)
m.e529 = Constraint(expr= m.b357 + m.b358 &lt;= 1)
m.e530 = Constraint(expr= m.b357 + m.b359 &lt;= 1)
m.e531 = Constraint(expr= m.b357 + m.b360 &lt;= 1)
m.e532 = Constraint(expr= m.b357 + m.b358 &lt;= 1)
m.e533 = Constraint(expr= m.b358 + m.b359 &lt;= 1)
m.e534 = Constraint(expr= m.b358 + m.b360 &lt;= 1)
m.e535 = Constraint(expr= m.b357 + m.b359 &lt;= 1)
m.e536 = Constraint(expr= m.b358 + m.b359 &lt;= 1)
m.e537 = Constraint(expr= m.b359 + m.b360 &lt;= 1)
m.e538 = Constraint(expr= m.b357 + m.b360 &lt;= 1)
m.e539 = Constraint(expr= m.b358 + m.b360 &lt;= 1)
m.e540 = Constraint(expr= m.b359 + m.b360 &lt;= 1)
m.e541 = Constraint(expr= m.b361 + m.b362 &lt;= 1)
m.e542 = Constraint(expr= m.b361 + m.b363 &lt;= 1)
m.e543 = Constraint(expr= m.b361 + m.b364 &lt;= 1)
m.e544 = Constraint(expr= m.b361 + m.b362 &lt;= 1)
m.e545 = Constraint(expr= m.b362 + m.b363 &lt;= 1)
m.e546 = Constraint(expr= m.b362 + m.b364 &lt;= 1)
m.e547 = Constraint(expr= m.b361 + m.b363 &lt;= 1)
m.e548 = Constraint(expr= m.b362 + m.b363 &lt;= 1)
m.e549 = Constraint(expr= m.b363 + m.b364 &lt;= 1)
m.e550 = Constraint(expr= m.b361 + m.b364 &lt;= 1)
m.e551 = Constraint(expr= m.b362 + m.b364 &lt;= 1)
m.e552 = Constraint(expr= m.b363 + m.b364 &lt;= 1)
m.e553 = Constraint(expr= m.b365 + m.b366 &lt;= 1)
m.e554 = Constraint(expr= m.b365 + m.b367 &lt;= 1)
m.e555 = Constraint(expr= m.b365 + m.b368 &lt;= 1)
m.e556 = Constraint(expr= m.b365 + m.b366 &lt;= 1)
m.e557 = Constraint(expr= m.b366 + m.b367 &lt;= 1)
m.e558 = Constraint(expr= m.b366 + m.b368 &lt;= 1)
m.e559 = Constraint(expr= m.b365 + m.b367 &lt;= 1)
m.e560 = Constraint(expr= m.b366 + m.b367 &lt;= 1)
m.e561 = Constraint(expr= m.b367 + m.b368 &lt;= 1)
m.e562 = Constraint(expr= m.b365 + m.b368 &lt;= 1)
m.e563 = Constraint(expr= m.b366 + m.b368 &lt;= 1)
m.e564 = Constraint(expr= m.b367 + m.b368 &lt;= 1)
m.e565 = Constraint(expr= m.b369 + m.b370 &lt;= 1)
m.e566 = Constraint(expr= m.b369 + m.b371 &lt;= 1)
m.e567 = Constraint(expr= m.b369 + m.b372 &lt;= 1)
m.e568 = Constraint(expr= m.b369 + m.b370 &lt;= 1)
m.e569 = Constraint(expr= m.b370 + m.b371 &lt;= 1)
m.e570 = Constraint(expr= m.b370 + m.b372 &lt;= 1)
m.e571 = Constraint(expr= m.b369 + m.b371 &lt;= 1)
m.e572 = Constraint(expr= m.b370 + m.b371 &lt;= 1)
m.e573 = Constraint(expr= m.b371 + m.b372 &lt;= 1)
m.e574 = Constraint(expr= m.b369 + m.b372 &lt;= 1)
m.e575 = Constraint(expr= m.b370 + m.b372 &lt;= 1)
m.e576 = Constraint(expr= m.b371 + m.b372 &lt;= 1)
m.e577 = Constraint(expr= m.b373 + m.b374 &lt;= 1)
m.e578 = Constraint(expr= m.b373 + m.b375 &lt;= 1)
m.e579 = Constraint(expr= m.b373 + m.b376 &lt;= 1)
m.e580 = Constraint(expr= m.b373 + m.b374 &lt;= 1)
m.e581 = Constraint(expr= m.b374 + m.b375 &lt;= 1)
m.e582 = Constraint(expr= m.b374 + m.b376 &lt;= 1)
m.e583 = Constraint(expr= m.b373 + m.b375 &lt;= 1)
m.e584 = Constraint(expr= m.b374 + m.b375 &lt;= 1)
m.e585 = Constraint(expr= m.b375 + m.b376 &lt;= 1)
m.e586 = Constraint(expr= m.b373 + m.b376 &lt;= 1)
m.e587 = Constraint(expr= m.b374 + m.b376 &lt;= 1)
m.e588 = Constraint(expr= m.b375 + m.b376 &lt;= 1)
m.e589 = Constraint(expr= m.b377 + m.b378 &lt;= 1)
m.e590 = Constraint(expr= m.b377 + m.b379 &lt;= 1)
m.e591 = Constraint(expr= m.b377 + m.b380 &lt;= 1)
m.e592 = Constraint(expr= m.b377 + m.b378 &lt;= 1)
m.e593 = Constraint(expr= m.b378 + m.b379 &lt;= 1)
m.e594 = Constraint(expr= m.b378 + m.b380 &lt;= 1)
m.e595 = Constraint(expr= m.b377 + m.b379 &lt;= 1)
m.e596 = Constraint(expr= m.b378 + m.b379 &lt;= 1)
m.e597 = Constraint(expr= m.b379 + m.b380 &lt;= 1)
m.e598 = Constraint(expr= m.b377 + m.b380 &lt;= 1)
m.e599 = Constraint(expr= m.b378 + m.b380 &lt;= 1)
m.e600 = Constraint(expr= m.b379 + m.b380 &lt;= 1)
m.e601 = Constraint(expr= m.b381 + m.b382 &lt;= 1)
m.e602 = Constraint(expr= m.b381 + m.b383 &lt;= 1)
m.e603 = Constraint(expr= m.b381 + m.b384 &lt;= 1)
m.e604 = Constraint(expr= m.b381 + m.b382 &lt;= 1)
m.e605 = Constraint(expr= m.b382 + m.b383 &lt;= 1)
m.e606 = Constraint(expr= m.b382 + m.b384 &lt;= 1)
m.e607 = Constraint(expr= m.b381 + m.b383 &lt;= 1)
m.e608 = Constraint(expr= m.b382 + m.b383 &lt;= 1)
m.e609 = Constraint(expr= m.b383 + m.b384 &lt;= 1)
m.e610 = Constraint(expr= m.b381 + m.b384 &lt;= 1)
m.e611 = Constraint(expr= m.b382 + m.b384 &lt;= 1)
m.e612 = Constraint(expr= m.b383 + m.b384 &lt;= 1)
m.e613 = Constraint(expr= m.b385 + m.b386 &lt;= 1)
m.e614 = Constraint(expr= m.b385 + m.b387 &lt;= 1)
m.e615 = Constraint(expr= m.b385 + m.b388 &lt;= 1)
m.e616 = Constraint(expr= m.b385 + m.b386 &lt;= 1)
m.e617 = Constraint(expr= m.b386 + m.b387 &lt;= 1)
m.e618 = Constraint(expr= m.b386 + m.b388 &lt;= 1)
m.e619 = Constraint(expr= m.b385 + m.b387 &lt;= 1)
m.e620 = Constraint(expr= m.b386 + m.b387 &lt;= 1)
m.e621 = Constraint(expr= m.b387 + m.b388 &lt;= 1)
m.e622 = Constraint(expr= m.b385 + m.b388 &lt;= 1)
m.e623 = Constraint(expr= m.b386 + m.b388 &lt;= 1)
m.e624 = Constraint(expr= m.b387 + m.b388 &lt;= 1)
m.e625 = Constraint(expr= m.b309 - m.b349 &lt;= 0)
m.e626 = Constraint(expr= -m.b309 + m.b310 - m.b350 &lt;= 0)
m.e627 = Constraint(expr= -m.b309 - m.b310 + m.b311 - m.b351 &lt;= 0)
m.e628 = Constraint(expr= -m.b309 - m.b310 - m.b311 + m.b312 - m.b352 &lt;= 0)
m.e629 = Constraint(expr= m.b313 - m.b353 &lt;= 0)
m.e630 = Constraint(expr= -m.b313 + m.b314 - m.b354 &lt;= 0)
m.e631 = Constraint(expr= -m.b313 - m.b314 + m.b315 - m.b355 &lt;= 0)
m.e632 = Constraint(expr= -m.b313 - m.b314 - m.b315 + m.b316 - m.b356 &lt;= 0)
m.e633 = Constraint(expr= m.b317 - m.b357 &lt;= 0)
m.e634 = Constraint(expr= -m.b317 + m.b318 - m.b358 &lt;= 0)
m.e635 = Constraint(expr= -m.b317 - m.b318 + m.b319 - m.b359 &lt;= 0)
m.e636 = Constraint(expr= -m.b317 - m.b318 - m.b319 + m.b320 - m.b360 &lt;= 0)
m.e637 = Constraint(expr= m.b321 - m.b361 &lt;= 0)
m.e638 = Constraint(expr= -m.b321 + m.b322 - m.b362 &lt;= 0)
m.e639 = Constraint(expr= -m.b321 - m.b322 + m.b323 - m.b363 &lt;= 0)
m.e640 = Constraint(expr= -m.b321 - m.b322 - m.b323 + m.b324 - m.b364 &lt;= 0)
m.e641 = Constraint(expr= m.b325 - m.b365 &lt;= 0)
m.e642 = Constraint(expr= -m.b325 + m.b326 - m.b366 &lt;= 0)
m.e643 = Constraint(expr= -m.b325 - m.b326 + m.b327 - m.b367 &lt;= 0)
m.e644 = Constraint(expr= -m.b325 - m.b326 - m.b327 + m.b328 - m.b368 &lt;= 0)
m.e645 = Constraint(expr= m.b329 - m.b369 &lt;= 0)
m.e646 = Constraint(expr= -m.b329 + m.b330 - m.b370 &lt;= 0)
m.e647 = Constraint(expr= -m.b329 - m.b330 + m.b331 - m.b371 &lt;= 0)
m.e648 = Constraint(expr= -m.b329 - m.b330 - m.b331 + m.b332 - m.b372 &lt;= 0)
m.e649 = Constraint(expr= m.b333 - m.b373 &lt;= 0)
m.e650 = Constraint(expr= -m.b333 + m.b334 - m.b374 &lt;= 0)
m.e651 = Constraint(expr= -m.b333 - m.b334 + m.b335 - m.b375 &lt;= 0)
m.e652 = Constraint(expr= -m.b333 - m.b334 - m.b335 + m.b336 - m.b376 &lt;= 0)
m.e653 = Constraint(expr= m.b337 - m.b377 &lt;= 0)
m.e654 = Constraint(expr= -m.b337 + m.b338 - m.b378 &lt;= 0)
m.e655 = Constraint(expr= -m.b337 - m.b338 + m.b339 - m.b379 &lt;= 0)
m.e656 = Constraint(expr= -m.b337 - m.b338 - m.b339 + m.b340 - m.b380 &lt;= 0)
m.e657 = Constraint(expr= m.b341 - m.b381 &lt;= 0)
m.e658 = Constraint(expr= -m.b341 + m.b342 - m.b382 &lt;= 0)
m.e659 = Constraint(expr= -m.b341 - m.b342 + m.b343 - m.b383 &lt;= 0)
m.e660 = Constraint(expr= -m.b341 - m.b342 - m.b343 + m.b344 - m.b384 &lt;= 0)
m.e661 = Constraint(expr= m.b345 - m.b385 &lt;= 0)
m.e662 = Constraint(expr= -m.b345 + m.b346 - m.b386 &lt;= 0)
m.e663 = Constraint(expr= -m.b345 - m.b346 + m.b347 - m.b387 &lt;= 0)
m.e664 = Constraint(expr= -m.b345 - m.b346 - m.b347 + m.b348 - m.b388 &lt;= 0)
m.e665 = Constraint(expr= m.b309 + m.b313 == 1)
m.e666 = Constraint(expr= m.b310 + m.b314 == 1)
m.e667 = Constraint(expr= m.b311 + m.b315 == 1)
m.e668 = Constraint(expr= m.b312 + m.b316 == 1)
m.e669 = Constraint(expr= -m.b317 + m.b329 + m.b333 &gt;= 0)
m.e670 = Constraint(expr= -m.b318 + m.b330 + m.b334 &gt;= 0)
m.e671 = Constraint(expr= -m.b319 + m.b331 + m.b335 &gt;= 0)
m.e672 = Constraint(expr= -m.b320 + m.b332 + m.b336 &gt;= 0)
m.e673 = Constraint(expr= -m.b321 + m.b337 &gt;= 0)
m.e674 = Constraint(expr= -m.b322 + m.b338 &gt;= 0)
m.e675 = Constraint(expr= -m.b323 + m.b339 &gt;= 0)
m.e676 = Constraint(expr= -m.b324 + m.b340 &gt;= 0)
m.e677 = Constraint(expr= m.b309 + m.b313 - m.b317 &gt;= 0)
m.e678 = Constraint(expr= m.b310 + m.b314 - m.b318 &gt;= 0)
m.e679 = Constraint(expr= m.b311 + m.b315 - m.b319 &gt;= 0)
m.e680 = Constraint(expr= m.b312 + m.b316 - m.b320 &gt;= 0)
m.e681 = Constraint(expr= m.b309 + m.b313 - m.b321 &gt;= 0)
m.e682 = Constraint(expr= m.b310 + m.b314 - m.b322 &gt;= 0)
m.e683 = Constraint(expr= m.b311 + m.b315 - m.b323 &gt;= 0)
m.e684 = Constraint(expr= m.b312 + m.b316 - m.b324 &gt;= 0)
m.e685 = Constraint(expr= m.b309 + m.b313 - m.b325 &gt;= 0)
m.e686 = Constraint(expr= m.b310 + m.b314 - m.b326 &gt;= 0)
m.e687 = Constraint(expr= m.b311 + m.b315 - m.b327 &gt;= 0)
m.e688 = Constraint(expr= m.b312 + m.b316 - m.b328 &gt;= 0)
m.e689 = Constraint(expr= m.b317 - m.b329 &gt;= 0)
m.e690 = Constraint(expr= m.b318 - m.b330 &gt;= 0)
m.e691 = Constraint(expr= m.b319 - m.b331 &gt;= 0)
m.e692 = Constraint(expr= m.b320 - m.b332 &gt;= 0)
m.e693 = Constraint(expr= m.b317 - m.b333 &gt;= 0)
m.e694 = Constraint(expr= m.b318 - m.b334 &gt;= 0)
m.e695 = Constraint(expr= m.b319 - m.b335 &gt;= 0)
m.e696 = Constraint(expr= m.b320 - m.b336 &gt;= 0)
m.e697 = Constraint(expr= m.b321 - m.b337 &gt;= 0)
m.e698 = Constraint(expr= m.b322 - m.b338 &gt;= 0)
m.e699 = Constraint(expr= m.b323 - m.b339 &gt;= 0)
m.e700 = Constraint(expr= m.b324 - m.b340 &gt;= 0)
m.e701 = Constraint(expr= m.b325 - m.b341 &gt;= 0)
m.e702 = Constraint(expr= m.b326 - m.b342 &gt;= 0)
m.e703 = Constraint(expr= m.b327 - m.b343 &gt;= 0)
m.e704 = Constraint(expr= m.b328 - m.b344 &gt;= 0)
m.e705 = Constraint(expr= m.b325 - m.b345 &gt;= 0)
m.e706 = Constraint(expr= m.b326 - m.b346 &gt;= 0)
m.e707 = Constraint(expr= m.b327 - m.b347 &gt;= 0)
m.e708 = Constraint(expr= m.b328 - m.b348 &gt;= 0)
</pre></body></html>