# NLP written by GAMS Convert at 02/17/22 17:20:13
#
# Equation counts
# Total E G L N X C B
# 60 2 36 22 0 0 0 0
#
# Variable counts
# x b i s1s s2s sc si
# Total cont binary integer sos1 sos2 scont sint
# 22 22 0 0 0 0 0 0
# FX 1
#
# Nonzero counts
# Total const NL
# 191 45 146
#
# Reformulation has removed 1 variable and 1 equation
from pyomo.environ import *
model = m = ConcreteModel()
m.x1 = Var(within=Reals, bounds=(4,72), initialize=4)
m.x2 = Var(within=Reals, bounds=(1.2,16.8), initialize=1.2)
m.x3 = Var(within=Reals, bounds=(1.2,2.8), initialize=1.2)
m.x4 = Var(within=Reals, bounds=(0.6,17.4), initialize=0.6)
m.x5 = Var(within=Reals, bounds=(0.6,3.4), initialize=0.6)
m.x6 = Var(within=Reals, bounds=(0.8,17.2), initialize=0.8)
m.x7 = Var(within=Reals, bounds=(0.8,3.2), initialize=0.8)
m.x8 = Var(within=Reals, bounds=(1.7,16.3), initialize=1.7)
m.x9 = Var(within=Reals, bounds=(1.7,2.3), initialize=1.7)
m.x10 = Var(within=Reals, bounds=(0.5,17.5), initialize=0.5)
m.x11 = Var(within=Reals, bounds=(0.5,3.5), initialize=0.5)
m.x12 = Var(within=Reals, bounds=(1.3,16.7), initialize=1.3)
m.x13 = Var(within=Reals, bounds=(1.3,2.7), initialize=1.3)
m.x14 = Var(within=Reals, bounds=(2,16), initialize=2)
m.x15 = Var(within=Reals, bounds=(2,2), initialize=2)
m.x16 = Var(within=Reals, bounds=(1.3,16.7), initialize=1.3)
m.x17 = Var(within=Reals, bounds=(1.3,2.7), initialize=1.3)
m.x18 = Var(within=Reals, bounds=(0.6,17.4), initialize=0.6)
m.x19 = Var(within=Reals, bounds=(0.6,3.4), initialize=0.6)
m.x20 = Var(within=Reals, bounds=(0,18), initialize=0)
m.x21 = Var(within=Reals, bounds=(0,4), initialize=0)
m.x22 = Var(within=Reals, bounds=(0,72), initialize=0)
m.obj = Objective(sense=minimize, expr= m.x22)
m.e1 = Constraint(expr= -m.x1 + m.x22 == -41.846014145816)
m.e2 = Constraint(expr= -m.x20 * m.x21 + m.x1 == 0)
m.e3 = Constraint(expr= (m.x2 - m.x4) * (m.x2 - m.x4) + (m.x3 - m.x5) * (m.x3
- m.x5) >= 3.24)
m.e4 = Constraint(expr= (m.x2 - m.x6) * (m.x2 - m.x6) + (m.x3 - m.x7) * (m.x3
- m.x7) >= 4)
m.e5 = Constraint(expr= (m.x2 - m.x8) * (m.x2 - m.x8) + (m.x3 - m.x9) * (m.x3
- m.x9) >= 8.41)
m.e6 = Constraint(expr= (m.x2 - m.x10) * (m.x2 - m.x10) + (m.x3 - m.x11) * (
m.x3 - m.x11) >= 2.89)
m.e7 = Constraint(expr= (m.x2 - m.x12) * (m.x2 - m.x12) + (m.x3 - m.x13) * (
m.x3 - m.x13) >= 6.25)
m.e8 = Constraint(expr= (m.x2 - m.x14) * (m.x2 - m.x14) + (m.x3 - m.x15) * (
m.x3 - m.x15) >= 10.24)
m.e9 = Constraint(expr= (m.x2 - m.x16) * (m.x2 - m.x16) + (m.x3 - m.x17) * (
m.x3 - m.x17) >= 6.25)
m.e10 = Constraint(expr= (m.x2 - m.x18) * (m.x2 - m.x18) + (m.x3 - m.x19) * (
m.x3 - m.x19) >= 3.24)
m.e11 = Constraint(expr= (m.x4 - m.x6) * (m.x4 - m.x6) + (m.x5 - m.x7) * (m.x5
- m.x7) >= 1.96)
m.e12 = Constraint(expr= (m.x4 - m.x8) * (m.x4 - m.x8) + (m.x5 - m.x9) * (m.x5
- m.x9) >= 5.29)
m.e13 = Constraint(expr= (m.x4 - m.x10) * (m.x4 - m.x10) + (m.x5 - m.x11) * (
m.x5 - m.x11) >= 1.21)
m.e14 = Constraint(expr= (m.x4 - m.x12) * (m.x4 - m.x12) + (m.x5 - m.x13) * (
m.x5 - m.x13) >= 3.61)
m.e15 = Constraint(expr= (m.x4 - m.x14) * (m.x4 - m.x14) + (m.x5 - m.x15) * (
m.x5 - m.x15) >= 6.76)
m.e16 = Constraint(expr= (m.x4 - m.x16) * (m.x4 - m.x16) + (m.x5 - m.x17) * (
m.x5 - m.x17) >= 3.61)
m.e17 = Constraint(expr= (m.x4 - m.x18) * (m.x4 - m.x18) + (m.x5 - m.x19) * (
m.x5 - m.x19) >= 1.44)
m.e18 = Constraint(expr= (m.x6 - m.x8) * (m.x6 - m.x8) + (m.x7 - m.x9) * (m.x7
- m.x9) >= 6.25)
m.e19 = Constraint(expr= (m.x6 - m.x10) * (m.x6 - m.x10) + (m.x7 - m.x11) * (
m.x7 - m.x11) >= 1.69)
m.e20 = Constraint(expr= (m.x6 - m.x12) * (m.x6 - m.x12) + (m.x7 - m.x13) * (
m.x7 - m.x13) >= 4.41)
m.e21 = Constraint(expr= (m.x6 - m.x14) * (m.x6 - m.x14) + (m.x7 - m.x15) * (
m.x7 - m.x15) >= 7.84)
m.e22 = Constraint(expr= (m.x6 - m.x16) * (m.x6 - m.x16) + (m.x7 - m.x17) * (
m.x7 - m.x17) >= 4.41)
m.e23 = Constraint(expr= (m.x6 - m.x18) * (m.x6 - m.x18) + (m.x7 - m.x19) * (
m.x7 - m.x19) >= 1.96)
m.e24 = Constraint(expr= (m.x8 - m.x10) * (m.x8 - m.x10) + (m.x9 - m.x11) * (
m.x9 - m.x11) >= 4.84)
m.e25 = Constraint(expr= (m.x8 - m.x12) * (m.x8 - m.x12) + (m.x9 - m.x13) * (
m.x9 - m.x13) >= 9)
m.e26 = Constraint(expr= (m.x8 - m.x14) * (m.x8 - m.x14) + (m.x9 - m.x15) * (
m.x9 - m.x15) >= 13.69)
m.e27 = Constraint(expr= (m.x8 - m.x16) * (m.x8 - m.x16) + (m.x9 - m.x17) * (
m.x9 - m.x17) >= 9)
m.e28 = Constraint(expr= (m.x8 - m.x18) * (m.x8 - m.x18) + (m.x9 - m.x19) * (
m.x9 - m.x19) >= 5.29)
m.e29 = Constraint(expr= (m.x10 - m.x12) * (m.x10 - m.x12) + (m.x11 - m.x13) *
(m.x11 - m.x13) >= 3.24)
m.e30 = Constraint(expr= (m.x10 - m.x14) * (m.x10 - m.x14) + (m.x11 - m.x15) *
(m.x11 - m.x15) >= 6.25)
m.e31 = Constraint(expr= (m.x10 - m.x16) * (m.x10 - m.x16) + (m.x11 - m.x17) *
(m.x11 - m.x17) >= 3.24)
m.e32 = Constraint(expr= (m.x10 - m.x18) * (m.x10 - m.x18) + (m.x11 - m.x19) *
(m.x11 - m.x19) >= 1.21)
m.e33 = Constraint(expr= (m.x12 - m.x14) * (m.x12 - m.x14) + (m.x13 - m.x15) *
(m.x13 - m.x15) >= 10.89)
m.e34 = Constraint(expr= (m.x12 - m.x16) * (m.x12 - m.x16) + (m.x13 - m.x17) *
(m.x13 - m.x17) >= 6.76)
m.e35 = Constraint(expr= (m.x12 - m.x18) * (m.x12 - m.x18) + (m.x13 - m.x19) *
(m.x13 - m.x19) >= 3.61)
m.e36 = Constraint(expr= (m.x14 - m.x16) * (m.x14 - m.x16) + (m.x15 - m.x17) *
(m.x15 - m.x17) >= 10.89)
m.e37 = Constraint(expr= (m.x14 - m.x18) * (m.x14 - m.x18) + (m.x15 - m.x19) *
(m.x15 - m.x19) >= 6.76)
m.e38 = Constraint(expr= (m.x16 - m.x18) * (m.x16 - m.x18) + (m.x17 - m.x19) *
(m.x17 - m.x19) >= 3.61)
m.e39 = Constraint(expr= m.x2 - m.x20 <= -1.2)
m.e40 = Constraint(expr= m.x3 - m.x21 <= -1.2)
m.e41 = Constraint(expr= m.x4 - m.x20 <= -0.6)
m.e42 = Constraint(expr= m.x5 - m.x21 <= -0.6)
m.e43 = Constraint(expr= m.x6 - m.x20 <= -0.8)
m.e44 = Constraint(expr= m.x7 - m.x21 <= -0.8)
m.e45 = Constraint(expr= m.x8 - m.x20 <= -1.7)
m.e46 = Constraint(expr= m.x9 - m.x21 <= -1.7)
m.e47 = Constraint(expr= m.x10 - m.x20 <= -0.5)
m.e48 = Constraint(expr= m.x11 - m.x21 <= -0.5)
m.e49 = Constraint(expr= m.x12 - m.x20 <= -1.3)
m.e50 = Constraint(expr= m.x13 - m.x21 <= -1.3)
m.e51 = Constraint(expr= m.x14 - m.x20 <= -2)
m.e52 = Constraint(expr= m.x15 - m.x21 <= -2)
m.e53 = Constraint(expr= m.x16 - m.x20 <= -1.3)
m.e54 = Constraint(expr= m.x17 - m.x21 <= -1.3)
m.e55 = Constraint(expr= m.x18 - m.x20 <= -0.6)
m.e56 = Constraint(expr= m.x19 - m.x21 <= -0.6)
m.e57 = Constraint(expr= m.x2 <= 9)
m.e58 = Constraint(expr= m.x3 <= 2)
m.e59 = Constraint(expr= m.x4 - m.x18 <= 0)
m.e60 = Constraint(expr= m.x12 - m.x16 <= 0)