I was trying to find out what the function
$$ f: x, y, x_0, y_0, \sigma_x, \sigma_y \mapsto \frac{1}{2\pi \sigma_x \sigma_y} \exp (\frac{-1}{2} ((\frac{x-x_0}{\sigma_x})^2 + (\frac{y-y_0}{\sigma_y})^2)) $$
is. (No domains were indicated for any of the six variables.)
So I set $x_0, y_0 = 0$, $\sigma_x, \sigma_y =1$ and viewed $f$ as a function in two variables.
$$ f: x, y \mapsto \frac{1}{2\pi} \exp(\frac{-1}{2}(x^2+y^2)) $$
In this constellation, $f$ takes its maximum at $(x,y) = (0,0)$ and decays rapidly as $(x,y) \to \infty$. It is fast-forward to prove that for general $(x_0, y_0) \in \mathbb{R}^2$ and $\sigma_x, \sigma_y > 0$, the implicitly defined function $f: \mathbb{R}^2 \to \mathbb{R}$ has a unique maximum at $(x_0, y_0)$ and decays rapidly as $(x, y) \to \infty$.
I suspect that such $f$ (with given $x_0, y_0, \sigma_x, \sigma_y$) is a density function. In order to prove this, I have to show that
$$ \int_{(x,y) \in \mathbb{R}^2} f(x,y) \, d(x,y) = 1. $$
This task can be reduced to the task of showing
$$ \int_{-\infty}^{+\infty} \exp( \frac{-1}{2} (\frac{x-x_0}{\sigma_x})^2) \, dx = \sqrt{2\pi} $$
I was confident to show this if $x_0 = 0$ and $\sigma_x = 1$. In this case the integral becomes the
$$ \int_{-\infty}^{+\infty} \exp(\frac{-1}{2} x^2) \, dx $$
from the title of this question.
Does someone have a hint for me? How can I prove (in the last, simplified case) that the integral is $\sqrt{2\pi}$? I thought that I would see an easy primitive, but I am missing a factor $x$ on the exponential for this trick.