Empezamos con la fórmula de Rodrigue \begin{eqnarray*} P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n} (x^2-1)^n. \end{eqnarray*} Utilizando la $n^{th}$ fórmula derivada tenemos la fórmula de Schlafli \begin{eqnarray*} P_n(x) = \frac{1}{2 \pi i} \int_C \frac{(t^2-1)^n}{2^n(t-x)^{n+1}} dt \end{eqnarray*} donde $C$ es cualquier contorno que rodea $t=x$ .
Ahora elige este contorno para que sea \begin{eqnarray*} t=x + (x^2-1)^{1/2} e^{i \phi} \end{eqnarray*} donde $\phi$ varía de $- \pi$ a $ \pi $ ... dando \begin{eqnarray*} P_n(x) = \frac{ 1}{ 2 \pi} \int_{- \pi }^{\pi} ( x+ (x^2-1)^{1/2} \cos( \phi ) )^{n} d \phi. \end{eqnarray*} Por último, obsérvese que el integrando es una función par de $\phi$ por lo que el intervalo se puede dividir por la mitad y se obtiene la fórmula deseada.
Edición: Tenga en cuenta que \begin{eqnarray*} t^2-1 &=& x^2-1 + 2x(x^2-1)^{1/2}e^{i \phi} + (x^2-1)e^{2 i \phi} \\ &=& (x^2-1)^{1/2}e^{i \phi}(2x + (x^2-1)^{1/2}(e^{ i \phi}+ e^{ -i \phi}))\\ &=& 2 (x^2-1)^{1/2}e^{i \phi}(x + (x^2-1)^{1/2}\cos \phi ).\\ \end{eqnarray*}