This paper presents a method of numerically computing zeros of an analytic function for the specific application of computing eigenvalues of the Sturm-Liouville problem. The Sturm-Liouville problem is an infinite dimensional eigenvalue problem that often arises in solving partial differential equations, including the heat and wave equations. To compute eigenvalues of the Sturm-Liouville problem...