How to prove
This is an old favorite of mine.
Now change to polar coordinates
The θ integral just gives 2π, while the r integral succumbs to the substitution u=r2
So I=√π and your integral is half this by symmetry
I have always wondered if somebody found it this way, or did it first using complex variables and noticed this would work.