How to solve two symbolic nonlinear equations in Julia?

Yes, exactly.
IntervalBox(X, n) gives the n-fold cartesian product.