Let $f(\epsilon)=\int_0^{\infty}\frac{e^{-t}}{(1+\epsilon t)^{1/2}}dt$. Derive an asymptotic expansion of the form $f(\epsilon)\approx \sum_{c=0}^{\infty}c_n\epsilon^n$ for $\epsilon\to 0+$.
Idea: This looks a bit like Watson's Lemma, which says that if $f(t)$ is continuous on $(0,\infty)$ with $f(t)=O(e^{\mu t})$ for some $\mu>0$ as $t\to\infty$ and have an asymptotic expansion of the form $f(t)\approx\sum_{n=0}^{\infty}c_nt^{\lambda_n}$ for $\lambda_n>-1$ a sequence of strictly increasing real numbers. Then \begin{equation}\int_0^{\infty}e^{-xt}f(t)dt\approx\sum_{n=0}^{\infty}c_n\frac{\Gamma(\lambda_n+1)}{x^{\lambda_n+1}}\end{equation} as $x\to\infty$. But I am not sure how to transform the original integral to the required form - we have an $e^{-t}$ term but not $e^{-xt}$. So we could, for instance, make the substitution $t=xs$, but I don't see how to proceed here.