I would like to show that given a function $f:\mathbb{R}\to \mathbb{R}$ such that is continuos, $|f(x)| < M$ and $\exists z \in \mathbb{R} $ such that $$-\infty < \lim_{\sigma \to +\infty} \int_{-\infty}^{+\infty} f(z + \varepsilon) \frac{1}{\sqrt{2\pi\sigma^2}} e^{-\frac{\varepsilon^2}{2\sigma^2}}\mathrm{d}\varepsilon < +\infty$$ then:
$$ \forall x,y \in \mathbb{R} \lim_{\sigma \to +\infty} \int_{-\infty}^{+\infty} f(x + \varepsilon) \frac{1}{\sqrt{2\pi\sigma^2}} e^{-\frac{\varepsilon^2}{2\sigma^2}}\mathrm{d}\varepsilon = \lim_{\sigma \to +\infty} \int_{-\infty}^{+\infty} f(y + \varepsilon) \frac{1}{\sqrt{2\pi\sigma^2}} e^{-\frac{\varepsilon^2}{2\sigma^2}}\mathrm{d}\varepsilon $$
I started from the following:
if $$\lim_{\sigma \to +\infty} \int_{-\infty}^{+\infty} f(x+\varepsilon) \frac{1}{2\pi\sigma^2} e^{-\frac{\varepsilon^2}{2\sigma^2}} \mathbb{d}\varepsilon = c$$
then we want to show that
$$\lim_{\sigma \to +\infty} \int_{-\infty}^{+\infty} (f(y+\varepsilon)-c) \frac{1}{2\pi\sigma^2} e^{-\frac{\varepsilon^2}{2\sigma^2}} \mathbb{d}\varepsilon$$
Let $\delta$ be arbitrary then there is $\theta = \theta_{\delta} > 0$ such that
$$\int_{-\theta}^{+\theta} (f(y+\varepsilon)-c) \frac{1}{2\pi\sigma^2} e^{-\frac{\varepsilon^2}{2\sigma^2}} \mathbb{d}\varepsilon < \int_{-\theta}^{+\theta} (M-c) \frac{1}{2\pi\sigma^2} e^{-\frac{\varepsilon^2}{2\sigma^2}} \mathbb{d}\varepsilon$$
(\varepsilon is bounded by -\theta and +\theta) and
$$\lim_{\sigma \to +\infty} \int_{-\theta}^{+\theta} (M-c) \frac{1}{2\pi\sigma^2} e^{-\frac{\varepsilon^2}{2\sigma^2}} \mathbb{d}\varepsilon = 0$$
Now It is more difficult to show that the integral between $+\theta$ and $+\infty$ (or the integral between $-\infty$ and $-\theta$ goes to 0. I think in someway the fact that with $x$ it was bounded.