Chebyshev subdivision and reduction methods for solving multivariable systems of equations

We present a new algorithm for finding isolated zeros of a system of real-valued functions in a bounded interval in $R^n$. It uses the Chebyshev proxy method combined with a mixture of subdivision, reduction methods, and elimination checks that leverage special properties of Chebyshev polynomials. We prove the method has quadratic convergence locally near simple zeros of the system. It also finds all nonsimple zeros, but convergence to those zeros is not guaranteed to be quadratic. We also analyze the arithmetic complexity and the numerical stability of the algorithm and provide numerical evidence in dimensions up to five that the method is both fast and accurate on a wide range of problems. Our tests show that the algorithm outperforms other standard methods on the problem of finding all real zeros in a bounded domain. Our Python implementation of the algorithm is publicly available at https://github.com/tylerjarvis/RootFinding.