How to define range of the variables when numerically solving a system of algebraic equations?

In other words, what is the multi-variable version of the following

>> evalin(symengine, 'numeric::solve(sin(x) = x^2 - 1, x = 0..2)')

I tried this and that..

This works

syms x y; S = solve(tan(x) == exp(x), cos(y) == y/x)'

but when I specify the range for x, it does not

syms x y; S = solve(tan(x) == exp(x), cos(y) == y/x, x = 0..2)'

and this gives empty sum.. even before I set any ranges

evalin(symengine, 'numeric::solve(tan(x) = exp(x), cos(y) = y/x)')

Can someone type me an example of numeric multi-variable please?