  The authors present a unified method for calculating the zeros of the
classical orthogonal polynomials based upon the electrostatic interpretation
and its connection to the energy minimization problem. Examples are given with
error estimates for three cases of the Jacobi polynomials, three cases of the
Laguerre polynomials, and the Hermite polynomials. In the case of the Chebyshev
polynomials, exact errors are given.
