Factorizing a bivariate polynomial

Question

I have a bivariate polynomial for each $n=0,1,2...$ $$ f_n(x,y)=\sum _{k=0}^n \frac{(-1)^k}{2 k+1} \binom{n}{k} \left(x ^2-y ^2\right)^{2 n-2 k}\left([y ( x^2 -1) +x(1 -y^2 )]^{2 k+1}\\ \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\quad-[y ( x^2 -1) -x(1 -y^2 )]^{2 k+1}\right) $$ looking at low orders explicitly in Mathematica reveals that this polynomial contains $n+1$ factors of $(1-y^2)$. I know that one factor comes from the difference of two odd powers, but cannot see where the other $n$ factors come from. The factorized expressions appear quite complex, for example: $$ f_3=\frac{2}{35} \left(y ^2-1\right)^4 x \left(-35 y^6+35 \left(5 y ^2-1\right) x ^{10}+7 \left(15 y ^4-73 y ^2+3\right) x ^8+(y -1) (y +1) \left(5 y^4-388 y^2+5\right) x ^6-7 y ^2 \left(3 y ^4-73 y ^2+15\right) x ^4+35 y ^4 \left(y^2-5\right) x ^2+35 x ^{12}\right) $$

Assuming $f_n$ contains $(1-y^2)^{n+1}$, How can I find the formula for $f_n$ in its factorized form, for general $n$?

If it is any use, the series coefficients are related to the following integral: $$ \sum _{k=0}^n \frac{(-1)^k}{2 k+1}\binom{n}{k} z^{2k+1}=\int_0^z(1-u^2)^n\mathrm{d}u $$

Answer 1

So, \begin{split} f_n(\xi,\eta) &= (\xi^2-\eta^2)^{2n+1} \left( F\big(\frac{\eta ( \xi^2 -1) +\xi(1 -\eta^2)}{\xi^2-\eta^2}\big) - F\big(\frac{\eta ( \xi^2 -1) -\xi(1 -\eta^2)}{\xi^2-\eta^2}\big)\right) \\ &= (\xi^2-\eta^2)^{2n+1} \left( F\big(\frac{\eta\xi + 1}{\xi+\eta}\big) - F\big(\frac{\eta\xi - 1}{\xi-\eta}\big)\right), \end{split} where $$F(x) := \int_0^x (1-u^2)^n\,du.$$

To show that $(1-\eta^2)^{n+1}$ is a factor $f_n(\xi,\eta)$, it's enough to show that $(1-\eta^2)^n$ divides \begin{split} \frac{\partial}{\partial \eta} \left( F\big(\frac{\eta\xi + 1}{\xi+\eta}\big) - F\big(\frac{\eta\xi - 1}{\xi-\eta}\big) \right) &= \bigg(1-\big(\frac{\eta\xi + 1}{\xi+\eta}\big)^2\bigg)^n \frac{\xi^2-1}{(\xi+\eta)^2} - \bigg(1-\big(\frac{\eta\xi - 1}{\xi-\eta}\big)^2\bigg)^n\frac{\xi^2-1}{(\xi-\eta)^2} \\ &= (1-\eta^2)^n (\xi^2-1)^{n+1} \left( \frac{1}{(\xi+\eta)^{2n+2}} - \frac{1}{(\xi-\eta)^{2n+2}}\right), \end{split} and we see that it does.