Enclosing all zeros of an analytic function — A rigorous approach

We present a method to find all zeros of an analytic function in a rectangular domain. The approach is based on finding guaranteed enclosures rather than approximations of the zeros. Well-isolated simple zeros are determined fast and with high accuracy. Clusters of zeros can in many cases be distinguished from multiple zeros by applying the argument principle to sufficiently high-order derivatives of the function. We illustrate the proposed method through five examples of varying levels of complexity.