  This paper develops a practical numerical method for locating the zeros and
poles of a meromorphic function, based upon the generalized argument principle.
Sensitivity is analysed carefully to assure the algorithm stability. A novel
approach to the error estimation is also proposed. Numerical examples are given
to illustrate the effectiveness and practicability, with emphasis on the
potential applications in plasma physics.
