The gamma function is an integral of the form $$\int_0^\infty e^{-t} t^n \,dt = n!$$
One can even throw in a coefficient without much trouble $$\int_0^\infty e^{-at} t^n \,dt = \frac{n!}{a^{n+1}}$$
But if we want to decay even faster, like $$\int_0^\infty e^{-e^t} t^n \,dt = ?$$ then a close approximation seems difficult, much less an exact answer. The new integral DOES decay far faster; for n=25 I get around 34,000 in comparison to 25! for the gamma function.
What is the best known approximation to this integral? My end goal is an approximation of $$\int_0^\infty e^{-\lambda e^{t/x}} t^{2n} \,dt = ?$$ where we fix $\lambda,x>0$ and $n\in\mathbb{N}$.