I'm not totally sure what "proof" you're talking about--I'm assuming it's the trick that the square of the integral can be written as a double integral in the familiar 2D $xy$-plane then converted to an easy integral using polar coordinates:
$$
\left(\int\limits_0^\infty e^{-4x^2}dx\right)^2 = \int_0^\infty e^{-4x^2}dx\int\limits_0^\infty e^{-4y^2}dy = \int\limits_0^\infty\int\limits_0^\infty e^{-4\left(x^2+y^2\right)}dxdy
$$
Before you convert to polar coordinates, you must understand that this region is only the first quadrant--which means the angle runs from $0$ to $\frac{\pi}{2}$ (as opposed to $[0,2\pi)$--as is the case when it's $-\infty$ to $\infty$).
Converting to polar coordinates gives:
$$
\int\limits_0^\infty\int\limits_0^\infty e^{-4\left(x^2+y^2\right)}dxdy = \int\limits_0^\frac{\pi}{2}d\theta\int\limits_0^\infty re^{-4r^2}dr
$$
The $r$-integral is easy through $u$-substitution and is $\frac{1}{8}$. The $\theta$-integral is independent of the $r$-integral so is just $\frac{\pi}{2}$. Giving:
$$
\left(\int\limits_0^\infty e^{-4x^2}dx\right)^2 = \frac{\pi}{16}
$$
And thus:
$$
\int\limits_0^\infty e^{-4x^2}dx = \frac{\sqrt{\pi}}{4}
$$
(just as can be retrieved from a u-substitution used from the result that $\int\limits_{-\infty}^\infty e^{-x^2}dx = \sqrt{\pi}$).