Orthogonal Polynomials and Sturm Liouville operators Classical Orthogonal polynomials (e.g., Hermite, Legendre) are eigenfunctions of Sturm Liouville operators. For example, define  $L[u]=u''-xu'$, then the $n$-th order Hermite polynomial satisfies $LH_n (x) = -nH_n (x)$. 
The only proof of this miracle that I know of goes through solving the SL problem by power-series method, and on the other hand obtaining the orthogonal polynomials via Gram-Schmidt algorithm.
My question: I am going to give a talk about this property for graduate students. Is there another, perhaps a more "natural" explanation for this miracle?
(this is related to a previous question of mine, but not entirely the same)
 A: It seems to be most usual to take a spectral theory viewpoint (and it would be odd to avoid it because it is also most natural). One can give a proof of the spectral theorem for compact self-adjoint operators and motivate the introduction of compact operators by talking briefly about Green's functions and their relation to differential operators (one can mention the ubiquitous example of the Laplacian with Dirichlet/Neumann BC's and its Green's function). After proving the spectral theorem one can spend a little time developing the corresponding operators for Sturm-Liouville equations to justify the properties needed in the theorem do in fact hold. Then one has demonstrated existence of an orthonormal basis of eigenfunctions for the differential operators (since they are eigenfunctions for the compact operator). It may take more than one talk to do it properly.
Edit: To show that certain orthogonal polynomials do in fact form an orthonormal basis of eigenfunctions is a bit more involving for many cases (demonstrating completeness for example). But it should be clear by inspection that they are in fact eigenfunctions, which leaves showing they are orthogonal and complete. This might be explained in a few computations for some specific examples. Of course, the above suggestions could be substituted with more computations if you are inclined to talk more about this. These suggestions are heavily motivated by some lecture notes from one of T. Ransford's mini-courses 
 for grad students here. Starting on page 19 compact operators are briefly described, then the spectral theorem, then SL equations.
A: Let $Lu=u''-xu'$, and suppose $Lu=\lambda u$. Let $v$ be such that $v'=u$. By integrating the equation $Lu=\lambda u$, you find $Lv=(\lambda-1)v+C$, where $C$ is a constant. You can make this constant zero by choosing the integration constant in $v$ appropriately. Now observe that $L(1)=0$ and integrate successively.
