<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 02/17/22 17:22:59
#
# Equation counts
#     Total        E        G        L        N        X        C        B
#       113       11       36       66        0        0        0        0
#
# Variable counts
#                  x        b        i      s1s      s2s       sc       si
#     Total     cont   binary  integer     sos1     sos2    scont     sint
#        65       45       20        0        0        0        0        0
# FX      0
#
# Nonzero counts
#     Total    const       NL
#       281      267       14
#
# Reformulation has removed 1 variable and 1 equation

from pyomo.environ import *

model = m = ConcreteModel()

m.x1 = Var(within=Reals, bounds=(0,10), initialize=0)
m.x2 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x3 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x4 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x5 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x6 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x7 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x8 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x9 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x10 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x11 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x12 = Var(within=Reals, bounds=(0,7), initialize=0)
m.x13 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x14 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x15 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x16 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x17 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x18 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x19 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x20 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x21 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x22 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x23 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x24 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x25 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x26 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x27 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x28 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x29 = Var(within=Reals, bounds=(0,5), initialize=0)
m.x30 = Var(within=Reals, bounds=(0,5), initialize=0)
m.x31 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x32 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x33 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x34 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x35 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x36 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x37 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x38 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x39 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x40 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x41 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x42 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x43 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x44 = Var(within=Reals, bounds=(0,None), initialize=0)
m.x45 = Var(within=Reals, bounds=(0,None), initialize=0)
m.b46 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b47 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b48 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b49 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b50 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b51 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b52 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b53 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b54 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b55 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b56 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b57 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b58 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b59 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b60 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b61 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b62 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b63 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b64 = Var(within=Binary, bounds=(0,1), initialize=0)
m.b65 = Var(within=Binary, bounds=(0,1), initialize=0)

m.obj = Objective(sense=maximize, expr= 5 * m.x7 + 200 * m.x37 + 250 * m.x38 +
    200 * m.x39 + 700 * m.x40 + 400 * m.x41 + 500 * m.x42 + 400 * m.x43 + 600 *
    m.x44 + 700 * m.x45 - 5 * m.b46 - 8 * m.b47 - 6 * m.b48 - 10 * m.b49 - 6 *
    m.b50 - 7 * m.b51 - 4 * m.b52 - 5 * m.b53 - 2 * m.b54 - 4 * m.b55 - 3 *
    m.b56 - 7 * m.b57 - 3 * m.b58 - 2 * m.b59 - 4 * m.b60 - 2 * m.b61 - 3 *
    m.b62 - 5 * m.b63 - 2 * m.b64 - 8 * m.b65)

