Quizá no sea la forma más directa. Dejemos que $x = \sin^2(\theta)$ . Entonces tenemos $dx = 2 \sin(\theta) \cos(\theta) d \theta$ . Por lo tanto, $$I_n = \int_0^{\pi/2} \sin^{2n}(\theta) \cos^{2n}(\theta) 2 \sin(\theta) \cos(\theta) d \theta$$ $$I_n = \frac1{2^{2n}} \int_0^{\pi/2} \sin^{2n+1}(2 \theta) d \theta = \frac1{2^{2n+1}} \int_0^{\pi} \sin^{2n+1}(\phi) d \phi = \frac1{2^{2n}} \int_0^{\pi/2} \sin^{2n+1}(\phi) d \phi$$ Denotemos $$J_{n} = \int_0^{\pi/2} \sin^{n}(\phi) d \phi.$$ Consulte aquí para la recurrencia que implica $J_n$ . Tenemos $J_{2n+1} = \frac{2n}{2n+1} J_{2n-1}$ . Por lo tanto, $$I_n = \frac{J_{2n+1}}{2^{2n}} = \frac14 \frac{2n}{2n+1} \frac{J_{2n-1}}{2^{2n-2}} = \frac14 \frac{2n}{2n+1} I_{n-1}.$$
EDITAR : La integralidad $I_n$ es $\beta(n+1,n+1)$ donde $\beta(x,y)$ es el $\beta$ función definida como $\displaystyle \beta(x,y) = \int_0^1 t^{x-1} (1-t)^{y-1} dt$ . $\beta(x,y)$ también tiene una representación en términos de $\Gamma$ función dada por $\displaystyle \beta(x,y) = \frac{\Gamma(x) \Gamma(y)}{\Gamma(x+y)}$ .
Por lo tanto, $$\begin{align} I_n & = \beta(n+1,n+1) = \frac{\Gamma(n+1) \Gamma(n+1)}{\Gamma(2n+2)} = \frac{n! n!}{(2n+1)!} = \frac{n^2}{(2n+1)(2n)}\frac{(n-1)! (n-1)!}{(2n-1)!}\\ & = \frac{n}{2(2n+1)} \frac{\Gamma(n) \Gamma(n)}{\Gamma(2n)} = \frac14 \frac{2n}{2n+1} \beta(n,n) = \frac14 \frac{2n}{2n+1} I_{n-1}. \end{align}$$
(Estaba añadiendo esta parte cuando m. k. publicó la respuesta).