Okay, take a deep breath, and fasten your seat belt!
Let [tex]I= \int_{0}^{\infty} e ^ {-2.x^2} dx[/tex].
Of course we also can write [tex]I= \int_{0}^{\infty} e ^ {-2.y^2} dy[/tex] since I does not depend on the "dummy" variable of integration.
Multiplying those together, [tex]I^2= \left(\int_{0}^{\infty} e ^ {-2.x^2} dx\right)\left( \int_{0}^{\infty} e ^ {-2.y^2} dy\right)[/tex].
Since one integral involves on x and the other only y, and the neither limit of integration depends on the other variable, we can combine them into a double integral:
[tex]I^2= \int_{0}^{\infty} \int_{0}^{\infty} \left(e ^ {-2.x^2}\right)\left(e^{-2y^2}\right) dxdy=\int_{0}^{\infty} \int_{0}^{\infty} \left(e ^ {-2(x^2+ y^2)}\right)dxdy[/tex]
We can interpret that as integrating over the first quadrant (x and y both positive) and convert to polar coordinates.
[tex]e^{-2(x^2+ y^2)}[/tex] is easy: its [itex]e^{-2r^2}[/itex] in polar coordinates. The r-integral is from 0 to [itex]\infty[/itex] and the [itex]\theta[/itex] integral is from 0 to [itex]\frac{\pi}{2}[/itex] of course. Finally, the "differential of area" in polar coordinates is [itex]rdrd\theta[/itex].
We have
[tex]I^2= \int_0^{\frac{\pi}{2}}\int_0^\infty e^{-2r^2}rdrd\theta[/tex].
Of course, the [itex]\theta[/itex] integral just give [itex]\frac{\pi}{2}[/itex] so
[tex]I^2= \frac{\pi}{2}\int_0^\infty e^{-2r^2} rdr[/tex]
The crucial point is that we now have "rdr"! Make the substitution u= 2r2. Then du= 4 rdr. or rdr= (1/4)du. Now we have
[tex]I^2= \frac{\pi}{8}\int_0^\infty e^{-u} du[/tex].
I'll leave it to you to finish that.