using JuMP, EAGO

                     m = Model()

                     EAGO.register_eago_operators!(m)

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

                     @variable(m, -6.103907844236108 <= q <= 8.378675476959124)

                     add_NL_constraint(m, :((-0.8297607126794757 + -0.5064858067335005*(-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))*(1 + erf((-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))/sqrt(2)))/2 + -0.8619658571743893*(-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))*(1 + erf((-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))/sqrt(2)))/2 + 0.8415580445654567*(-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))*(1 + erf((-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))/sqrt(2)))/2)*(1 + erf((-0.8297607126794757 + -0.5064858067335005*(-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))*(1 + erf((-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))/sqrt(2)))/2 + -0.8619658571743893*(-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))*(1 + erf((-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))/sqrt(2)))/2 + 0.8415580445654567*(-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))*(1 + erf((-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))/sqrt(2)))/2)/sqrt(2)))/2 + (-0.0809505831099564 + 0.6185835236357913*(-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))*(1 + erf((-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))/sqrt(2)))/2 + -0.5288114683363174*(-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))*(1 + erf((-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))/sqrt(2)))/2 + -0.42070009030061817*(-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))*(1 + erf((-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))/sqrt(2)))/2)*(1 + erf((-0.0809505831099564 + 0.6185835236357913*(-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))*(1 + erf((-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))/sqrt(2)))/2 + -0.5288114683363174*(-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))*(1 + erf((-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))/sqrt(2)))/2 + -0.42070009030061817*(-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))*(1 + erf((-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))/sqrt(2)))/2)/sqrt(2)))/2 + (0.16011867292674786 + -0.14898607305663658*(-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))*(1 + erf((-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))/sqrt(2)))/2 + -0.8220978033005912*(-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))*(1 + erf((-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))/sqrt(2)))/2 + -0.8518875556726564*(-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))*(1 + erf((-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))/sqrt(2)))/2)*(1 + erf((0.16011867292674786 + -0.14898607305663658*(-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))*(1 + erf((-0.6188655598844357 + -0.38627203797486676*$(x[1]) + 0.4878933200630442*$(x[2]))/sqrt(2)))/2 + -0.8220978033005912*(-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))*(1 + erf((-0.7914291135758744 + 0.8682104207735972*$(x[1]) + 0.3684602934879275*$(x[2]))/sqrt(2)))/2 + -0.8518875556726564*(-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))*(1 + erf((-0.26405078449033637 + -0.6452849815457378*$(x[1]) + 0.9586530591287103*$(x[2]))/sqrt(2)))/2)/sqrt(2)))/2 - $q <= 0.0))

                     @objective(m, Min, q)

                     return m

                    