Electrostatic Interpretation of Zeros of Orthogonal Polynomials

We study the differential equation $ - (p(x) y')' + q(x) y' = \lambda y,$ where $p(x)$ is a polynomial of degree at most 2 and $q(x)$ is a polynomial of degree at most 1. This includes the classical Jacobi polynomials, Hermite polynomials, Legendre polynomials, Chebychev polynomials and Laguerre polynomials. We provide a general electrostatic interpretation of zeros of such polynomials: the set of real numbers $\left\{x_1, \dots, x_n\right\}$ satisfies $$ p(x_i) \sum_{k = 1 \atop k \neq i}^{n}{\frac{2}{x_k - x_i}} = q(x_i) - p'(x_i) \qquad \mbox{for all}~ 1\leq i \leq n$$ if and only if they are zeros of a polynomial solving the differential equation. We also derive a system of ODEs depending on $p(x),q(x)$ whose solutions converge to the zeros of the orthogonal polynomial at an exponential rate.