I want to find the full asymptotic behaviour of the integral
$$I = \int_1^\infty (t^2-1)^{-1/2}e^{-xt}\,{\rm d}t$$
at $x\to\infty$. So firstly we note that the integrand doesn't exist at $t=1$ so we must attempt to remove the "singularity". Letting $t=\cosh{(u)}$ gives
$$I = \int_0^\infty e^{-x\cosh{u}}\,{\rm d}u$$
then recall that
$$e^x = \sum_{n=0}^\infty \frac{x^n}{n!}\qquad \cosh{x}=\sum_{m=0}^\infty \frac{x^{2m}}{(2m)!}$$
which both converge for all $x\in\mathbb{R}$. Using these we find that
$$I=\int_0^\infty \sum_{n=0}^\infty\frac{(-x)^n}{n!}\left(\sum_{m=0}^\infty \frac{u^{2m}}{(2m)!}\right)^n\,{\rm d}x$$
Now how would I progress?
Sub $t=\cosh{u}$; then
$$I(x) = \int_0^{\infty} du \, e^{-x \cosh{u}} $$
Apply Laplace's method. Note that, while
$$\cosh{u} = 1 + \frac{u^2}{2!} + \sum_{k=2}^{\infty} \frac{u^{2 k}}{(2 k)!} $$
we Taylor expand the exponential only for the higher order terms in the sum in the limit as $x \to \infty$. Thus, for example, to order $2 N$:
$$I(x) = e^{-x} \int_0^{\infty} du \, e^{-x u^2/2} \left [1+ a_4 u^4 + a_6 u^6 + \cdots + O(u^{2 N}) \right ] $$
Now, the coefficients $a_4$, $a_6$, and so on would need to be computed via Taylor expansion of the exponential of the above sum. Use the fact that
$$\int_0^{\infty} du \, u^{2 k} e^{-x u^2/2} = \frac12 \left (\frac{2}{x} \right )^{k+1/2} \Gamma \left ( k+\frac12 \right ) $$
and you can derive the full asymptotic expansion to as many terms as you need.