Básicamente, usted está considerando la integral
$$ I := \frac{1}{3}\int_{-\infty}^{\infty}x^n f(x) dx= \frac{1}{3}\int_{-\infty}^{\infty}x^n \int_{x}^{x+3}e^{-t^2}dt \,dx.$$
Cambiando el orden de integración de los rendimientos
$$ I = \frac{1}{3}\int_{-\infty}^{\infty}e^{-t^2} \int_{t-3}^{t}x^ndx \,dt $$
$$ = \frac{1}{3}\int_{-\infty}^{\infty}e^{-t^2}\left(\frac{t^{n+1}-(t-3)^{n+1}}{n+1} \right) \,dt $$
$$ = \frac{1}{3(n+1)}\int_{-\infty}^{\infty}t^{n+1} e^{-t^2} dt - \frac{1}{3(n+1)}\int_{-\infty}^{\infty}(t-3)^{n+1} e^{-t^2} dt dt $$
$$\implies I = \frac{1}{3(n+1)} I_1 + \frac{1}{3(n+1)} I_2. $$
Ahora, tenga en cuenta esto, para $n\in \mathbb{N} \cup \left\{ 0\right\}$, $I_1$ es convergente (ver el análisis para el caso de $I_2$ al final de la respuesta) y pueden ser evaluados en términos de la función gamma como
$$ I_1 = \frac{1}{2}\, \left( 1+ \left( -1 \right)^{n+1} \right) \Gamma \left( \frac{n}{2}
+1 \right).$$
El examen de la convergencia de $I_2$ $n\in \mathbb{N} \cup \left\{ 0\right\}$, tenemos
$$ I_2 = \int_{-\infty}^{\infty}(t-3)^{n+1} e^{-t^2} dt = \int_{-\infty}^{\infty}y^{n+1} e^{-(y+3)^2} dy $$
$$ = (-1)^n \int_{0}^{\infty}y^{n+1} e^{-(y-3)^2} dy + \int_{0}^{\infty}y^{n+1} e^{-(y+3)^2} dy $$
$$ I_2 = I_{21} + I_{22}. $$
Ahora, tanto de las integrales de $I_{21}$ $I_{22}$ convergen, desde
$$ y^{n+1}e^{-(y-3)^2} \leq e^{y} e^{-(y-3)^2}, \quad y^{n+1}e^{-(y+3)^2}\leq e^{y} e^{-(y+3)^2} $$
y
$$ \int_{0}^{\infty} e^{y} e^{-(y-3)^2}dy <\infty,\quad \int_{0}^{\infty} e^{y} e^{-(y+3)^2}dy <\infty. $$
Usted puede comprobar Gaussiano Integrales para las dos últimas desigualdades.
Nota: Para el caso de $n$ es entero negativo, tenemos
1) Si $n=-1$, $I_1$ $I_2$ convergen.
2) Si $n$ es un entero negativo, incluso, a continuación, $I_1$ $I_2$ son indefinidos.
3) Si $n$ es negativo entero impar, entonces $I_1$ $I_2$ diverge a infinito.