Oh, sí, sustituyendo $t=-6x^2+x^4$ en el este y la "simplificación" se ve como un montón de diversión, de hecho. =D es broma. Recordemos la generación de la función de los polinomios de Legendre: $$\frac1{\sqrt{1-2xt+t^2}}=\sum^\infty_{n=0}P_n(x)\,t^n.$$ Multiplying with $t-x$ gives $$\frac{t-x}{\sqrt{1-2xt+t^2}}=-x+\sum^\infty_{n=1}[P_{n-1}(x)-x\,P_n(x)]\,t^n, \tag{1}$$ y la relación de recurrencia (http://mathworld.wolfram.com/LegendrePolynomial.html, (43)) da
$$P_{n-1}(x)-x\,P_n(x)=(n+1)\frac{x\,P_n(x)-P_{n-1}(x)}n,$$ so integrating this from $0$ to $t$, obtenemos
$$\sqrt{1-2xt+t^2}=1-xt+\sum^\infty_{n=2}\frac{x\,P_{n-1}(x)-P_n(x)}{n-1}t^n \tag{2}.$$ Substituting $x=3$ en (2), se llega a
$$\sqrt{1-6t+t^2}=1-3t+\sum^\infty_{n=2}\frac{3\,P_{n-1}(3)-P_n(3)}{n-1}t^n \tag{3}.$$ For the final result, we substitute $t=z^2$ para obtener
$$\sqrt{1-6z^2+z^4}=1-3z^2+\sum^\infty_{n=2}\frac{3\,P_{n-1}(3)-P_n(3)}{n-1}z^{2n}.$$ BTW, $P_n(3)$ es siempre un número entero, porque (http://mathworld.wolfram.com/LegendrePolynomial.html, (33))
$$P_n(3)=\sum^n_{k=0}\binom{n}{k}^2 2^k.$$ Los coeficientes en (3) parece ser enteros, demasiado, al menos eso es lo que los resultados numéricos sugieren (http://swift.sandbox.bluemix.net/#/repl/597d0432c3917c7502ddfa97).