Uno puede usar el truco con coordenadas esféricas:
$$
\Gamma\left(\frac{1}{2}\right) = \int_0^\infty \mathrm{e}^{-x} \frac{\mathrm{d} x}{\sqrt{x}} = \int_0^\infty \mathrm{e}^{-x} \mathrm{d} (2 \sqrt{x}) = \int_{-\infty}^\infty \mathrm{e}^{-u^2} \, \mathrm{d} u
$$
Entonces:
$$
\Gamma\left(\frac{1}{2}\right) \Gamma\left(\frac{1}{2}\right) = \int_{-\infty}^\infty \int_{-\infty}^\infty \mathrm{e}^{-x^2} \mathrm{e}^{-y^2} \, \mathrm{d}x \, \mathrm{d} y = \underbrace{\int_0^\infty r \mathrm{e}^{-r^2} \, \mathrm{d} r}_{\frac{1}{2}} \cdot \underbrace{\int_{0}^{2\pi} \, \mathrm{d} \phi}_{2 \pi} = \pi
$$
Alternativamente, se puede usar el segundo de Euler integral:
$$
\frac{\Gamma\left(\frac{1}{2}\right) \Gamma\left(\frac{1}{2}\right) }{\Gamma(1)} = \int_0^1 t^{-1/2} (1-t)^{-1/2} \mathrm{d}t = \int_0^1 \mathrm{d} \left(2\arcsin\left(\sqrt{t}\right)\right) = \pi
$$
Ahora, el uso de $\Gamma(1) = 1$ el resultado de la siguiente manera.
Otro método es utilizar la duplicación de identidad:
$$
\Gamma(2s) = \frac{2^{2s-1}}{\sqrt{\pi}} \Gamma(s) \Gamma\left(s+\frac{1}{2}\right)
$$
en $ s= \frac{1}{2}$.