Abstract. The author proposes a method to calculate all the roots of a system of nonlinear equations inside a multidimensional interval. The main idea of the method is to subdivide the original interval into subintervals, which either do not have roots or satisfy the root uniqueness criterion based on the Krawczyk operator. An algorithm performing such a subdivision is presented.

Keywords: system of nonlinear equations, localization of roots, Krawczyk operator.