using JuMP, EAGO

                     m = Model()

                     EAGO.register_eago_operators!(m)

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

                     @variable(m, -11.981958290664823 <= q <= 14.886364381032921)

                     add_NL_constraint(m, :((0.5346214578242647 + 0.7328398696189291*$(x[1]) + 0.5275606331718761*$(x[2]) + -0.05282758308498137*$(x[3]) + -0.6801317780174543*$(x[4]) + 0.6177845475163379*$(x[5]))*(1 + erf((0.5346214578242647 + 0.7328398696189291*$(x[1]) + 0.5275606331718761*$(x[2]) + -0.05282758308498137*$(x[3]) + -0.6801317780174543*$(x[4]) + 0.6177845475163379*$(x[5]))/sqrt(2)))/2 + (-0.36399156729020854 + 0.9651821996491203*$(x[1]) + -0.1221090890096228*$(x[2]) + -0.0037390408554469623*$(x[3]) + -0.25471778237955167*$(x[4]) + 0.004904328924779211*$(x[5]))*(1 + erf((-0.36399156729020854 + 0.9651821996491203*$(x[1]) + -0.1221090890096228*$(x[2]) + -0.0037390408554469623*$(x[3]) + -0.25471778237955167*$(x[4]) + 0.004904328924779211*$(x[5]))/sqrt(2)))/2 + (0.6792969956775154 + -0.9847763250624375*$(x[1]) + -0.7653951229931786*$(x[2]) + 0.7304491787436271*$(x[3]) + -0.81259603694894*$(x[4]) + -0.9953733056123499*$(x[5]))*(1 + erf((0.6792969956775154 + -0.9847763250624375*$(x[1]) + -0.7653951229931786*$(x[2]) + 0.7304491787436271*$(x[3]) + -0.81259603694894*$(x[4]) + -0.9953733056123499*$(x[5]))/sqrt(2)))/2 + (-0.022503806788797487 + -0.877492362258455*$(x[1]) + 0.9093714868360583*$(x[2]) + -0.36771946270068234*$(x[3]) + 0.3633871862886231*$(x[4]) + -0.6803298907372901*$(x[5]))*(1 + erf((-0.022503806788797487 + -0.877492362258455*$(x[1]) + 0.9093714868360583*$(x[2]) + -0.36771946270068234*$(x[3]) + 0.3633871862886231*$(x[4]) + -0.6803298907372901*$(x[5]))/sqrt(2)))/2 + (0.6247799657612747 + 0.2803258134848612*$(x[1]) + 0.15247610330367944*$(x[2]) + -0.6290763102374144*$(x[3]) + -0.40592776579339773*$(x[4]) + -0.5176681326197756*$(x[5]))*(1 + erf((0.6247799657612747 + 0.2803258134848612*$(x[1]) + 0.15247610330367944*$(x[2]) + -0.6290763102374144*$(x[3]) + -0.40592776579339773*$(x[4]) + -0.5176681326197756*$(x[5]))/sqrt(2)))/2 - $q <= 0.0))

                     @objective(m, Min, q)

                     return m

                    