 Weighted orthogonality of Legendre polynomial derivatives

Differentiating the Legendre differential equation gives $((1-x^2{)}^2{R}_n^{\prime }{)}^{\prime }+(n-1)(n+2)(1-x^2)R_n=0$ for $R_n=P_n^{\prime }$. Multiply the equations for two distinct indices by the opposite polynomial and subtract. The endpoint term vanishes because of $(1-x^2{)}^2$, giving the displayed weighted orthogonality. The zero polynomial $R_0$ is not a nonzero eigenfunction.