using JuMP, EAGO

                     m = Model()

                     EAGO.register_eago_operators!(m)

                     @variable(m, -1 <= x[i=1:2] <= 1)

                     @variable(m, -2.6308206135564163 <= q <= 3.997051239074799)

                     add_NL_constraint(m, :((0.8157606255059098 + 0.2914814987498171*$(x[1]) + -0.962674206030298*$(x[2]))/(1 + exp(-(0.8157606255059098 + 0.2914814987498171*$(x[1]) + -0.962674206030298*$(x[2])))) + (0.7907664778929373 + 0.4020245577506052*$(x[1]) + -0.1936195018990734*$(x[2]))/(1 + exp(-(0.7907664778929373 + 0.4020245577506052*$(x[1]) + -0.1936195018990734*$(x[2])))) + (-0.5840361484149175 + 0.5201746807648506*$(x[1]) + 0.3251106634438119*$(x[2]))/(1 + exp(-(-0.5840361484149175 + 0.5201746807648506*$(x[1]) + 0.3251106634438119*$(x[2])))) + (-0.33937564222473826 + -0.3486039500062561*$(x[1]) + 0.27024686767089534*$(x[2]))/(1 + exp(-(-0.33937564222473826 + -0.3486039500062561*$(x[1]) + 0.27024686767089534*$(x[2])))) - $q <= 0.0))

                     @objective(m, Min, q)

                     return m

                    