Closed form for this integral $\int_0^\infty \frac{e^{-x} dx}{\sqrt{(x+a)^2+b^2}}$ This is a Laplace transform, however I couldn't find it in the tables and Wolfram doesn't know the answer either:
$$I(a,b)=\int_0^\infty \frac{e^{-x} dx}{\sqrt{(x+a)^2+b^2}}$$
Some kind of closed form (even in terms of special functions) will help me very much with my problem, as this is the simple part of a double integral which will need to be evaluated many times.
With some simple substitutions we can write:
$$I(a,b)=e^a \int_a^\infty \frac{e^{-y} dy}{\sqrt{y^2+b^2}}=e^a \int_{a/b}^\infty \frac{e^{-b z} dz}{\sqrt{z^2+1}}=e^a \int_{\sinh^{-1} a/b}^\infty e^{-b \sinh u} du$$
Looks like some kind of incomplete Bessel function to me, but is there a standart closed form expression?

I've been trying to find a possible way to deal with this integral, an here's one promising attempt:
$$I(a,b)=\int_0^\infty \frac{e^{-x} dx}{\sqrt{x^2+2ax+a^2+b^2}}$$
Let's consider two separate cases:
$$1) \qquad x<\sqrt{a^2+b^2}$$
$$r_1=\frac{x}{\sqrt{a^2+b^2}} <1 \\ q= -\frac{a}{\sqrt{a^2+b^2}}, \quad |q| <1$$
$$x^2+2ax+a^2+b^2=(a^2+b^2) (r_1^2-2 q r_1+1)$$
$$2) \qquad x > \sqrt{a^2+b^2}$$
$$r_2=\frac{\sqrt{a^2+b^2}}{x} <1$$
$$x^2+2ax+a^2+b^2=x^2 (r_2^2-2 q r_2+1)$$
Which means:
$$I(a,b)=\int_0^\sqrt{a^2+b^2} \frac{e^{-x} dx}{\sqrt{a^2+b^2}} \sum_{n=0}^\infty P_n (q) r_1^n+\int_\sqrt{a^2+b^2}^\infty  \frac{e^{-x} dx}{x} \sum_{n=0}^\infty P_n (q) r_2^n$$
Where $P_n$ are Legendre polynomials.
Let's rename $$c=\sqrt{a^2+b^2}, \qquad r=r_1$$
$$I(a,b)= \sum_{n=0}^\infty P_n (q) \left( \int_0^1 r^n e^{-c r} dr+ \int_1^\infty  \frac{e^{-c r} dr}{r^{n+1}} \right)=$$
$$=\sum_{n=0}^\infty P_n (q) \left( \frac{n!-\Gamma(n+1,c)}{c^{n+1}}+ \operatorname{Ei}_{n+1} (c) \right)$$
This is rather ugly, so if anyone has a better option, please share.
