Double integral for gamma function I have the following double integral $$\int_0^{2\pi}\int_0^\infty r^{p+q-1}\cos^{p-1}({\theta})\sin^{q-1}({\theta})e^{-r(\sin{\theta}+{\cos{\theta}})}drd\theta $$ However, I can't seem to be able to isolate the required terms to prove the identity $$\Gamma(p)\Gamma(q)=B(p,q)\Gamma(p+q)$$ where $$ \Gamma(p) \equiv \int_0^\infty x^{p-1}e^{-x}dx $$ and $$ B(p,q) \equiv \int_0^1 x^{p-1}(1-x)^{q-1}dx = 2\int_0^\frac{\pi}{2}\sin^{2p-1}({2\theta})\cos^{2q-1}({2\theta})d\theta$$
Many thanks.
 A: Write the gamma function $\Gamma(p)$ as
$$
\Gamma(p)=\int_0^{\infty}t^{p-1}\operatorname{e}^{-t}\operatorname{d}\!t=\int_0^{\infty}x^{2p-2}\operatorname{e}^{-x^2}2x\operatorname{d}\!x=2\int_0^{\infty}x^{2p-1}\operatorname{e}^{-x^2}\operatorname{d}\!x
$$
changing the integration variable from $t$ to $x^2$.
Write the product of $\Gamma(p)\Gamma(q)$ as
$$
\Gamma(p)\Gamma(q)=4\int_0^{\infty}\int_0^{\infty}x^{2p-1}y^{2q-1}\operatorname{e}^{-x^2-y^2}\operatorname{d}\!x\operatorname{d}\!y.\tag 1
$$
We now rewrite the (1) in polar coordinates $(r,\theta)$ with $x=r\cos\theta,\,y=r\sin\theta$, $x^2+y^2=r^2$ and Jacobian $J=r$; in the polar coordnates the integration is over the entire first quadrant, i.e. $0\le r <\infty,\,0\le\theta\le \pi/2$. We get
$$\begin{align}
\Gamma(p)\Gamma(q)&=4\int_0^{\pi/2}\int_0^{\infty}(r\cos\theta)^{2p-1}(r\sin\theta)^{2q-1}\operatorname{e}^{-r^2}r\operatorname{d}\!r\operatorname{d}\!\theta\\
&=4\int_0^{\pi/2}\int_0^{\infty}(\cos\theta)^{2p-1}(\sin\theta)^{2q-1}\operatorname{e}^{-r^2}r^{2(p+q)-1}\operatorname{d}\!r\operatorname{d}\!\theta\\
&=4\int_0^{\infty}\operatorname{e}^{-r^2}r^{2(p+q)-1}\operatorname{d}\!r\int_0^{\pi/2}(\cos\theta)^{2p-1}(\sin\theta)^{2q-1}\operatorname{d}\!\theta.
\end{align}
$$
Making the substitution $r^2=t$ in the $r$ integration, we find
$$
\int_0^{\infty}\operatorname{e}^{-r^2}r^{2(p+q)-1}\operatorname{d}\!r=\frac{1}{2}\int_0^{\infty}\operatorname{e}^{-t}t^{(p+q)-1}\operatorname{d}\!t=\frac{1}{2}\Gamma(p+q)
$$
so we have
$$
\Gamma(p)\Gamma(q)=\frac{1}{2}\Gamma(p+q)4\int_0^{\pi/2}(\cos\theta)^{2p-1}(\sin\theta)^{2q-1}\operatorname{d}\!\theta
$$
that is
$$
\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}=B(p,q)=2\int_0^{\pi/2}(\cos\theta)^{2p-1}(\sin\theta)^{2q-1}\operatorname{d}\!\theta.
$$
