# NLP written by GAMS Convert at 02/17/22 17:19:07
#
# Equation counts
#     Total        E        G        L        N        X        C        B
#         2        2        0        0        0        0        0        0
#
# Variable counts
#                  x        b        i      s1s      s2s       sc       si
#     Total     cont   binary  integer     sos1     sos2    scont     sint
#         4        4        0        0        0        0        0        0
# FX      0
#
# Nonzero counts
#     Total    const       NL
#         4        4        0
#
# Reformulation has removed 1 variable and 1 equation

from pyomo.environ import *

model = m = ConcreteModel()

m.x1 = Var(within=Reals, bounds=(1e-07,0.5), initialize=0.4998)
m.x2 = Var(within=Reals, bounds=(1e-07,0.5), initialize=0.0002)
m.x3 = Var(within=Reals, bounds=(1e-07,0.5), initialize=0.0451)
m.x4 = Var(within=Reals, bounds=(1e-07,0.5), initialize=0.4549)

m.obj = Objective(sense=minimize, expr= log(4.8274 * m.x1 + 0.92 * m.x3) * (
    31.4830434782609 * m.x1 + 6 * m.x3) - 1.36551138119385 * m.x1 +
    2.8555953099828 * m.x3 + 11.5030434782609 * log(m.x1 / (4.8274 * m.x1 +
    0.92 * m.x3)) * m.x1 + 20.98 * log(m.x1 / (4.196 * m.x1 + 1.4 * m.x3)) *
    m.x1 + 7 * log(m.x3 / (4.196 * m.x1 + 1.4 * m.x3)) * m.x3 + log(4.196 *
    m.x1 + 1.4 * m.x3) * (4.196 * m.x1 + 1.4 * m.x3) + 1.62 * log(m.x1 / (
    7.52678200680961 * m.x1 + 0.443737968424621 * m.x3)) * m.x1 + 0.848 * log(
    m.x1 / (7.52678200680961 * m.x1 + 0.443737968424621 * m.x3)) * m.x1 + 1.728
    * log(m.x1 / (1.82245052351472 * m.x1 + 1.4300083598626 * m.x3)) * m.x1 +
    1.4 * log(m.x3 / (0.504772348000588 * m.x1 + 1.4 * m.x3)) * m.x3 + log(
    4.8274 * m.x2 + 0.92 * m.x4) * (31.4830434782609 * m.x2 + 6 * m.x4) -
    1.36551138119385 * m.x2 + 2.8555953099828 * m.x4 + 11.5030434782609 * log(
    m.x2 / (4.8274 * m.x2 + 0.92 * m.x4)) * m.x2 + 20.98 * log(m.x2 / (4.196 *
    m.x2 + 1.4 * m.x4)) * m.x2 + 7 * log(m.x4 / (4.196 * m.x2 + 1.4 * m.x4)) *
    m.x4 + log(4.196 * m.x2 + 1.4 * m.x4) * (4.196 * m.x2 + 1.4 * m.x4) + 1.62
    * log(m.x2 / (7.52678200680961 * m.x2 + 0.443737968424621 * m.x4)) * m.x2
    + 0.848 * log(m.x2 / (7.52678200680961 * m.x2 + 0.443737968424621 * m.x4))
    * m.x2 + 1.728 * log(m.x2 / (1.82245052351472 * m.x2 + 1.4300083598626 *
    m.x4)) * m.x2 + 1.4 * log(m.x4 / (0.504772348000588 * m.x2 + 1.4 * m.x4))
    * m.x4 - 35.6790434782609 * log(m.x1) * m.x1 - 7.4 * log(m.x3) * m.x3 -
    35.6790434782609 * log(m.x2) * m.x2 - 7.4 * log(m.x4) * m.x4)

m.e1 = Constraint(expr= m.x1 + m.x2 == 0.5)
m.e2 = Constraint(expr= m.x3 + m.x4 == 0.5)