using JuMP, EAGO

                     m = Model()

                     EAGO.register_eago_operators!(m)

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

                     @variable(m, -5.118547195289687 <= q <= 7.5851414059354285)

                     add_NL_constraint(m, :(sigmoid(0.14385330157018572 + -0.4061277131705876*$(x[1]) + 0.38692435149352855*$(x[2]) + 0.23112211611402778*$(x[3])) + sigmoid(0.5778882819806599 + 0.8181171649048586*$(x[1]) + -0.37932684839331854*$(x[2]) + 0.7567763507538823*$(x[3])) + sigmoid(0.9224951510136457 + -0.8061937532687207*$(x[1]) + 0.8183596504529116*$(x[2]) + -0.17766965347054509*$(x[3])) + sigmoid(-0.41093962924162053 + -0.25704902180989464*$(x[1]) + 0.523731990676354*$(x[2]) + 0.7904456861039284*$(x[3])) - $q <= 0.0))

                     @objective(m, Min, q)

                     return m

                    