Aloha :)
$$I^2=I\cdot I=\int\limits_{-\infty}^\infty e^{-x^2}dx\cdot\int\limits_{-\infty}^\infty e^{-y^2}dy=\int\limits_{-\infty}^\infty\int\limits_{-\infty}^\infty e^{-x^2}e^{-y^2}\,dx\,dy=\int\limits_{-\infty}^\infty\int\limits_{-\infty}^\infty e^{-(x^2+y^2)}\,dx\,dy$$
Wir gehen nun über zu Polarkoordinaten:
$$\binom{x}{y}=\binom{r\cos\varphi}{r\sin\varphi}\quad;\quad dx\,dy=r\,dr\,d\varphi\quad;\quad r\in[0;\infty)\quad;\quad \varphi\in[0;2\pi]$$Die haben nämlich den Charme, dass \(x^2+y^2=r^2\) gilt und wir durch die Verzerrung des Flächenelements \(dx\,dy=r\,dr\,d\varphi\) den Faktor \(r\) vor die \(e\)-Funktion schreiben können, wodurch sich das Integral sofort hinschreiben lässt:
$$I^2=\int\limits_0^{2\pi}d\varphi\int\limits_0^\infty r\,e^{-r^2}\,dr=2\pi\cdot\left[-\frac{e^{-r^2}}{2}\right]_0^\infty=2\pi\left(-0+\frac{e^{0}}{2}\right)=2\pi\cdot\frac{1}{2}=\pi$$Damit ist:$$\int\limits_{-\infty}^\infty e^{-x^2}dx=\sqrt\pi$$