A long while ago (at university) I learned the following (fairly standard) method for evaluating $I:=\int_{-\infty}^{\infty}e^{-z^2}dz$ by squaring $I$ then using polar coordinates:
$$ \begin{align} I^2 &= \int_{-\infty}^{\infty}e^{-z^2}dz \int_{-\infty}^{\infty}e^{-z^2}dz \\ &= 4\int_{0}^{\infty}e^{-x^2}dx \int_{0}^{\infty}e^{-y^2}dy \\ &= 4\int_{0}^{\infty}\int_{0}^{\infty}e^{-x^2}e^{-y^2}dxdy \\ &= 4\int_{0}^{\infty}\int_{0}^{\infty}e^{-(x^2+y^2)}dxdy \\ &= 4\int_{0}^{\frac{\pi}{2}}\int_{0}^{\infty}e^{-r^2}rdrd\theta \\ &= 2\pi\int_{0}^{\infty}re^{-r^2}dr \\ &= \pi \left[ -e^{-r^2} \right]_{0}^{\infty} \\ &= \pi \end{align} \\ \\ \implies I=\sqrt{\pi} $$
I remember being amazed at the ingenuity of that first line: squaring $I$ in order to turn the single integral into a double integral.
I am curious whether any other methods are known for evaluating $I$. Especially methods that do not begin by squaring $I$.
I have not been able to do anything worthwhile with it myself. My research has yielded this method that does not require polar coordinates but still begins by squaring $I$. I have also considered complex integration using a clever contour but have not managed to complete a method.
Does anybody know of any other methods?