One dimension:

\[
G(x) = \frac{1}{\sqrt{2\pi \sigma^2}} e^{-\frac{x^2}{2\sigma^2}}
\]

Two dimensions:

\[
G(x, y) = \frac{1}{2\pi \sigma^2} e^{-\frac{x^2 + y^2}{2\sigma^2}}
\]

This is because:

\[
I(t)  = \int_{-t}^{t} e^{-x^2}dx \\
I^2(t)  = (\int_{-t}^{t} e^{-x^2}dx) \cdot (\int_{-t}^{t} e^{-y^2}dy) = \int_{-t}^{t} \int_{-t}^{t} e^{-(x^2+y^2)}dx dy \\
\int_{0}^{t} \int_{0}^{2\pi}e^{-r^2}r dr d\theta < I^2(t) < \int_{0}^{t\sqrt 2} \int_{0}^{2\pi}e^{-r^2}r dr d\theta \\
\int_{0}^{t} 2\pi e^{-r^2}  r dr d\theta < I^2(t) < \int_{0}^{t\sqrt 2} 2\pi e^{-r^2}r dr d\theta \\
\pi(1-e^{-t^2}) < I^2(t) < \pi(1-e^{-2t^2}) \\
\lim_{t \to +\infty} \pi(1-e^{-t^2}) < \lim_{t \to +\infty} I^2(t) < \lim_{t \to +\infty} \pi(1-e^{-2t^2}) \\
\pi < \lim_{t \to +\infty} I^2(t) < \pi
\]

so we obtain:

\[
\lim_{t \to +\infty} I^2(t) = \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e^{-(x^2+y^2)}dx dy = \pi \\
\lim_{t \to +\infty} I(t) = \int_{-\infty}^{\infty} e^{-x^2} dx = \sqrt \pi \\
\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{1}{2\pi \sigma^2} e^{-\frac{x^2+y^2}{2\sigma^2 }} dx dy = 1
\]