Integración por partes de los rendimientos
$$
\frac{1}{k!}\int_x^\infty e^{-t}\,t^k\,\mathrm{d}t=\frac{1}{k!}x^ke^{-x}+\frac{1}{(k-1)!}\int_x^\infty e^{-t}\,t^{k-1}\,\mathrm{d}t\etiqueta{1}
$$
La iteración $(1)$ da
$$
\frac{1}{n!}\int_x^\infty e^{-t}\,t^n\,\mathrm{d}t=e^{-x}\sum_{k=0}^n\frac{x^k}{k!}\la etiqueta{2}
$$
Por lo tanto, obtenemos
$$
e^{-n}\sum_{k=0}^n\frac{n^k}{k!}=\frac{1}{n!}\int_n^\infty e^{-t}\,t^n\,\mathrm{d}t\etiqueta{3}
$$
Ahora, voy a reproducir parte del argumento que doy aquí, que desarrolla un completo asintótica de expansión. Además, se incluyen algunas estimaciones de error que anteriormente estaban desaparecidos.
$$
\begin{align}
\int_n^\infty e^{-t}\,t^n\,\mathrm{d}t
&=n^{n+1}e^{-n}\int_0^\infty e^{ns}\,(s+1)^n\,\mathrm{d}\\
&=n^{n+1}e^{-n}\int_0^\infty e^{-n(s-\log(1+s)}\,\mathrm{d}\\
&=n^{n+1}e^{-n}\int_0^\infty e^{-nu^2/2}\,\,\mathrm{d}u\etiqueta{4}
\end{align}
$$
donde $t=n(s+1)$ and $u^2/2=s-\log(1+s)$.
Tenga en cuenta que $\frac{ss'}{1+s}=u$; thus, when $s\ge1$, $s'\le2u$. Esto lleva a que el obligado
$$
\begin{align}
\int_{s\ge1} e^{-nu^2/2}\,\,\mathrm{d}u
&\le\int_{3/4}^\infty e^{-nu^2/2}\,2u\,\mathrm{d}u\\
&=\frac2ne^{-\frac98n}\etiqueta{5}
\end{align}
$$
$(5)$ también muestran que
$$
\int_{s\ge1}e^{-nu^2/2}\,\mathrm{d}u\le\frac2ne^{-\frac98n}\etiqueta{6}
$$
Para $|s|<1$, obtenemos
$$
u^2/2=s-\log(1+s)=s^2/2-s^3/3+s^4/4-\dots\etiqueta{7}
$$
Podemos invertir la serie para obtener $s'=1+\frac23u+O(u^2)$. Por lo tanto,
$$
\begin{align}
\int_0^\infty e^{-nu^2/2}\,\,\mathrm{d}u
&=\int_{s\in[0,1]} e^{-nu^2/2}\,\,\mathrm{d}u+\color{red}{\int_{s>1} e^{-nu^2/2}\,\,\mathrm{d}u}\\
&=\int_0^\infty\left(1+\frac23u\right)e^{-nu^2/2}\,\mathrm{d}u-\color{darkorange}{\int_{s>1}\left(1+\frac23u\right)e^{-nu^2/2}\,\mathrm{d}u}\\
&+\int_0^\infty e^{-nu^2/2}\,O(u^2)\,\mathrm{d}u-\color{darkorange}{\int_{s>1} e^{-nu^2/2}\,O(u^2)\,\mathrm{d}u}\\
&+\color{red}{\int_{s>1} e^{-nu^2/2}\,\,\mathrm{d}u}\\
&=\sqrt{\frac{\pi}{2n}}+\frac2{3n}+O\a la izquierda(n^{-3/2}\right)\etiqueta{8}
\end{align}
$$
El rojo y el naranja de las integrales disminuir exponencialmente por $(5)$ and $(6)$.
Conectar $(8)$ into $(4)$ rendimientos
$$
\int_n^\infty e^{-t}\,t^n\,\mathrm{d}t=\left(\sqrt{\frac{\pi n}{2}}+\frac23\right)\,n^ne^{-n}+O(n^{n-1/2}e^{-n})\etiqueta{9}
$$
El argumento anterior puede ser utilizado para comprobar la aproximación de Stirling, que dice que
$$
n!=\sqrt{2\pi n}\,n^ne^{-n}+O(n^{n-1/2}e^{-n})\etiqueta{10}
$$
La combinación de $(9)$ and $(10)$ rendimientos
$$
\begin{align}
e^{-n}\sum_{k=0}^n\frac{n^k}{k!}
&=\frac{1}{n!}\int_n^\infty e^{-t}\,t^n\,\mathrm{d}t\\
&=\frac12+\frac{2/3}{\sqrt{2\pi n}}+O(n^{-1})\etiqueta{11}
\end{align}
$$