using JuMP, EAGO

                     m = Model()

                     EAGO.register_eago_operators!(m)

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

                     @variable(m, -6.482349417201714 <= q <= 6.592803073996197)

                     add_NL_constraint(m, :(sigmoid(-0.7946398993851154 + 0.38713964342476537*sigmoid(0.024809189775459206 + 0.8902991643012754*$(x[1]) + 0.015536230415071905*$(x[2]) + -0.8172113269588737*$(x[3]) + 0.9746036973952403*$(x[4])) + -0.9438792568465639*sigmoid(-0.017899887293421735 + -0.33137455646820024*$(x[1]) + 0.5594894653083848*$(x[2]) + -0.34768764708208755*$(x[3]) + 0.43788143230656695*$(x[4]))) + sigmoid(0.8346867556655706 + -0.9900146410670576*sigmoid(0.024809189775459206 + 0.8902991643012754*$(x[1]) + 0.015536230415071905*$(x[2]) + -0.8172113269588737*$(x[3]) + 0.9746036973952403*$(x[4])) + -0.7397521454164151*sigmoid(-0.017899887293421735 + -0.33137455646820024*$(x[1]) + 0.5594894653083848*$(x[2]) + -0.34768764708208755*$(x[3]) + 0.43788143230656695*$(x[4]))) - $q <= 0.0))

                     @objective(m, Min, q)

                     return m

                    