Show that $$\begin{align*} \forall x \in [-1,1]: \int_0^{\pi} \frac{\sin^n \theta}{(1+x^2-2x \cdot \cos \theta)^{\frac{n}{2}}} \, d\theta &= c_n \tag{1} \\ \int_0^{\pi} \frac{\sin^{n+2} \theta}{(1+x^2-2x \cdot \cos \theta)^{\frac{n}{2}}} \, d\theta &= a_n \cdot x^2+b_n \tag{2} \end{align*}$$ where $a_n, b_n, c_n$ are constants (which do not depend on $x \in [-1,1]$), $n \geq 3$.
I tried several approaches (differentiation to show that the derivative of $(1)$ is equal to $0$, Weierstraß substitution, ...), but always got stuck. For example, one can show that the integrand in $(1)$ equals
$$\begin{align*} \frac{\sin^n \theta}{(1+x^2-2x \cdot \cos \theta)^{\frac{n}{2}}} &= \left( \sqrt{ \left( \frac{x-\cos \theta}{\sin \theta} \right)^2+1} \right)^{-n} \\ &= \left( \sqrt{4 \left( \frac{\sin \frac{\theta+\varrho}{2} \cdot \sin \frac{\theta-\varrho}{2}}{\sin \theta} \right)^2+1} \right)^{-n} \end{align*}$$ where $x=\cos \varrho$. I hoped to get some kind of symmetrization out of it, but (as far as I can see) it doesn't work.
Any ideas? (The aim is to find a rather quick or direct proof - a lengthy one is already known, using a recursive/inductive approach.)
Thanks!
The Gegenbauer polynomials were built to solve problems like this. Suppose $m\in\mathbb{N}$. We have $$\begin{eqnarray*} \int_0^{\pi} d\theta \, \frac{\sin^{n+2m} \theta}{(1-2x \cos \theta+x^2)^{n/2}} &=& \int_0^\pi d\theta \, \sin^{n+2m} \theta \sum_{k=0}^\infty C_k^{(n/2)}(\cos\theta) x^k \\ &=& \int_{-1}^1 d u \, (1-u^2)^{n/2-1/2} (1-u^2)^m \sum_{k=0}^\infty C_k^{(n/2)}(u) x^k \hspace{5ex} (u=\cos\theta) \\ &=& \int_{-1}^1 d u \, (1-u^2)^{n/2-1/2} \sum_{l=0}^m \beta_{2l} C_{2l}^{(n/2)}(u) \sum_{k=0}^\infty C_k^{(n/2)}(u) x^k \\ &=& \sum_{l=0}^m \sum_{k=0}^\infty \beta_{2l} x^k \underbrace{\int_{-1}^1 d u \, (1-u^2)^{n/2-1/2} C_{2l}^{(n/2)}(u) C_k^{(n/2)}(u)}_{\gamma_{2l}^{(n/2)} \delta_{2l,k}} \\ &=& \sum_{l=0}^m \beta_{2l} \gamma_{2l}^{(n/2)} x^{2l}. \end{eqnarray*}$$ This proves the claim for any $m\in\mathbb{N}$. The restriction $n\ge 3$ is too loose. The result holds for any $n$ such that $n>-1$.
For a given $m$ it is a straightforward exercise to find the $\beta$s and $\gamma$s. For example, using the fact that $1=C_0^{(n/2)}(u)$ and the normalization relation for the Gegenbauer polynomials we find $$\int_0^{\pi} d\theta \, \frac{\sin^{n} \theta}{(1-2x \cos \theta+x^2)^{n/2}} = \frac{\sqrt{\pi}\Gamma\left(\frac{n+1}{2}\right)}{\Gamma\left(\frac{n}{2}+1\right)}.$$
Some details: