Consider the following Gaussian Integral $$I = \int_{-\infty}^{\infty} e^{-x^2} \ dx$$

The usual trick to calculate this is to consider $$I^2 = \left(\int_{-\infty}^{\infty} e^{-x^2} \ dx \right) \left(\int_{-\infty}^{\infty} e^{-y^{2}} \ dy \right)$$

and convert to polar coordinates. We get $\sqrt{\pi}$ as the answer.

Is it possible to get the same answer by considering $I^{3}, I^{4}, \dots, I^{n}$?


Define $$ I_n = \prod_{i=1}^n \int^{\infty}_{-\infty} e^{-x_i^2}\,dx_i = \int_{\mathbb{R}^n} e^{-|x|^2}\,dx, $$ where $x = (x_1,\ldots,x_n)$ and $|x| = \sqrt{x_1^2 + \cdots + x_n^2}$.

By spherical coordinate integration transform in $\mathbb{R}^n$:

$$ I_n = \int_0^\infty\left\{\int_{\partial B(0,r)} e^{-|x|^2}\,dS\right\}\,dr. $$ For $e^{-|x|^2} = e^{-r^2}$ on $\partial B(0,r)$ which is a constant when the radius is fixed, hence above integral reads: $$ I_n = \int_0^\infty e^{-r^2}\left\{\int_{\partial B(0,r)} 1\,dS\right\}\,dr. $$ Now $$ \int_{\partial B(0,r)} 1\,dS = |\partial B(0,r)|= \omega_n r^{n-1}, $$ which is the surface area of the $(n-1)$-sphere with radius $r$, and $\omega_n$ is the surface area of the unit $(n-1)$-sphere: $$ \omega_n = \frac{n\pi^{n/2}}{\Gamma\left(1+\frac{n}{2}\right)}, $$ notice this can be computed by a recursive relation or taking derivative of the volume element of the $n$-ball.

Hence $$ I_n = \omega_n\int_0^\infty e^{-r^2}r^{n-1}\,dr = \frac{1}{2}\omega_n \Gamma\left(\frac{n}{2}\right) . $$

Now by Gamma function's property (proved by integration by parts): $$ \Gamma\left(1+\frac{n}{2}\right) = \frac{n}{2}\Gamma\left(\frac{n}{2}\right). $$ There goes your desired result: $$ I_n = \pi^{n/2}. $$


For three dimensions (if I figure out higher dimensions I will edit), one could try a similar argument by using spherical coordinates. If you can perform the integral, $$\int_0^\infty r^2 e^{-r^2} dr$$, then we can compute the desired integral. But This question addresses this. For this you need the spherical change of co-ordinates see this. You simply mimic the proof for the two dimensional case. I will leave the work to you.

Edit It seems in higher dimensions if one knows how to compute the integral, $$\int_0^\infty r^{n-1} e^{-r^2} dr$$ then we can compute the integral in essentially the same way. To see this you need higher dimensional spherical co-ordinates. For this see here.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.