As mentioned this problem can be solved with Legendre Polynomials and their identities, but there is no need. There are many nice identities among the Legendre Polynomials and not a single one is used here.
The trick is to write the differencial equation as (in Sturm-Louiville form),
Thus,
![[(1-x^2)P_n'(x)]'+n(n+1)P_n(x)=0 [(1-x^2)P_n'(x)]'+n(n+1)P_n(x)=0](http://www.mathhelpforum.com/math-help/latex2/img/09749436fbc23be7c427fad817f4b1f3-1.gif)
Multiply by

to get,
![[(1-x^2)P_n'(x)]'P_m(x)+n(n+1)P_n(x)P_m(x)=0 [(1-x^2)P_n'(x)]'P_m(x)+n(n+1)P_n(x)P_m(x)=0](http://www.mathhelpforum.com/math-help/latex2/img/63811478048d2c0aec8d34ce219a0093-1.gif)
And integrate,
![\int_{-1}^1 [(1-x^2)P_n'(x)]'P_m(x) dx + n(n+1)\int_{-1}^1 P_n(x)P_m(x)dx=0 \int_{-1}^1 [(1-x^2)P_n'(x)]'P_m(x) dx + n(n+1)\int_{-1}^1 P_n(x)P_m(x)dx=0](http://www.mathhelpforum.com/math-help/latex2/img/9cd83d38656e2d701e77def453efbab3-1.gif)
Use integration by part on the first integral,

The first term vanishes,

Doing the similar steps with

instead we have,

Now subtract these two equations to get,
![[n(n+1)-m(m+1)]\int_{-1}^1 P_n(x)P_m(x) dx = 0 [n(n+1)-m(m+1)]\int_{-1}^1 P_n(x)P_m(x) dx = 0](http://www.mathhelpforum.com/math-help/latex2/img/b314287c0786955d1417f59c43b9efa2-1.gif)
Since

and they are integers we have,

.
Q.E.D.