m.e1 = Constraint(expr= m.x1 - m.x2 - m.x3 == 0)
m.e2 = Constraint(expr= -m.x4 - m.x5 + m.x6 == 0)
m.e3 = Constraint(expr= m.x6 - m.x7 - m.x8 == 0)
m.e4 = Constraint(expr= m.x8 - m.x9 - m.x10 - m.x11 == 0)
m.e5 = Constraint(expr= m.x13 - m.x16 - m.x17 == 0)
m.e6 = Constraint(expr= m.x15 - m.x18 - m.x19 - m.x20 == 0)
m.e7 = Constraint(expr= m.x23 - m.x27 - m.x28 == 0)
m.e8 = Constraint(expr= -m.x24 - m.x30 + m.x31 == 0)
m.e9 = Constraint(expr= m.x25 - m.x32 - m.x33 == 0)
m.e10 = Constraint(expr= m.x26 - m.x34 - m.x35 - m.x36 == 0)
m.e11 = Constraint(expr= -log(1 + m.x2) + m.x4 + m.b46 &lt;= 1)
m.e12 = Constraint(expr= m.x2 - 10 * m.b46 &lt;= 0)
m.e13 = Constraint(expr= m.x4 - 2.39789527279837 * m.b46 &lt;= 0)
m.e14 = Constraint(expr= -1.2 * log(1 + m.x3) + m.x5 + m.b47 &lt;= 1)
m.e15 = Constraint(expr= m.x3 - 10 * m.b47 &lt;= 0)
m.e16 = Constraint(expr= m.x5 - 2.87747432735804 * m.b47 &lt;= 0)
m.e17 = Constraint(expr= -0.75 * m.x9 + m.x13 + m.b48 &lt;= 1)
m.e18 = Constraint(expr= -0.75 * m.x9 + m.x13 - m.b48 &gt;= -1)
m.e19 = Constraint(expr= m.x9 - 2.87747432735804 * m.b48 &lt;= 0)
m.e20 = Constraint(expr= m.x13 - 2.15810574551853 * m.b48 &lt;= 0)
m.e21 = Constraint(expr= -1.5 * log(1 + m.x10) + m.x14 + m.b49 &lt;= 1)
m.e22 = Constraint(expr= m.x10 - 2.87747432735804 * m.b49 &lt;= 0)
m.e23 = Constraint(expr= m.x14 - 2.03277599268042 * m.b49 &lt;= 0)
m.e24 = Constraint(expr= -m.x11 + m.x15 + m.b50 &lt;= 1)
m.e25 = Constraint(expr= -m.x11 + m.x15 - m.b50 &gt;= -1)
m.e26 = Constraint(expr= -0.5 * m.x12 + m.x15 + m.b50 &lt;= 1)
m.e27 = Constraint(expr= -0.5 * m.x12 + m.x15 - m.b50 &gt;= -1)
m.e28 = Constraint(expr= m.x11 - 2.87747432735804 * m.b50 &lt;= 0)
m.e29 = Constraint(expr= m.x12 - 7 * m.b50 &lt;= 0)
m.e30 = Constraint(expr= m.x15 - 3.5 * m.b50 &lt;= 0)
m.e31 = Constraint(expr= -1.25 * log(1 + m.x16) + m.x21 + m.b51 &lt;= 1)
m.e32 = Constraint(expr= m.x16 - 2.15810574551853 * m.b51 &lt;= 0)
m.e33 = Constraint(expr= m.x21 - 1.43746550029693 * m.b51 &lt;= 0)
m.e34 = Constraint(expr= -0.9 * log(1 + m.x17) + m.x22 + m.b52 &lt;= 1)
m.e35 = Constraint(expr= m.x17 - 2.15810574551853 * m.b52 &lt;= 0)
m.e36 = Constraint(expr= m.x22 - 1.03497516021379 * m.b52 &lt;= 0)
m.e37 = Constraint(expr= -log(1 + m.x14) + m.x23 + m.b53 &lt;= 1)
m.e38 = Constraint(expr= m.x14 - 2.03277599268042 * m.b53 &lt;= 0)
m.e39 = Constraint(expr= m.x23 - 1.10947836929589 * m.b53 &lt;= 0)
m.e40 = Constraint(expr= -0.9 * m.x18 + m.x24 + m.b54 &lt;= 1)
m.e41 = Constraint(expr= -0.9 * m.x18 + m.x24 - m.b54 &gt;= -1)
m.e42 = Constraint(expr= m.x18 - 3.5 * m.b54 &lt;= 0)
m.e43 = Constraint(expr= m.x24 - 3.15 * m.b54 &lt;= 0)
m.e44 = Constraint(expr= -0.6 * m.x19 + m.x25 + m.b55 &lt;= 1)
m.e45 = Constraint(expr= -0.6 * m.x19 + m.x25 - m.b55 &gt;= -1)
m.e46 = Constraint(expr= m.x19 - 3.5 * m.b55 &lt;= 0)
m.e47 = Constraint(expr= m.x25 - 2.1 * m.b55 &lt;= 0)
m.e48 = Constraint(expr= -1.1 * log(1 + m.x20) + m.x26 + m.b56 &lt;= 1)
m.e49 = Constraint(expr= m.x20 - 3.5 * m.b56 &lt;= 0)
m.e50 = Constraint(expr= m.x26 - 1.6544851364539 * m.b56 &lt;= 0)
m.e51 = Constraint(expr= -0.9 * m.x21 + m.x37 + m.b57 &lt;= 1)
m.e52 = Constraint(expr= -0.9 * m.x21 + m.x37 - m.b57 &gt;= -1)
m.e53 = Constraint(expr= -m.x29 + m.x37 + m.b57 &lt;= 1)
m.e54 = Constraint(expr= -m.x29 + m.x37 - m.b57 &gt;= -1)
m.e55 = Constraint(expr= m.x21 - 1.43746550029693 * m.b57 &lt;= 0)
m.e56 = Constraint(expr= m.x29 - 5 * m.b57 &lt;= 0)
m.e57 = Constraint(expr= m.x37 - 5 * m.b57 &lt;= 0)
m.e58 = Constraint(expr= -log(1 + m.x22) + m.x38 + m.b58 &lt;= 1)
m.e59 = Constraint(expr= m.x22 - 1.03497516021379 * m.b58 &lt;= 0)
m.e60 = Constraint(expr= m.x38 - 0.710483612536911 * m.b58 &lt;= 0)
m.e61 = Constraint(expr= -0.7 * log(1 + m.x27) + m.x39 + m.b59 &lt;= 1)
m.e62 = Constraint(expr= m.x27 - 1.10947836929589 * m.b59 &lt;= 0)
m.e63 = Constraint(expr= m.x39 - 0.522508489006913 * m.b59 &lt;= 0)
m.e64 = Constraint(expr= -0.65 * log(1 + m.x28) + m.x40 + m.b60 &lt;= 1)
m.e65 = Constraint(expr= -0.65 * log(1 + m.x31) + m.x40 + m.b60 &lt;= 1)
m.e66 = Constraint(expr= m.x28 - 1.10947836929589 * m.b60 &lt;= 0)
m.e67 = Constraint(expr= m.x31 - 8.15 * m.b60 &lt;= 0)
m.e68 = Constraint(expr= m.x40 - 1.43894002153683 * m.b60 &lt;= 0)
m.e69 = Constraint(expr= -m.x32 + m.x41 + m.b61 &lt;= 1)
m.e70 = Constraint(expr= -m.x32 + m.x41 - m.b61 &gt;= -1)
m.e71 = Constraint(expr= m.x32 - 2.1 * m.b61 &lt;= 0)
m.e72 = Constraint(expr= m.x41 - 2.1 * m.b61 &lt;= 0)
m.e73 = Constraint(expr= -m.x33 + m.x42 + m.b62 &lt;= 1)
m.e74 = Constraint(expr= -m.x33 + m.x42 - m.b62 &gt;= -1)
m.e75 = Constraint(expr= m.x33 - 2.1 * m.b62 &lt;= 0)
m.e76 = Constraint(expr= m.x42 - 2.1 * m.b62 &lt;= 0)
m.e77 = Constraint(expr= -0.75 * log(1 + m.x34) + m.x43 + m.b63 &lt;= 1)
m.e78 = Constraint(expr= m.x34 - 1.6544851364539 * m.b63 &lt;= 0)
m.e79 = Constraint(expr= m.x43 - 0.732188035236726 * m.b63 &lt;= 0)
m.e80 = Constraint(expr= -0.8 * log(1 + m.x35) + m.x44 + m.b64 &lt;= 1)
m.e81 = Constraint(expr= m.x35 - 1.6544851364539 * m.b64 &lt;= 0)
m.e82 = Constraint(expr= m.x44 - 0.781000570919175 * m.b64 &lt;= 0)
m.e83 = Constraint(expr= -0.85 * log(1 + m.x36) + m.x45 + m.b65 &lt;= 1)
m.e84 = Constraint(expr= m.x36 - 1.6544851364539 * m.b65 &lt;= 0)
m.e85 = Constraint(expr= m.x45 - 0.829813106601623 * m.b65 &lt;= 0)
m.e86 = Constraint(expr= m.b46 + m.b47 == 1)
m.e87 = Constraint(expr= -m.b48 + m.b51 + m.b52 &gt;= 0)
m.e88 = Constraint(expr= -m.b51 + m.b57 &gt;= 0)
m.e89 = Constraint(expr= -m.b52 + m.b58 &gt;= 0)
m.e90 = Constraint(expr= -m.b49 + m.b53 &gt;= 0)
m.e91 = Constraint(expr= -m.b53 + m.b59 + m.b60 &gt;= 0)
m.e92 = Constraint(expr= -m.b50 + m.b54 + m.b55 + m.b56 &gt;= 0)
m.e93 = Constraint(expr= -m.b54 + m.b60 &gt;= 0)
m.e94 = Constraint(expr= -m.b55 + m.b61 + m.b62 &gt;= 0)
m.e95 = Constraint(expr= -m.b56 + m.b63 + m.b64 + m.b65 &gt;= 0)
m.e96 = Constraint(expr= m.b46 + m.b47 - m.b48 &gt;= 0)
m.e97 = Constraint(expr= m.b46 + m.b47 - m.b49 &gt;= 0)
m.e98 = Constraint(expr= m.b46 + m.b47 - m.b50 &gt;= 0)
m.e99 = Constraint(expr= m.b48 - m.b51 &gt;= 0)
m.e100 = Constraint(expr= m.b48 - m.b52 &gt;= 0)
m.e101 = Constraint(expr= m.b49 - m.b53 &gt;= 0)
m.e102 = Constraint(expr= m.b50 - m.b54 &gt;= 0)
m.e103 = Constraint(expr= m.b50 - m.b55 &gt;= 0)
m.e104 = Constraint(expr= m.b50 - m.b56 &gt;= 0)
m.e105 = Constraint(expr= m.b51 - m.b57 &gt;= 0)
m.e106 = Constraint(expr= m.b52 - m.b58 &gt;= 0)
m.e107 = Constraint(expr= m.b53 - m.b59 &gt;= 0)
m.e108 = Constraint(expr= m.b53 - m.b60 &gt;= 0)
m.e109 = Constraint(expr= m.b55 - m.b61 &gt;= 0)
m.e110 = Constraint(expr= m.b55 - m.b62 &gt;= 0)
m.e111 = Constraint(expr= m.b56 - m.b63 &gt;= 0)
m.e112 = Constraint(expr= m.b56 - m.b64 &gt;= 0)
m.e113 = Constraint(expr= m.b56 - m.b65 &gt;= 0)
</pre></body></html>