# MINLP written by GAMS Convert at 02/17/22 17:18:33 # # Equation counts # Total E G L N X C B # 284 122 96 66 0 0 0 0 # # Variable counts # x b i s1s s2s sc si # Total cont binary integer sos1 sos2 scont sint # 296 208 88 0 0 0 0 0 # FX 17 # # Nonzero counts # Total const NL # 1240 750 490 # # Reformulation has removed 1 variable and 1 equation from pyomo.environ import * model = m = ConcreteModel() m.x1 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=2.99573227355399) m.x2 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=2.99573227355399) m.x3 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=2.94443897916644) m.x4 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=2.94443897916644) m.x5 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=0) m.x6 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=1.9677597942808) m.x7 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=1.27461261372086) m.x8 = Var(within=Reals, bounds=(-6.90775527898214,3.68887945411394) , initialize=-1.1298234508907) m.x9 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=0) m.x10 = Var(within=Reals, bounds=(-1.6094379124341,3.68887945411394) , initialize=2.0433602879393) m.x11 = Var(within=Reals, bounds=(-6.90775527898214,-0.693147180559945) , initialize=-0.693147180559945) m.x12 = Var(within=Reals, bounds=(0.693147180559945,6.90775527898214) , initialize=0.693147180559945) m.x13 = Var(within=Reals, bounds=(-0.693147180559945,2.99573227355399) , initialize=0) m.x14 = Var(within=Reals, bounds=(0.693147180559945,6.90775527898214) , initialize=0.693147180559945) m.x15 = Var(within=Reals, bounds=(0.693147180559945,6.90775527898214) , initialize=0.693147180559945) m.x16 = Var(within=Reals, bounds=(0.693147180559945,6.90775527898214) , initialize=0.693147180559945) m.x17 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x18 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x19 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x20 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x21 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x22 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x23 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x24 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x25 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x26 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x27 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x28 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x29 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x30 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x31 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x32 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x33 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x34 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x35 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x36 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x37 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x38 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x39 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x40 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.61614353244073) m.x41 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x42 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x43 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x44 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x45 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x46 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x47 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x48 = Var(within=Reals, bounds=(0,5.29831736654804) , initialize=2.94443897916644) m.x49 = Var(within=Reals, bounds=(None,None), initialize=3.17756893714648) m.x50 = Var(within=Reals, bounds=(None,None), initialize=0.233622804898907) m.x51 = Var(within=Reals, bounds=(None,None), initialize=0.399415336353673) m.x52 = Var(within=Reals, bounds=(None,None), initialize=0.577807612471274) m.x53 = Var(within=Reals, bounds=(None,None), initialize=0.0827635536142284) m.x54 = Var(within=Reals, bounds=(0.587786664902119,0.587786664902119) , initialize=0.587786664902119) m.x55 = Var(within=Reals, bounds=(None,None), initialize=-1.07858898697236) m.x56 = Var(within=Reals, bounds=(None,None), initialize=-0.604447984103336) m.x57 = Var(within=Reals, bounds=(None,None), initialize=3.17756893714648) m.x58 = Var(within=Reals, bounds=(None,None), initialize=0.591132698406824) m.x59 = Var(within=Reals, bounds=(None,None), initialize=0) m.x60 = Var(within=Reals, bounds=(None,None), initialize=0) m.x61 = Var(within=Reals, bounds=(None,None), initialize=0.516209904652203) m.x62 = Var(within=Reals, bounds=(0.587786664902119,0.587786664902119) , initialize=0.587786664902119) m.x63 = Var(within=Reals, bounds=(None,None), initialize=-1.1469821150246) m.x64 = Var(within=Reals, bounds=(None,None), initialize=-0.689448619982603) m.x65 = Var(within=Reals, bounds=(None,None), initialize=3.17756893714648) m.x66 = Var(within=Reals, bounds=(None,None), initialize=0.647103242058539) m.x67 = Var(within=Reals, bounds=(None,None), initialize=0) m.x68 = Var(within=Reals, bounds=(None,None), initialize=0) m.x69 = Var(within=Reals, bounds=(None,None), initialize=0.917954032114244) m.x70 = Var(within=Reals, bounds=(0.587786664902119,0.587786664902119) , initialize=0.587786664902119) m.x71 = Var(within=Reals, bounds=(None,None), initialize=-1.01757906523769) m.x72 = Var(within=Reals, bounds=(None,None), initialize=-0.603587003016137) m.x73 = Var(within=Reals, bounds=(None,None), initialize=3.17756893714648) m.x74 = Var(within=Reals, bounds=(None,None), initialize=0.259570495557779) m.x75 = Var(within=Reals, bounds=(None,None), initialize=0.538673267899503) m.x76 = Var(within=Reals, bounds=(None,None), initialize=0.594337876584314) m.x77 = Var(within=Reals, bounds=(None,None), initialize=0.165264958927998) m.x78 = Var(within=Reals, bounds=(0.587786664902119,0.587786664902119) , initialize=0.587786664902119) m.x79 = Var(within=Reals, bounds=(None,None), initialize=-0.938029907934514) m.x80 = Var(within=Reals, bounds=(None,None), initialize=-0.456539109459263) m.x81 = Var(within=Reals, bounds=(0,1.6094379124341) , initialize=1.38629436111989) m.x82 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x83 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x84 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x85 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x86 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x87 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x88 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x89 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x90 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x91 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x92 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x93 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x94 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x95 = Var(within=Reals, bounds=(0,1.6094379124341), initialize=0) m.x96 = Var(within=Reals, bounds=(0,1.6094379124341) , initialize=1.38629436111989) m.x97 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=3.68887945411394) m.x98 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=4.38202663467388) m.x99 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=2.0433602879393) m.x100 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=2.73650746849925) m.x101 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=1.9677597942808) m.x102 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=1.9677597942808) m.x103 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=0.949618090789135) m.x104 = Var(within=Reals, bounds=(-2.30258509299405,6.90775527898214) , initialize=0) m.x105 = Var(within=Reals, bounds=(-6.90775527898214,2.99573227355399) , initialize=-1.15316440313985) m.x106 = Var(within=Reals, bounds=(-6.90775527898214,2.99573227355399) , initialize=-1.15316440313985) m.x107 = Var(within=Reals, bounds=(-6.90775527898214,2.99573227355399) , initialize=-0.824868956414142) m.x108 = Var(within=Reals, bounds=(-6.90775527898214,2.99573227355399) , initialize=-1.15316440313985) m.b109 = Var(within=Binary, bounds=(0,1), initialize=0) m.b110 = Var(within=Binary, bounds=(0,1), initialize=0) m.b111 = Var(within=Binary, bounds=(0,1), initialize=0) m.b112 = Var(within=Binary, bounds=(0,1), initialize=0) m.b113 = Var(within=Binary, bounds=(0,1), initialize=0) m.b114 = Var(within=Binary, bounds=(0,1), initialize=0) m.b115 = Var(within=Binary, bounds=(0,1), initialize=0) m.b116 = Var(within=Binary, bounds=(0,0), initialize=0) m.b117 = Var(within=Binary, bounds=(0,1), initialize=0) m.b118 = Var(within=Binary, bounds=(0,1), initialize=0) m.b119 = Var(within=Binary, bounds=(0,1), initialize=0) m.b120 = Var(within=Binary, bounds=(0,1), initialize=0) m.b121 = Var(within=Binary, bounds=(0,1), initialize=0) m.b122 = Var(within=Binary, bounds=(0,1), initialize=0) m.b123 = Var(within=Binary, bounds=(0,1), initialize=0) m.b124 = Var(within=Binary, bounds=(0,1), initialize=0) m.b125 = Var(within=Binary, bounds=(0,1), initialize=0) m.b126 = Var(within=Binary, bounds=(0,1), initialize=0) m.b127 = Var(within=Binary, bounds=(0,1), initialize=0) m.b128 = Var(within=Binary, bounds=(0,1), initialize=0) m.b129 = Var(within=Binary, bounds=(0,1), initialize=0) m.b130 = Var(within=Binary, bounds=(0,1), initialize=0) m.b131 = Var(within=Binary, bounds=(0,1), initialize=0) m.b132 = Var(within=Binary, bounds=(0,1), initialize=0) m.b133 = Var(within=Binary, bounds=(0,1), initialize=0) m.b134 = Var(within=Binary, bounds=(0,1), initialize=0) m.b135 = Var(within=Binary, bounds=(0,1), initialize=0) m.b136 = Var(within=Binary, bounds=(0,1), initialize=0) m.b137 = Var(within=Binary, bounds=(0,1), initialize=0) m.b138 = Var(within=Binary, bounds=(0,1), initialize=0) m.b139 = Var(within=Binary, bounds=(0,1), initialize=0) m.b140 = Var(within=Binary, bounds=(0,1), initialize=0) m.b141 = Var(within=Binary, bounds=(0,1), initialize=0) m.b142 = Var(within=Binary, bounds=(0,1), initialize=0) m.b143 = Var(within=Binary, bounds=(0,1), initialize=0) m.b144 = Var(within=Binary, bounds=(0,1), initialize=0) m.b145 = Var(within=Binary, bounds=(0,1), initialize=0) m.b146 = Var(within=Binary, bounds=(0,1), initialize=0) m.b147 = Var(within=Binary, bounds=(0,1), initialize=0) m.b148 = Var(within=Binary, bounds=(0,1), initialize=0) m.b149 = Var(within=Binary, bounds=(0,1), initialize=0) m.b150 = Var(within=Binary, bounds=(0,1), initialize=0) m.b151 = Var(within=Binary, bounds=(0,1), initialize=0) m.b152 = Var(within=Binary, bounds=(0,1), initialize=0) m.b153 = Var(within=Binary, bounds=(0,1), initialize=0) m.b154 = Var(within=Binary, bounds=(0,1), initialize=0) m.b155 = Var(within=Binary, bounds=(0,1), initialize=0) m.b156 = Var(within=Binary, bounds=(0,1), initialize=0) m.b157 = Var(within=Binary, bounds=(0,1), initialize=0) m.b158 = Var(within=Binary, bounds=(0,1), initialize=0) m.b159 = Var(within=Binary, bounds=(0,1), initialize=0) m.b160 = Var(within=Binary, bounds=(0,1), initialize=0) m.b161 = Var(within=Binary, bounds=(0,1), initialize=0) m.b162 = Var(within=Binary, bounds=(0,1), initialize=0) m.b163 = Var(within=Binary, bounds=(0,1), initialize=0) m.b164 = Var(within=Binary, bounds=(0,1), initialize=0) m.b165 = Var(within=Binary, bounds=(0,1), initialize=0) m.b166 = Var(within=Binary, bounds=(0,1), initialize=0) m.b167 = Var(within=Binary, bounds=(0,1), initialize=0) m.b168 = Var(within=Binary, bounds=(0,1), initialize=0) m.b169 = Var(within=Binary, bounds=(0,1), initialize=0) m.b170 = Var(within=Binary, bounds=(0,1), initialize=0) m.b171 = Var(within=Binary, bounds=(0,1), initialize=0) m.b172 = Var(within=Binary, bounds=(0,1), initialize=0) m.b173 = Var(within=Binary, bounds=(0,1), initialize=0) m.b174 = Var(within=Binary, bounds=(0,1), initialize=0) m.b175 = Var(within=Binary, bounds=(0,1), initialize=0) m.b176 = Var(within=Binary, bounds=(0,1), initialize=0) m.b177 = Var(within=Binary, bounds=(0,1), initialize=0) m.b178 = Var(within=Binary, bounds=(0,1), initialize=0) m.b179 = Var(within=Binary, bounds=(0,1), initialize=0) m.b180 = Var(within=Binary, bounds=(0,1), initialize=0) m.b181 = Var(within=Binary, bounds=(0,1), initialize=0) m.b182 = Var(within=Binary, bounds=(0,1), initialize=0) m.b183 = Var(within=Binary, bounds=(0,1), initialize=0) m.b184 = Var(within=Binary, bounds=(0,1), initialize=0) m.b185 = Var(within=Binary, bounds=(0,1), initialize=0) m.b186 = Var(within=Binary, bounds=(0,1), initialize=0) m.b187 = Var(within=Binary, bounds=(0,1), initialize=0) m.b188 = Var(within=Binary, bounds=(0,1), initialize=0) m.b189 = Var(within=Binary, bounds=(0,1), initialize=0) m.b190 = Var(within=Binary, bounds=(0,1), initialize=0) m.b191 = Var(within=Binary, bounds=(0,1), initialize=0) m.b192 = Var(within=Binary, bounds=(0,1), initialize=0) m.b193 = Var(within=Binary, bounds=(0,1), initialize=0) m.b194 = Var(within=Binary, bounds=(0,1), initialize=0) m.b195 = Var(within=Binary, bounds=(0,1), initialize=0) m.b196 = Var(within=Binary, bounds=(0,1), initialize=0) m.x197 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.862636840383206) m.x198 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.862636840383206) m.x199 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.2489312015031) m.x200 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.2489312015031) m.x201 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.55578402094315) m.x202 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.63138451460165) m.x203 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.32453169516159) m.x204 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-3.09889665499776) m.x205 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.752600105692548) m.x206 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.752600105692548) m.x207 = Var(within=Reals, bounds=(None,1.6094379124341), initialize=0) m.x208 = Var(within=Reals, bounds=(None,1.6094379124341), initialize=0) m.x209 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.0594529251326021) m.x210 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.17980108705855) m.x211 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.8729482676185) m.x212 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-3.39972974065073) m.x213 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=0.379588741113263) m.x214 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=0.379588741113263) m.x215 = Var(within=Reals, bounds=(None,1.6094379124341), initialize=0) m.x216 = Var(within=Reals, bounds=(None,1.6094379124341), initialize=0) m.x217 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=1.07273592167321) m.x218 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.04921352161972) m.x219 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.74236070217966) m.x220 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.76785883450285) m.x221 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.207931510667194) m.x222 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.207931510667194) m.x223 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.59422587178708) m.x224 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.59422587178708) m.x225 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.901078691227139) m.x226 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-0.976679184885637) m.x227 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-1.66982636544558) m.x228 = Var(within=Reals, bounds=(None,1.6094379124341) , initialize=-2.68796806893725) m.x229 = Var(within=Reals, bounds=(None,None), initialize=-1.55578402094315) m.x230 = Var(within=Reals, bounds=(-2.25129179860649,-2.25129179860649) , initialize=-2.25129179860649) m.x231 = Var(within=Reals, bounds=(None,None), initialize=-0.0594529251326021) m.x232 = Var(within=Reals, bounds=(-2.25129179860649,-2.25129179860649) , initialize=-2.25129179860649) m.x233 = Var(within=Reals, bounds=(None,None), initialize=1.07273592167321) m.x234 = Var(within=Reals, bounds=(-2.25129179860649,-2.25129179860649) , initialize=-2.25129179860649) m.x235 = Var(within=Reals, bounds=(None,None), initialize=-0.901078691227139) m.x236 = Var(within=Reals, bounds=(-2.25129179860649,-2.25129179860649) , initialize=-2.25129179860649) m.x237 = Var(within=Reals, bounds=(None,None), initialize=0.422047741265856) m.x238 = Var(within=Reals, bounds=(None,None), initialize=0.105511935316464) m.x239 = Var(within=Reals, bounds=(None,None), initialize=0.105511935316464) m.x240 = Var(within=Reals, bounds=(None,None), initialize=0.211023870632928) m.x241 = Var(within=Reals, bounds=(None,None), initialize=0.0978292472354046) m.x242 = Var(within=Reals, bounds=(None,None), initialize=0.0978292472354046) m.x243 = Var(within=Reals, bounds=(None,None), initialize=0.0450989346367396) m.x244 = Var(within=Reals, bounds=(None,None), initialize=0.471139945124032) m.x245 = Var(within=Reals, bounds=(None,None), initialize=0.942279890248065) m.x246 = Var(within=Reals, bounds=(0,0), initialize=0) m.x247 = Var(within=Reals, bounds=(0,0), initialize=0) m.x248 = Var(within=Reals, bounds=(None,None), initialize=0.0565320091501422) m.x249 = Var(within=Reals, bounds=(None,None), initialize=0.0565320091501422) m.x250 = Var(within=Reals, bounds=(None,None), initialize=0.033382290617452) m.x251 = Var(within=Reals, bounds=(None,None), initialize=1.46168333554586) m.x252 = Var(within=Reals, bounds=(None,None), initialize=2.92336667109171) m.x253 = Var(within=Reals, bounds=(0,0), initialize=0) m.x254 = Var(within=Reals, bounds=(0,0), initialize=0) m.x255 = Var(within=Reals, bounds=(None,None), initialize=0.175106537920172) m.x256 = Var(within=Reals, bounds=(None,None), initialize=0.175106537920172) m.x257 = Var(within=Reals, bounds=(None,None), initialize=0.0627963182070535) m.x258 = Var(within=Reals, bounds=(None,None), initialize=0.812262666134973) m.x259 = Var(within=Reals, bounds=(None,None), initialize=0.203065666533743) m.x260 = Var(within=Reals, bounds=(None,None), initialize=0.203065666533743) m.x261 = Var(within=Reals, bounds=(None,None), initialize=0.406131333067486) m.x262 = Var(within=Reals, bounds=(None,None), initialize=0.188279754672) m.x263 = Var(within=Reals, bounds=(None,None), initialize=0.188279754672) m.x264 = Var(within=Reals, bounds=(None,None), initialize=0.0680190089875055) m.x265 = Var(within=Reals, bounds=(0.1,1), initialize=0.903778326897858) m.x266 = Var(within=Reals, bounds=(0.1,1), initialize=0.975106465816068) m.x267 = Var(within=Reals, bounds=(0.1,1), initialize=0.882989317955773) m.x268 = Var(within=Reals, bounds=(0.95,0.95), initialize=0.95) m.x269 = Var(within=Reals, bounds=(0.1,1), initialize=0.998315513250229) m.x270 = Var(within=Reals, bounds=(0.1,1), initialize=0.931004925149313) m.x271 = Var(within=Reals, bounds=(0.95,0.95), initialize=0.95) m.x272 = Var(within=Reals, bounds=(0.1,1), initialize=0.998315513250229) m.x273 = Var(within=Reals, bounds=(0.1,1), initialize=0.901706690780938) m.x274 = Var(within=Reals, bounds=(0.95,0.95), initialize=0.95) m.x275 = Var(within=Reals, bounds=(0.1,1), initialize=0.903778326897858) m.x276 = Var(within=Reals, bounds=(0.1,1), initialize=0.975106465816068) m.x277 = Var(within=Reals, bounds=(0.1,1), initialize=0.917593932953673) m.x278 = Var(within=Reals, bounds=(0.95,0.95), initialize=0.95) m.x279 = Var(within=Reals, bounds=(10,54.5), initialize=50) m.x280 = Var(within=Reals, bounds=(10,54.5), initialize=50) m.x281 = Var(within=Reals, bounds=(10,54.5), initialize=50) m.x282 = Var(within=Reals, bounds=(10,54.5), initialize=50) m.x283 = Var(within=Reals, bounds=(10,250), initialize=200) m.x284 = Var(within=Reals, bounds=(10,250), initialize=200) m.x285 = Var(within=Reals, bounds=(10,250), initialize=200) m.x286 = Var(within=Reals, bounds=(10,250), initialize=200) m.x287 = Var(within=Reals, bounds=(0.1,10), initialize=1.25) m.x288 = Var(within=Reals, bounds=(0.1,10), initialize=1.25) m.x289 = Var(within=Reals, bounds=(0.5,10), initialize=1.5) m.x290 = Var(within=Reals, bounds=(0.5,10), initialize=1.5) m.x291 = Var(within=Reals, bounds=(0.5,6), initialize=3) m.x292 = Var(within=Reals, bounds=(0.5,6), initialize=3) m.x293 = Var(within=Reals, bounds=(0.1,4), initialize=1) m.x294 = Var(within=Reals, bounds=(0.1,4), initialize=1) m.x295 = Var(within=Reals, bounds=(0.1,4), initialize=1) m.x296 = Var(within=Reals, bounds=(0.1,4), initialize=1) m.obj = Objective(sense=minimize, expr= 63400 * exp(0.6 * m.x1 + m.x81 + m.x89) + 5750 * exp(0.6 * m.x2 + m.x82 + m.x90) + 5750 * exp(0.6 * m.x3 + m.x83 + m.x91) + 5750 * exp(0.6 * m.x4 + m.x84 + m.x92) + 5750 * exp(0.6 * m.x5 + m.x85 + m.x93) + 23100 * exp(0.65 * m.x6 + m.x86 + m.x94) + 5750 * exp( 0.6 * m.x7 + m.x87 + m.x95) + 5750 * exp(0.6 * m.x8 + m.x88 + m.x96) + 5750 * exp(0.6 * m.x9 + m.x82 + m.x90) + 5750 * exp(0.6 * m.x10 + m.x84 + m.x92) + 3.6e+05 * exp(0.975 * m.x11 + m.x88 + m.x96) + 2900 * exp(0.85 * m.x12 + m.x82 + m.x90) + 12100 * exp(0.75 * m.x13 + m.x83 + m.x91) + 2900 * exp( 0.85 * m.x14 + m.x84 + m.x92) + 2900 * exp(0.85 * m.x15 + m.x85 + m.x93) + 2900 * exp(0.85 * m.x16 + m.x87 + m.x95) + 5750 * (exp(0.6 * m.x97) + exp( 0.6 * m.x98) + exp(0.6 * m.x99) + exp(0.6 * m.x100) + exp(0.6 * m.x101) + exp(0.6 * m.x102) + exp(0.6 * m.x103) + exp(0.6 * m.x104))) m.e1 = Constraint(expr= m.x1 - m.x17 + m.x89 - m.x197 >= 0) m.e2 = Constraint(expr= m.x2 - m.x18 + m.x90 - m.x198 >= 0) m.e3 = Constraint(expr= m.x3 - m.x19 + m.x91 - m.x199 >= 0) m.e4 = Constraint(expr= m.x4 - m.x20 + m.x92 - m.x200 >= 0) m.e5 = Constraint(expr= m.x5 - m.x21 + m.x93 - m.x201 >= 0) m.e6 = Constraint(expr= m.x6 - m.x22 + m.x94 - m.x202 >= 0) m.e7 = Constraint(expr= m.x7 - m.x23 + m.x95 - m.x203 >= 0) m.e8 = Constraint(expr= m.x8 - m.x24 + m.x96 - m.x204 >= 0) m.e9 = Constraint(expr= m.x1 - m.x25 + m.x89 - m.x205 >= 0) m.e10 = Constraint(expr= m.x2 - m.x26 + m.x90 - m.x206 >= 0) m.e11 = Constraint(expr= m.x3 - m.x27 + m.x91 - m.x207 >= 0) m.e12 = Constraint(expr= m.x4 - m.x28 + m.x92 - m.x208 >= 0) m.e13 = Constraint(expr= m.x5 - m.x29 + m.x93 - m.x209 >= 0) m.e14 = Constraint(expr= m.x6 - m.x30 + m.x94 - m.x210 >= 0) m.e15 = Constraint(expr= m.x7 - m.x31 + m.x95 - m.x211 >= 0) m.e16 = Constraint(expr= m.x8 - m.x32 + m.x96 - m.x212 >= 0) m.e17 = Constraint(expr= m.x1 - m.x33 + m.x89 - m.x213 >= 0) m.e18 = Constraint(expr= m.x2 - m.x34 + m.x90 - m.x214 >= 0) m.e19 = Constraint(expr= m.x3 - m.x35 + m.x91 - m.x215 >= 0) m.e20 = Constraint(expr= m.x4 - m.x36 + m.x92 - m.x216 >= 0) m.e21 = Constraint(expr= m.x5 - m.x37 + m.x93 - m.x217 >= 0) m.e22 = Constraint(expr= m.x6 - m.x38 + m.x94 - m.x218 >= 0) m.e23 = Constraint(expr= m.x7 - m.x39 + m.x95 - m.x219 >= 0) m.e24 = Constraint(expr= m.x8 - m.x40 + m.x96 - m.x220 >= 0) m.e25 = Constraint(expr= m.x1 - m.x41 + m.x89 - m.x221 >= 0) m.e26 = Constraint(expr= m.x2 - m.x42 + m.x90 - m.x222 >= 0) m.e27 = Constraint(expr= m.x3 - m.x43 + m.x91 - m.x223 >= 0) m.e28 = Constraint(expr= m.x4 - m.x44 + m.x92 - m.x224 >= 0) m.e29 = Constraint(expr= m.x5 - m.x45 + m.x93 - m.x225 >= 0) m.e30 = Constraint(expr= m.x6 - m.x46 + m.x94 - m.x226 >= 0) m.e31 = Constraint(expr= m.x7 - m.x47 + m.x95 - m.x227 >= 0) m.e32 = Constraint(expr= m.x8 - m.x48 + m.x96 - m.x228 >= 0) m.e33 = Constraint(expr= m.x10 - m.x20 + m.x92 - m.x229 >= 0) m.e34 = Constraint(expr= m.x11 - m.x24 + m.x96 - m.x230 >= 0) m.e35 = Constraint(expr= m.x9 - m.x26 + m.x90 - m.x231 >= 0) m.e36 = Constraint(expr= m.x11 - m.x32 + m.x96 - m.x232 >= 0) m.e37 = Constraint(expr= m.x9 - m.x34 + m.x90 - m.x233 >= 0) m.e38 = Constraint(expr= m.x11 - m.x40 + m.x96 - m.x234 >= 0) m.e39 = Constraint(expr= m.x10 - m.x44 + m.x92 - m.x235 >= 0) m.e40 = Constraint(expr= m.x11 - m.x48 + m.x96 - m.x236 >= 0) m.e41 = Constraint(expr= -log(15.6 / (m.x279 * m.x265 * m.x266 * m.x267 * m.x268)) + m.x197 == 0) m.e42 = Constraint(expr= -log(15.6 / (m.x279 * m.x265 * m.x266 * m.x267 * m.x268)) + m.x198 == 0) m.e43 = Constraint(expr= -log(15.6 / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268)) + m.x199 == 0) m.e44 = Constraint(expr= -log(15.6 / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268)) + m.x200 == 0) m.e45 = Constraint(expr= -log((7.8 + 15.6 * m.x289) / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268)) + m.x201 == 0) m.e46 = Constraint(expr= -log((0.075 + 0.075 * m.x293) / (m.x268 * exp(-0.03 * m.x291) * m.x267)) + m.x202 == 0) m.e47 = Constraint(expr= -log(0.075 / (m.x268 * exp(-0.03 * m.x291) * m.x267)) + m.x203 == 0) m.e48 = Constraint(expr= -log((0.05 * m.x293 / (1 + m.x293)**2 + 0.025 * m.x267) / (m.x267 * m.x268 * exp(-0.03 * m.x291))) + m.x204 == 0) m.e49 = Constraint(expr= -log(20.8 / (m.x280 * m.x271 * m.x269 * m.x270)) + m.x205 == 0) m.e50 = Constraint(expr= -log(20.8 / (m.x280 * m.x271 * m.x269 * m.x270)) + m.x206 == 0) m.e51 = Constraint(expr= -log((-20.8 * m.x280 / m.x284 + 20.8 + 20.8 * m.x287) / (m.x280 * m.x269 * m.x270 * m.x271)) + m.x209 == 0) m.e52 = Constraint(expr= -log((0.05 + 0.05 * m.x294) / (m.x271 * m.x270)) + m.x210 == 0) m.e53 = Constraint(expr= -log(0.05 / (m.x271 * m.x270)) + m.x211 == 0) m.e54 = Constraint(expr= -log((0.025 * m.x294 / (1 + m.x294)**2 + 0.025 * m.x270) / (m.x270 * m.x271)) + m.x212 == 0) m.e55 = Constraint(expr= -log(62.5 / (m.x281 * m.x274 * m.x272 * m.x273)) + m.x213 == 0) m.e56 = Constraint(expr= -log(62.5 / (m.x281 * m.x274 * m.x272 * m.x273)) + m.x214 == 0) m.e57 = Constraint(expr= -log((-62.5 * m.x281 / m.x285 + 62.5 + 62.5 * m.x288) / (m.x281 * m.x272 * m.x273 * m.x274)) + m.x217 == 0) m.e58 = Constraint(expr= -log((0.15 + 0.15 * m.x295) / (m.x274 * m.x273)) + m.x218 == 0) m.e59 = Constraint(expr= -log(0.15 / (m.x274 * m.x273)) + m.x219 == 0) m.e60 = Constraint(expr= -log((0.125 * m.x295 / (1 + m.x295)**2 + 0.025 * m.x273) / (m.x273 * m.x274)) + m.x220 == 0) m.e61 = Constraint(expr= -log(31.2 / (m.x282 * m.x275 * m.x276 * m.x277 * m.x278)) + m.x221 == 0) m.e62 = Constraint(expr= -log(31.2 / (m.x282 * m.x275 * m.x276 * m.x277 * m.x278)) + m.x222 == 0) m.e63 = Constraint(expr= -log(31.2 / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278)) + m.x223 == 0) m.e64 = Constraint(expr= -log(31.2 / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278)) + m.x224 == 0) m.e65 = Constraint(expr= -log((15.6 + 31.2 * m.x290) / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278)) + m.x225 == 0) m.e66 = Constraint(expr= -log((0.15 + 0.15 * m.x296) / (m.x278 * exp(-0.03 * m.x292) * m.x277)) + m.x226 == 0) m.e67 = Constraint(expr= -log(0.15 / (m.x278 * exp(-0.03 * m.x292) * m.x277)) + m.x227 == 0) m.e68 = Constraint(expr= -log((0.125 * m.x296 / (1 + m.x296)**2 + 0.025 * m.x277) / (m.x277 * m.x278 * exp(-0.03 * m.x292))) + m.x228 == 0) m.e69 = Constraint(expr= -log((7.8 + 15.6 * m.x289) / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268)) + m.x229 == 0) m.e70 = Constraint(expr= -log((-20.8 * m.x280 / m.x284 + 20.8 + 20.8 * m.x287) / (m.x280 * m.x269 * m.x270 * m.x271)) + m.x231 == 0) m.e71 = Constraint(expr= -log((-62.5 * m.x281 / m.x285 + 62.5 + 62.5 * m.x288) / (m.x281 * m.x272 * m.x273 * m.x274)) + m.x233 == 0) m.e72 = Constraint(expr= -log((15.6 + 31.2 * m.x290) / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278)) + m.x235 == 0) m.e73 = Constraint(expr= -exp(m.x198) + m.x237 == 0) m.e74 = Constraint(expr= -15.6 / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268) + m.x238 == 0) m.e75 = Constraint(expr= -15.6 / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268) + m.x239 == 0) m.e76 = Constraint(expr= -(7.8 + 15.6 * m.x289) / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268) + m.x240 == 0) m.e77 = Constraint(expr= -0.075 / (m.x268 * exp(-0.03 * m.x291) * m.x267) + m.x241 == 0) m.e78 = Constraint(expr= -0.075 / (m.x268 * exp(-0.03 * m.x291) * m.x267) + m.x242 == 0) m.e79 = Constraint(expr= -exp(m.x204) + m.x243 == 0) m.e80 = Constraint(expr= -exp(m.x206) + m.x244 == 0) m.e81 = Constraint(expr= -(-20.8 * m.x280 / m.x284 + 20.8 + 20.8 * m.x287) / ( m.x280 * m.x269 * m.x270 * m.x271) + m.x245 == 0) m.e82 = Constraint(expr= -0.05 / (m.x271 * m.x270) + m.x248 == 0) m.e83 = Constraint(expr= -0.05 / (m.x271 * m.x270) + m.x249 == 0) m.e84 = Constraint(expr= -exp(m.x212) + m.x250 == 0) m.e85 = Constraint(expr= -exp(m.x214) + m.x251 == 0) m.e86 = Constraint(expr= -(-62.5 * m.x281 / m.x285 + 62.5 + 62.5 * m.x288) / ( m.x281 * m.x272 * m.x273 * m.x274) + m.x252 == 0) m.e87 = Constraint(expr= -0.15 / (m.x274 * m.x273) + m.x255 == 0) m.e88 = Constraint(expr= -0.15 / (m.x274 * m.x273) + m.x256 == 0) m.e89 = Constraint(expr= -exp(m.x220) + m.x257 == 0) m.e90 = Constraint(expr= -exp(m.x222) + m.x258 == 0) m.e91 = Constraint(expr= -31.2 / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278) + m.x259 == 0) m.e92 = Constraint(expr= -31.2 / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278) + m.x260 == 0) m.e93 = Constraint(expr= -(15.6 + 31.2 * m.x290) / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278) + m.x261 == 0) m.e94 = Constraint(expr= -0.15 / (m.x278 * exp(-0.03 * m.x292) * m.x277) + m.x262 == 0) m.e95 = Constraint(expr= -0.15 / (m.x278 * exp(-0.03 * m.x292) * m.x277) + m.x263 == 0) m.e96 = Constraint(expr= -exp(m.x228) + m.x264 == 0) m.e97 = Constraint(expr= -log(3.8 * log(0.35 * m.x279 / (1 - 0.0181818181818182 * m.x279)) + 4) + m.x49 == 0) m.e98 = Constraint(expr= -log(m.x265 * (-62.5 * m.x279 / m.x283 + 62.5) / m.x279 * m.x266 * m.x267 * m.x268 * exp(m.x18) / exp(m.x12) + 1.25) + m.x50 == 0) m.e99 = Constraint(expr= -log(12.5 * m.x291 / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268) * exp(m.x19) / exp(m.x13) + 1.25) + m.x51 == 0) m.e100 = Constraint(expr= -log((62.5 + 125 * m.x289) / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268) * exp(m.x20) / exp(m.x14) + 1.75) + m.x52 == 0) m.e101 = Constraint(expr= -log((312.5 + 625 * m.x289) / (m.x283 * m.x265 * m.x266 * m.x267 * m.x268) * (-m.x283 * (-0.24 * exp(-1.5 * m.x291) + 0.24) * (-0.5 * exp(-2 * m.x289) + 1) / (25 + 50 * m.x289) + 1) * exp(m.x21) / exp(m.x15) + 1) + m.x53 == 0) m.e102 = Constraint(expr= -log((-(2 * m.x293) / (1 + m.x293)**2 + 3 - m.x267) / (m.x267 * m.x268 * exp(-0.03 * m.x291)) * exp(m.x23) / exp(m.x16) + 0.3) + m.x55 == 0) m.e103 = Constraint(expr= -log((0.005 * m.x293 / (1 + m.x293)**2 + 0.0025 * m.x267) / (m.x267 * m.x268 * exp(-0.03 * m.x291)) * exp(m.x24) / exp(m.x11) + 0.375) + m.x56 == 0) m.e104 = Constraint(expr= -log(3.8 * log(0.35 * m.x280 / (1 - 0.0181818181818182 * m.x280)) + 4) + m.x57 == 0) m.e105 = Constraint(expr= -log(m.x269 * (-83.5 * m.x280 / m.x284 + 83.5 + 83.5 * m.x287) / m.x280 * m.x270 * m.x271 * exp(m.x26) / exp(m.x12) + 1.75) + m.x58 == 0) m.e106 = Constraint(expr= -log((-835 * m.x280 / m.x284 + 835 + 835 * m.x287) / (m.x280 * m.x269 * m.x270 * m.x271) * (-0.12 * m.x280 * m.x269 / (-50 * m.x280 / m.x284 + 50 + 50 * m.x287) + 1) * exp(m.x29) / exp(m.x15) + 1) + m.x61 == 0) m.e107 = Constraint(expr= -log((-m.x294 / (1 + m.x294)**2 + 2 - m.x270) / ( m.x270 * m.x271) * exp(m.x31) / exp(m.x16) + 0.3) + m.x63 == 0) m.e108 = Constraint(expr= -log((0.0025 * m.x294 / (1 + m.x294)**2 + 0.0025 * m.x270) / (m.x270 * m.x271) * exp(m.x32) / exp(m.x11) + 0.375) + m.x64 == 0) m.e109 = Constraint(expr= -log(3.8 * log(0.35 * m.x281 / (1 - 0.0181818181818182 * m.x281)) + 4) + m.x65 == 0) m.e110 = Constraint(expr= -log((-250 * m.x281 / m.x285 + 250 + 250 * m.x288) / (m.x281 * m.x272 * m.x273 * m.x274) * exp(m.x34) / exp(m.x12) + 1.75) + m.x66 == 0) m.e111 = Constraint(expr= -log((-2500 * m.x281 / m.x285 + 2500 + 2500 * m.x288) / (m.x281 * m.x272 * m.x273 * m.x274) * (-0.12 * m.x281 * m.x272 / (-50 * m.x281 / m.x285 + 50 + 50 * m.x288) + 1) * exp(m.x37) / exp(m.x15) + 1) + m.x69 == 0) m.e112 = Constraint(expr= -log((-(5 * m.x295) / (1 + m.x295)**2 + 6 - m.x273) / (m.x273 * m.x274) * exp(m.x39) / exp(m.x16) + 0.3) + m.x71 == 0) m.e113 = Constraint(expr= -log((0.0125 * m.x295 / (1 + m.x295)**2 + 0.0025 * m.x273) / (m.x273 * m.x274) * exp(m.x40) / exp(m.x11) + 0.375) + m.x72 == 0) m.e114 = Constraint(expr= -log(3.8 * log(0.35 * m.x282 / (1 - 0.0181818181818182 * m.x282)) + 4) + m.x73 == 0) m.e115 = Constraint(expr= -log((-125 * m.x282 / m.x286 + 125) / (m.x282 * m.x275 * m.x276 * m.x277 * m.x278) * exp(m.x42) / exp(m.x12) + 1.25) + m.x74 == 0) m.e116 = Constraint(expr= -log(25 * m.x292 / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278) * exp(m.x43) / exp(m.x13) + 1.25) + m.x75 == 0) m.e117 = Constraint(expr= -log((125 + 250 * m.x290) / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278) * exp(m.x44) / exp(m.x14) + 1.75) + m.x76 == 0) m.e118 = Constraint(expr= -log((625 + 1250 * m.x290) / (m.x286 * m.x275 * m.x276 * m.x277 * m.x278) * (-m.x286 * (-0.24 * exp(-1.5 * m.x292) + 0.24) * (-0.5 * exp(-2 * m.x290) + 1) / (25 + 50 * m.x290) + 1) * exp(m.x45) / exp(m.x15) + 1) + m.x77 == 0) m.e119 = Constraint(expr= -log((-(5 * m.x296) / (1 + m.x296)**2 + 6 - m.x277) / (m.x277 * m.x278 * exp(-0.03 * m.x292)) * exp(m.x47) / exp(m.x16) + 0.3) + m.x79 == 0) m.e120 = Constraint(expr= -log((0.0125 * m.x296 / (1 + m.x296)**2 + 0.0025 * m.x277) / (m.x277 * m.x278 * exp(-0.03 * m.x292)) * exp(m.x48) / exp(m.x11) + 0.375) + m.x80 == 0) m.e121 = Constraint(expr= -(-exp(-1.5 * m.x291) + 1) * exp(-0.03 * m.x291) + m.x265 == 0) m.e122 = Constraint(expr= 0.5 * exp(-2 * m.x289) + m.x266 == 1) m.e123 = Constraint(expr= -(25.1 * m.x293) / ((m.x293 * 10**(4.9 * m.x293 / (1 + m.x293) - 3.5) + 1) * (1 + 25.1 * m.x293)) + m.x267 == 0) m.e124 = Constraint(expr= m.x280 / m.x284 * exp(-m.x287 * m.x284 / m.x280) + m.x269 == 1) m.e125 = Constraint(expr= -(50.1 * m.x294) / ((m.x294 * 10**(5.95 * m.x294 / (1 + m.x294) - 4.25) + 1) * (1 + 50.1 * m.x294)) + m.x270 == 0) m.e126 = Constraint(expr= m.x281 / m.x285 * exp(-m.x288 * m.x285 / m.x281) + m.x272 == 1) m.e127 = Constraint(expr= -(31.6 * m.x295) / ((m.x295 * 10**(5.25 * m.x295 / (1 + m.x295) - 3.75) + 1) * (1 + 31.6 * m.x295)) + m.x273 == 0) m.e128 = Constraint(expr= -(-exp(-1.5 * m.x292) + 1) * exp(-0.03 * m.x292) + m.x275 == 0) m.e129 = Constraint(expr= 0.5 * exp(-2 * m.x290) + m.x276 == 1) m.e130 = Constraint(expr= -(39.8 * m.x296) / ((m.x296 * 10**(5.6 * m.x296 / (1 + m.x296) - 4) + 1) * (1 + 39.8 * m.x296)) + m.x277 == 0) m.e131 = Constraint(expr= -m.x17 + m.x49 - m.x81 - m.x105 <= 0) m.e132 = Constraint(expr= -m.x18 + m.x50 - m.x82 - m.x105 <= 0) m.e133 = Constraint(expr= -m.x19 + m.x51 - m.x83 - m.x105 <= 0) m.e134 = Constraint(expr= -m.x20 + m.x52 - m.x84 - m.x105 <= 0) m.e135 = Constraint(expr= -m.x21 + m.x53 - m.x85 - m.x105 <= 0) m.e136 = Constraint(expr= -m.x22 + m.x54 - m.x86 - m.x105 <= 0) m.e137 = Constraint(expr= -m.x23 + m.x55 - m.x87 - m.x105 <= 0) m.e138 = Constraint(expr= -m.x24 + m.x56 - m.x88 - m.x105 <= 0) m.e139 = Constraint(expr= -m.x25 + m.x57 - m.x81 - m.x106 <= 0) m.e140 = Constraint(expr= -m.x26 + m.x58 - m.x82 - m.x106 <= 0) m.e141 = Constraint(expr= -m.x27 + m.x59 - m.x83 - m.x106 <= 0) m.e142 = Constraint(expr= -m.x28 + m.x60 - m.x84 - m.x106 <= 0) m.e143 = Constraint(expr= -m.x29 + m.x61 - m.x85 - m.x106 <= 0) m.e144 = Constraint(expr= -m.x30 + m.x62 - m.x86 - m.x106 <= 0) m.e145 = Constraint(expr= -m.x31 + m.x63 - m.x87 - m.x106 <= 0) m.e146 = Constraint(expr= -m.x32 + m.x64 - m.x88 - m.x106 <= 0) m.e147 = Constraint(expr= -m.x33 + m.x65 - m.x81 - m.x107 <= 0) m.e148 = Constraint(expr= -m.x34 + m.x66 - m.x82 - m.x107 <= 0) m.e149 = Constraint(expr= -m.x35 + m.x67 - m.x83 - m.x107 <= 0) m.e150 = Constraint(expr= -m.x36 + m.x68 - m.x84 - m.x107 <= 0) m.e151 = Constraint(expr= -m.x37 + m.x69 - m.x85 - m.x107 <= 0) m.e152 = Constraint(expr= -m.x38 + m.x70 - m.x86 - m.x107 <= 0) m.e153 = Constraint(expr= -m.x39 + m.x71 - m.x87 - m.x107 <= 0) m.e154 = Constraint(expr= -m.x40 + m.x72 - m.x88 - m.x107 <= 0) m.e155 = Constraint(expr= -m.x41 + m.x73 - m.x81 - m.x108 <= 0) m.e156 = Constraint(expr= -m.x42 + m.x74 - m.x82 - m.x108 <= 0) m.e157 = Constraint(expr= -m.x43 + m.x75 - m.x83 - m.x108 <= 0) m.e158 = Constraint(expr= -m.x44 + m.x76 - m.x84 - m.x108 <= 0) m.e159 = Constraint(expr= -m.x45 + m.x77 - m.x85 - m.x108 <= 0) m.e160 = Constraint(expr= -m.x46 + m.x78 - m.x86 - m.x108 <= 0) m.e161 = Constraint(expr= -m.x47 + m.x79 - m.x87 - m.x108 <= 0) m.e162 = Constraint(expr= -m.x48 + m.x80 - m.x88 - m.x108 <= 0) m.e163 = Constraint(expr= 6000 * exp(m.x105) + 3000 * exp(m.x106) + 1500 * exp( m.x107) + 1000 * exp(m.x108) <= 6000) m.e164 = Constraint(expr= m.x17 - m.x18 - 2.99573227355399 * m.b109 <= 0) m.e165 = Constraint(expr= m.x18 - m.x19 - 2.99573227355399 * m.b110 <= 0) m.e166 = Constraint(expr= m.x19 - m.x20 - 2.99573227355399 * m.b111 <= 0) m.e167 = Constraint(expr= m.x20 - m.x21 - 2.99573227355399 * m.b112 <= 0) m.e168 = Constraint(expr= m.x21 - m.x22 - 2.99573227355399 * m.b113 <= 0) m.e169 = Constraint(expr= m.x22 - m.x23 - 2.99573227355399 * m.b114 <= 0) m.e170 = Constraint(expr= m.x23 - m.x24 - 2.99573227355399 * m.b115 <= 0) m.e171 = Constraint(expr= m.x25 - m.x26 - 2.99573227355399 * m.b109 <= 0) m.e172 = Constraint(expr= m.x26 - m.x27 - 2.99573227355399 * m.b110 <= 0) m.e173 = Constraint(expr= m.x27 - m.x28 - 2.99573227355399 * m.b111 <= 0) m.e174 = Constraint(expr= m.x28 - m.x29 - 2.99573227355399 * m.b112 <= 0) m.e175 = Constraint(expr= m.x29 - m.x30 - 2.99573227355399 * m.b113 <= 0) m.e176 = Constraint(expr= m.x30 - m.x31 - 2.99573227355399 * m.b114 <= 0) m.e177 = Constraint(expr= m.x31 - m.x32 - 2.99573227355399 * m.b115 <= 0) m.e178 = Constraint(expr= m.x33 - m.x34 - 2.99573227355399 * m.b109 <= 0) m.e179 = Constraint(expr= m.x34 - m.x35 - 2.99573227355399 * m.b110 <= 0) m.e180 = Constraint(expr= m.x35 - m.x36 - 2.99573227355399 * m.b111 <= 0) m.e181 = Constraint(expr= m.x36 - m.x37 - 2.99573227355399 * m.b112 <= 0) m.e182 = Constraint(expr= m.x37 - m.x38 - 2.99573227355399 * m.b113 <= 0) m.e183 = Constraint(expr= m.x38 - m.x39 - 2.99573227355399 * m.b114 <= 0) m.e184 = Constraint(expr= m.x39 - m.x40 - 2.99573227355399 * m.b115 <= 0) m.e185 = Constraint(expr= m.x41 - m.x42 - 2.99573227355399 * m.b109 <= 0) m.e186 = Constraint(expr= m.x42 - m.x43 - 2.99573227355399 * m.b110 <= 0) m.e187 = Constraint(expr= m.x43 - m.x44 - 2.99573227355399 * m.b111 <= 0) m.e188 = Constraint(expr= m.x44 - m.x45 - 2.99573227355399 * m.b112 <= 0) m.e189 = Constraint(expr= m.x45 - m.x46 - 2.99573227355399 * m.b113 <= 0) m.e190 = Constraint(expr= m.x46 - m.x47 - 2.99573227355399 * m.b114 <= 0) m.e191 = Constraint(expr= m.x47 - m.x48 - 2.99573227355399 * m.b115 <= 0) m.e192 = Constraint(expr= m.x17 - m.x18 + 2.99573227355399 * m.b109 >= 0) m.e193 = Constraint(expr= m.x18 - m.x19 + 2.99573227355399 * m.b110 >= 0) m.e194 = Constraint(expr= m.x19 - m.x20 + 2.99573227355399 * m.b111 >= 0) m.e195 = Constraint(expr= m.x20 - m.x21 + 2.99573227355399 * m.b112 >= 0) m.e196 = Constraint(expr= m.x21 - m.x22 + 2.99573227355399 * m.b113 >= 0) m.e197 = Constraint(expr= m.x22 - m.x23 + 2.99573227355399 * m.b114 >= 0) m.e198 = Constraint(expr= m.x23 - m.x24 + 2.99573227355399 * m.b115 >= 0) m.e199 = Constraint(expr= m.x25 - m.x26 + 2.99573227355399 * m.b109 >= 0) m.e200 = Constraint(expr= m.x26 - m.x27 + 2.99573227355399 * m.b110 >= 0) m.e201 = Constraint(expr= m.x27 - m.x28 + 2.99573227355399 * m.b111 >= 0) m.e202 = Constraint(expr= m.x28 - m.x29 + 2.99573227355399 * m.b112 >= 0) m.e203 = Constraint(expr= m.x29 - m.x30 + 2.99573227355399 * m.b113 >= 0) m.e204 = Constraint(expr= m.x30 - m.x31 + 2.99573227355399 * m.b114 >= 0) m.e205 = Constraint(expr= m.x31 - m.x32 + 2.99573227355399 * m.b115 >= 0) m.e206 = Constraint(expr= m.x33 - m.x34 + 2.99573227355399 * m.b109 >= 0) m.e207 = Constraint(expr= m.x34 - m.x35 + 2.99573227355399 * m.b110 >= 0) m.e208 = Constraint(expr= m.x35 - m.x36 + 2.99573227355399 * m.b111 >= 0) m.e209 = Constraint(expr= m.x36 - m.x37 + 2.99573227355399 * m.b112 >= 0) m.e210 = Constraint(expr= m.x37 - m.x38 + 2.99573227355399 * m.b113 >= 0) m.e211 = Constraint(expr= m.x38 - m.x39 + 2.99573227355399 * m.b114 >= 0) m.e212 = Constraint(expr= m.x39 - m.x40 + 2.99573227355399 * m.b115 >= 0) m.e213 = Constraint(expr= m.x41 - m.x42 + 2.99573227355399 * m.b109 >= 0) m.e214 = Constraint(expr= m.x42 - m.x43 + 2.99573227355399 * m.b110 >= 0) m.e215 = Constraint(expr= m.x43 - m.x44 + 2.99573227355399 * m.b111 >= 0) m.e216 = Constraint(expr= m.x44 - m.x45 + 2.99573227355399 * m.b112 >= 0) m.e217 = Constraint(expr= m.x45 - m.x46 + 2.99573227355399 * m.b113 >= 0) m.e218 = Constraint(expr= m.x46 - m.x47 + 2.99573227355399 * m.b114 >= 0) m.e219 = Constraint(expr= m.x47 - m.x48 + 2.99573227355399 * m.b115 >= 0) m.e220 = Constraint(expr= -m.x237 * (exp(m.x17 - m.x97) + exp(m.x18 - m.x97)) - 10 * m.b109 >= -11) m.e221 = Constraint(expr= -m.x238 * (exp(m.x18 - m.x98) + exp(m.x19 - m.x98)) - 10 * m.b110 >= -11) m.e222 = Constraint(expr= -m.x239 * (exp(m.x19 - m.x99) + exp(m.x20 - m.x99)) - 10 * m.b111 >= -11) m.e223 = Constraint(expr= -m.x240 * (exp(m.x20 - m.x100) + exp(m.x21 - m.x100)) - 10 * m.b112 >= -11) m.e224 = Constraint(expr= -m.x241 * (exp(m.x21 - m.x101) + exp(m.x22 - m.x101)) - 10 * m.b113 >= -11) m.e225 = Constraint(expr= -m.x242 * (exp(m.x22 - m.x102) + exp(m.x23 - m.x102)) - 10 * m.b114 >= -11) m.e226 = Constraint(expr= -m.x243 * (exp(m.x23 - m.x103) + exp(m.x24 - m.x103)) - 10 * m.b115 >= -11) m.e227 = Constraint(expr= -m.x244 * (exp(m.x25 - m.x97) + exp(m.x26 - m.x97)) - 10 * m.b109 >= -11) m.e228 = Constraint(expr= -m.x245 * (exp(m.x26 - m.x98) + exp(m.x27 - m.x98)) - 10 * m.b110 >= -11) m.e229 = Constraint(expr= -m.x246 * (exp(m.x27 - m.x99) + exp(m.x28 - m.x99)) - 10 * m.b111 >= -11) m.e230 = Constraint(expr= -m.x247 * (exp(m.x28 - m.x100) + exp(m.x29 - m.x100)) - 10 * m.b112 >= -11) m.e231 = Constraint(expr= -m.x248 * (exp(m.x29 - m.x101) + exp(m.x30 - m.x101)) - 10 * m.b113 >= -11) m.e232 = Constraint(expr= -m.x249 * (exp(m.x30 - m.x102) + exp(m.x31 - m.x102)) - 10 * m.b114 >= -11) m.e233 = Constraint(expr= -m.x250 * (exp(m.x31 - m.x103) + exp(m.x32 - m.x103)) - 10 * m.b115 >= -11) m.e234 = Constraint(expr= -m.x251 * (exp(m.x33 - m.x97) + exp(m.x34 - m.x97)) - 10 * m.b109 >= -11) m.e235 = Constraint(expr= -m.x252 * (exp(m.x34 - m.x98) + exp(m.x35 - m.x98)) - 10 * m.b110 >= -11) m.e236 = Constraint(expr= -m.x253 * (exp(m.x35 - m.x99) + exp(m.x36 - m.x99)) - 10 * m.b111 >= -11) m.e237 = Constraint(expr= -m.x254 * (exp(m.x36 - m.x100) + exp(m.x37 - m.x100)) - 10 * m.b112 >= -11) m.e238 = Constraint(expr= -m.x255 * (exp(m.x37 - m.x101) + exp(m.x38 - m.x101)) - 10 * m.b113 >= -11) m.e239 = Constraint(expr= -m.x256 * (exp(m.x38 - m.x102) + exp(m.x39 - m.x102)) - 10 * m.b114 >= -11) m.e240 = Constraint(expr= -m.x257 * (exp(m.x39 - m.x103) + exp(m.x40 - m.x103)) - 10 * m.b115 >= -11) m.e241 = Constraint(expr= -m.x258 * (exp(m.x41 - m.x97) + exp(m.x42 - m.x97)) - 10 * m.b109 >= -11) m.e242 = Constraint(expr= -m.x259 * (exp(m.x42 - m.x98) + exp(m.x43 - m.x98)) - 10 * m.b110 >= -11) m.e243 = Constraint(expr= -m.x260 * (exp(m.x43 - m.x99) + exp(m.x44 - m.x99)) - 10 * m.b111 >= -11) m.e244 = Constraint(expr= -m.x261 * (exp(m.x44 - m.x100) + exp(m.x45 - m.x100)) - 10 * m.b112 >= -11) m.e245 = Constraint(expr= -m.x262 * (exp(m.x45 - m.x101) + exp(m.x46 - m.x101)) - 10 * m.b113 >= -11) m.e246 = Constraint(expr= -m.x263 * (exp(m.x46 - m.x102) + exp(m.x47 - m.x102)) - 10 * m.b114 >= -11) m.e247 = Constraint(expr= -m.x264 * (exp(m.x47 - m.x103) + exp(m.x48 - m.x103)) - 10 * m.b115 >= -11) m.e248 = Constraint(expr= m.x81 - 0.693147180559945 * m.b125 - 1.09861228866811 * m.b133 - 1.38629436111989 * m.b141 - 1.6094379124341 * m.b149 == 0) m.e249 = Constraint(expr= m.x82 - 0.693147180559945 * m.b126 - 1.09861228866811 * m.b134 - 1.38629436111989 * m.b142 - 1.6094379124341 * m.b150 == 0) m.e250 = Constraint(expr= m.x83 - 0.693147180559945 * m.b127 - 1.09861228866811 * m.b135 - 1.38629436111989 * m.b143 - 1.6094379124341 * m.b151 == 0) m.e251 = Constraint(expr= m.x84 - 0.693147180559945 * m.b128 - 1.09861228866811 * m.b136 - 1.38629436111989 * m.b144 - 1.6094379124341 * m.b152 == 0) m.e252 = Constraint(expr= m.x85 - 0.693147180559945 * m.b129 - 1.09861228866811 * m.b137 - 1.38629436111989 * m.b145 - 1.6094379124341 * m.b153 == 0) m.e253 = Constraint(expr= m.x86 - 0.693147180559945 * m.b130 - 1.09861228866811 * m.b138 - 1.38629436111989 * m.b146 - 1.6094379124341 * m.b154 == 0) m.e254 = Constraint(expr= m.x87 - 0.693147180559945 * m.b131 - 1.09861228866811 * m.b139 - 1.38629436111989 * m.b147 - 1.6094379124341 * m.b155 == 0) m.e255 = Constraint(expr= m.x88 - 0.693147180559945 * m.b132 - 1.09861228866811 * m.b140 - 1.38629436111989 * m.b148 - 1.6094379124341 * m.b156 == 0) m.e256 = Constraint(expr= m.b117 + m.b125 + m.b133 + m.b141 + m.b149 == 1) m.e257 = Constraint(expr= m.b118 + m.b126 + m.b134 + m.b142 + m.b150 == 1) m.e258 = Constraint(expr= m.b119 + m.b127 + m.b135 + m.b143 + m.b151 == 1) m.e259 = Constraint(expr= m.b120 + m.b128 + m.b136 + m.b144 + m.b152 == 1) m.e260 = Constraint(expr= m.b121 + m.b129 + m.b137 + m.b145 + m.b153 == 1) m.e261 = Constraint(expr= m.b122 + m.b130 + m.b138 + m.b146 + m.b154 == 1) m.e262 = Constraint(expr= m.b123 + m.b131 + m.b139 + m.b147 + m.b155 == 1) m.e263 = Constraint(expr= m.b124 + m.b132 + m.b140 + m.b148 + m.b156 == 1) m.e264 = Constraint(expr= m.x89 - 0.693147180559945 * m.b165 - 1.09861228866811 * m.b173 - 1.38629436111989 * m.b181 - 1.6094379124341 * m.b189 == 0) m.e265 = Constraint(expr= m.x90 - 0.693147180559945 * m.b166 - 1.09861228866811 * m.b174 - 1.38629436111989 * m.b182 - 1.6094379124341 * m.b190 == 0) m.e266 = Constraint(expr= m.x91 - 0.693147180559945 * m.b167 - 1.09861228866811 * m.b175 - 1.38629436111989 * m.b183 - 1.6094379124341 * m.b191 == 0) m.e267 = Constraint(expr= m.x92 - 0.693147180559945 * m.b168 - 1.09861228866811 * m.b176 - 1.38629436111989 * m.b184 - 1.6094379124341 * m.b192 == 0) m.e268 = Constraint(expr= m.x93 - 0.693147180559945 * m.b169 - 1.09861228866811 * m.b177 - 1.38629436111989 * m.b185 - 1.6094379124341 * m.b193 == 0) m.e269 = Constraint(expr= m.x94 - 0.693147180559945 * m.b170 - 1.09861228866811 * m.b178 - 1.38629436111989 * m.b186 - 1.6094379124341 * m.b194 == 0) m.e270 = Constraint(expr= m.x95 - 0.693147180559945 * m.b171 - 1.09861228866811 * m.b179 - 1.38629436111989 * m.b187 - 1.6094379124341 * m.b195 == 0) m.e271 = Constraint(expr= m.x96 - 0.693147180559945 * m.b172 - 1.09861228866811 * m.b180 - 1.38629436111989 * m.b188 - 1.6094379124341 * m.b196 == 0) m.e272 = Constraint(expr= m.b157 + m.b165 + m.b173 + m.b181 + m.b189 == 1) m.e273 = Constraint(expr= m.b158 + m.b166 + m.b174 + m.b182 + m.b190 == 1) m.e274 = Constraint(expr= m.b159 + m.b167 + m.b175 + m.b183 + m.b191 == 1) m.e275 = Constraint(expr= m.b160 + m.b168 + m.b176 + m.b184 + m.b192 == 1) m.e276 = Constraint(expr= m.b161 + m.b169 + m.b177 + m.b185 + m.b193 == 1) m.e277 = Constraint(expr= m.b162 + m.b170 + m.b178 + m.b186 + m.b194 == 1) m.e278 = Constraint(expr= m.b163 + m.b171 + m.b179 + m.b187 + m.b195 == 1) m.e279 = Constraint(expr= m.b164 + m.b172 + m.b180 + m.b188 + m.b196 == 1) m.e280 = Constraint(expr= m.b109 + m.b110 + m.b111 + m.b112 + m.b113 + m.b114 + m.b115 + m.b116 <= 7) m.e281 = Constraint(expr= m.x279 - m.x283 <= 0) m.e282 = Constraint(expr= m.x280 - m.x284 <= 0) m.e283 = Constraint(expr= m.x281 - m.x285 <= 0) m.e284 = Constraint(expr= m.x282 - m.x286 <= 0)