So, I've been studying the power series $$f(x)=\sum_{n=0}^\infty n!x^n$$ and through algebraic manipulation I've found two ways to get a closed form for it.
The first notices that $$1+xf(x)+x^2f'(x)=f(x)$$ and then proceeds to solve this linear ODE assuming that $\displaystyle\lim_{x\to\pm 0}f(x)=1$.
The answer here is $f(x)=\frac{1}{x}e^{-\frac{1}{x}}\text{Ei}(x)$, where Ei is defined here.
The second method notices that $$n!=\int_0^\infty e^{-t}t^n dt$$ so $$f(x)=\int_0^\infty e^{-t}\sum_{n=0}^\infty (tx)^n dt=\int_0^\infty \frac{e^{-t}}{1-xt}dt$$ which again eventually resolves to $$f(x)=\frac{1}{x}e^{-\frac{1}{x}}\text{Ei}(x)$$ (assuming $\int_{-a}^a \frac{1}{t} dt=0)$.
It seems that the power series $\sum_{n=0}^\infty n!x^n$ is truly equal to $\frac{1}{x}e^{-\frac{1}{x}}\text{Ei}(x)$, in an analytic continuation kind of sense. However, the series has a radius of convergence equal to $0$, so it shouldn't be possible to fit any function to it. What is going on here? Is it possible to find another function that fits this power series (that doesn't just add a term with power series zero like $e^{-x^{-2}}$)?