Finding roots of multivariate function with given bounds?

[cvanaret]

Go for IntervalRootFinding.jl. It uses interval arithmetic and is able to separate and bound all the roots, while being robust to roundoff errors.