A change of variables reduces it to evaluating the Gaussian integral (integral of e^{-x^2} over the real line). This is well known to be equal to sqrt(pi), and is most commonly computed using polar coordinates as in https://en.wikipedia.org/wiki/Gaussian_integral#By_polar_coordinates ).