using JuMP, EAGO

                     m = Model()

                     EAGO.register_eago_operators!(m)

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

                     @variable(m, -9.75097970893739 <= q <= 8.436779502208893)

                     add_NL_constraint(m, :(gelu(0.5456914178120993 + -0.994137029231358*gelu(0.4043679177252084 + -0.9159920740102909*$(x[1]) + 0.630569400327222*$(x[2]) + 0.15548796895771444*$(x[3])) + 0.38451384631825425*gelu(-0.7218630242098185 + -0.7615826451198555*$(x[1]) + 0.556650647741431*$(x[2]) + -0.7074053991868854*$(x[3])) + 0.9673084611668079*gelu(-0.4849043710003591 + 0.9895364493832828*$(x[1]) + -0.17212002157434236*$(x[2]) + -0.9444649970468406*$(x[3]))) + gelu(0.6143810270772536 + 0.32603207961941205*gelu(0.4043679177252084 + -0.9159920740102909*$(x[1]) + 0.630569400327222*$(x[2]) + 0.15548796895771444*$(x[3])) + 0.7660966505536861*gelu(-0.7218630242098185 + -0.7615826451198555*$(x[1]) + 0.556650647741431*$(x[2]) + -0.7074053991868854*$(x[3])) + 0.30918160231413605*gelu(-0.4849043710003591 + 0.9895364493832828*$(x[1]) + -0.17212002157434236*$(x[2]) + -0.9444649970468406*$(x[3]))) + gelu(-0.4079478457750749 + 0.06443443960646311*gelu(0.4043679177252084 + -0.9159920740102909*$(x[1]) + 0.630569400327222*$(x[2]) + 0.15548796895771444*$(x[3])) + 0.09351478982242645*gelu(-0.7218630242098185 + -0.7615826451198555*$(x[1]) + 0.556650647741431*$(x[2]) + -0.7074053991868854*$(x[3])) + -0.7258007123919756*gelu(-0.4849043710003591 + 0.9895364493832828*$(x[1]) + -0.17212002157434236*$(x[2]) + -0.9444649970468406*$(x[3]))) - $q <= 0.0))

                     @objective(m, Min, q)

                     return m

                    