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

from pyomo.environ import *

model = m = ConcreteModel()

m.x1 = Var(within=Reals, bounds=(1e-07,40.30707), initialize=31.459)
m.x2 = Var(within=Reals, bounds=(1e-07,40.30707), initialize=0.901998)
m.x3 = Var(within=Reals, bounds=(1e-07,40.30707), initialize=1e-07)
m.x4 = Var(within=Reals, bounds=(1e-07,5.14979), initialize=3.10348)
m.x5 = Var(within=Reals, bounds=(1e-07,5.14979), initialize=9.6e-06)
m.x6 = Var(within=Reals, bounds=(1e-07,5.14979), initialize=1e-07)
m.x7 = Var(within=Reals, bounds=(1e-07,54.54314), initialize=26.1669)
m.x8 = Var(within=Reals, bounds=(1e-07,54.54314), initialize=15.0141)
m.x9 = Var(within=Reals, bounds=(1e-07,54.54314), initialize=1e-07)

m.obj = Objective(sense=minimize, expr= m.x3 * (log(m.x3 / (m.x3 + m.x6 + m.x9))
    + 0.156969560191053) + m.x6 * (log(m.x6 / (m.x3 + m.x6 + m.x9)) +
    0.156969560191053) + m.x9 * (log(m.x9 / (m.x3 + m.x6 + m.x9)) +
    0.156969560191053) + log(3.9235 * m.x1 + 6.0909 * m.x4 + 0.92 * m.x7) * (
    26.9071667605344 * m.x1 + 41.7710875549227 * m.x4 + 6.30931398488382 * m.x7)
    + 0.113370955614741 * m.x1 - 2.43897680885457 * m.x4 + 2.8555953099828 *
    m.x7 + 9.58716676053442 * log(m.x1) * m.x1 + 16.9310875549227 * log(m.x4)
    * m.x4 + 0.309313984883821 * log(m.x7) * m.x7 - 9.58716676053442 * log(
    3.9235 * m.x1 + 6.0909 * m.x4 + 0.92 * m.x7) * m.x1 - 16.9310875549227 *
    log(3.9235 * m.x1 + 6.0909 * m.x4 + 0.92 * m.x7) * m.x4 - 0.309313984883821
    * log(3.9235 * m.x1 + 6.0909 * m.x4 + 0.92 * m.x7) * m.x7 + 18.32 * log(
    m.x1) * m.x1 + 25.84 * log(m.x4) * m.x4 + 7 * log(m.x7) * m.x7 - 18.32 *
    log(3.664 * m.x1 + 5.168 * m.x4 + 1.4 * m.x7) * m.x1 - 25.84 * log(3.664 *
    m.x1 + 5.168 * m.x4 + 1.4 * m.x7) * m.x4 - 7 * log(3.664 * m.x1 + 5.168 *
    m.x4 + 1.4 * m.x7) * m.x7 + log(4.0643 * m.x1 + 5.7409 * m.x4 + 1.6741 *
    m.x7) * (4.0643 * m.x1 + 5.7409 * m.x4 + 1.6741 * m.x7) + 4.0643 * log(m.x1)
    * m.x1 + 5.7409 * log(m.x4) * m.x4 + 1.6741 * log(m.x7) * m.x7 - 4.0643 *
    log(4.0643 * m.x1 + 3.22644664511275 * m.x4 + 1.44980651607875 * m.x7) *
    m.x1 - 5.7409 * log(5.31147575751424 * m.x1 + 5.7409 * m.x4 +
    0.00729924451284409 * m.x7) * m.x4 - 1.6741 * log(2.25846661774355 * m.x1
    + 3.70876916588753 * m.x4 + 1.6741 * m.x7) * m.x7 + log(3.9235 * m.x2 +
    6.0909 * m.x5 + 0.92 * m.x8) * (26.9071667605344 * m.x2 + 41.7710875549227
    * m.x5 + 6.30931398488382 * m.x8) + 0.113370955614741 * m.x2 -
    2.43897680885457 * m.x5 + 2.8555953099828 * m.x8 + 9.58716676053442 * log(
    m.x2) * m.x2 + 16.9310875549227 * log(m.x5) * m.x5 + 0.309313984883821 *
    log(m.x8) * m.x8 - 9.58716676053442 * log(3.9235 * m.x2 + 6.0909 * m.x5 +
    0.92 * m.x8) * m.x2 - 16.9310875549227 * log(3.9235 * m.x2 + 6.0909 * m.x5
    + 0.92 * m.x8) * m.x5 - 0.309313984883821 * log(3.9235 * m.x2 + 6.0909 *
    m.x5 + 0.92 * m.x8) * m.x8 + 18.32 * log(m.x2) * m.x2 + 25.84 * log(m.x5)
    * m.x5 + 7 * log(m.x8) * m.x8 - 18.32 * log(3.664 * m.x2 + 5.168 * m.x5 +
    1.4 * m.x8) * m.x2 - 25.84 * log(3.664 * m.x2 + 5.168 * m.x5 + 1.4 * m.x8)
    * m.x5 - 7 * log(3.664 * m.x2 + 5.168 * m.x5 + 1.4 * m.x8) * m.x8 + log(
    4.0643 * m.x2 + 5.7409 * m.x5 + 1.6741 * m.x8) * (4.0643 * m.x2 + 5.7409 *
    m.x5 + 1.6741 * m.x8) + 4.0643 * log(m.x2) * m.x2 + 5.7409 * log(m.x5) *
    m.x5 + 1.6741 * log(m.x8) * m.x8 - 4.0643 * log(4.0643 * m.x2 +
    3.22644664511275 * m.x5 + 1.44980651607875 * m.x8) * m.x2 - 5.7409 * log(
    5.31147575751424 * m.x2 + 5.7409 * m.x5 + 0.00729924451284409 * m.x8) *
    m.x5 - 1.6741 * log(2.25846661774355 * m.x2 + 3.70876916588753 * m.x5 +
    1.6741 * m.x8) * m.x8 - 0.3658348 * m.x1 - 0.3658348 * m.x2 - 0.9825555 *
    m.x4 - 0.9825555 * m.x5 - 0.3663657 * m.x7 - 0.3663657 * m.x8 -
    30.9714667605344 * log(m.x1) * m.x1 - 47.5119875549227 * log(m.x4) * m.x4
    - 7.98341398488382 * log(m.x7) * m.x7 - 30.9714667605344 * log(m.x2) *
    m.x2 - 47.5119875549227 * log(m.x5) * m.x5 - 7.98341398488382 * log(m.x8)
    * m.x8)

m.e1 = Constraint(expr= m.x1 + m.x2 + m.x3 == 40.30707)
m.e2 = Constraint(expr= m.x4 + m.x5 + m.x6 == 5.14979)
m.e3 = Constraint(expr= m.x7 + m.x8 + m.x9 == 54.54314)