= Solution
Multiply the two displayed equations by $e^{-x}$ and subtract the one with indices reversed. The product terms combine into the <Lagrange identity for a Sturm-Liouville operator>:
$$
\frac{d}{dx}\left[xe^{-x}(L_mL_n'-L_nL_m')\right]
+(n-m)e^{-x}L_mL_n=0.
$$
Integrating from zero to infinity gives
$$
\left[xe^{-x}(L_mL_n'-L_nL_m')\right]_0^\infty
+(n-m)\int_0^\infty e^{-x}L_mL_n\,dx=0.
$$
All factors $L_j,L_j'$ are <polynomials>. The endpoint term at zero vanishes because of the factor $x$, and the one at infinity vanishes because the <exponential function> dominates every <polynomial>. Therefore, if $m\ne n$,
$$
\boxed{\int_0^\infty e^{-x}L_m(x)L_n(x)\,dx=0}.
$$
The <integral> is absolutely convergent. This is <orthogonality> in the <weighted inner product> with weight $e^{-x}$; the <boundary> cancellation is the self-adjoint mechanism behind the <Laguerre differential operator on polynomials>.
Back to article page
