Calculating Inverse Laplace Transform of ((1+s/wh)/(1+s/wb))^n I try to compute the inverse Laplace Transform of $\left( \frac{1+s/\omega_{h}}{1+s/\omega_{b}}\right) ^{n}$ with $0<\omega_{b}<\omega_{h}$ and $0<n<1$.
As this function has 2 branch points, one at $-\omega_{b}$ and the other at $-\omega_{h}$, I considered the following integration path:considered path to apply residue theorem.
I found that all the residue are null (perhaps wrong) and there is a contribution only of path $\gamma_{3}$ and $\gamma_{4}$. Thus the result obtained is
$\frac{sin(n\pi)}{\pi}\int_{\omega_{b}}^{\omega_{h}} \left( \frac{1-u/\omega_{h}}{u/\omega_{b}-1}\right) ^{n}e^{-tu} du$
but this is wrong, something is missing in this result, somethink like a dirac function multiplied by $\frac{\omega_{b}}{\omega_{h}}$ that I do not know how to make appear.
Thank you in advance to those who could help me.

First, pull a factor of $(\omega_b/\omega_h)^n$ outside and consider the integral
$$\oint_C dz \, (z+\omega_h)^n (z+\omega_b)^{-n} e^{z t} $$
where $C$ is the modified Bromwich integral that detours about the negative axis around the branch points at $z=-\omega_b$ and $z=-\omega_h$ when $t \gt 0$.  We may then neglect the integrals about the large circular arcs of radius $R$ which vanish as $R \to \infty$ and the small circular arcs about the branch points of radius $\epsilon$ as $\epsilon \to 0$.  The contour integral is then
$$\int_{c-i \infty}^{c+i \infty} ds \, (s+\omega_h)^n (s+\omega_b)^{-n} e^{s t} + e^{i \pi} \int_{\infty}^{\omega_b} dx \, (x-\omega_h)^n (x-\omega_b)^{-n} e^{-x t} \\ + e^{i \pi} \int_{\omega_h}^{\omega_b} (\omega_h-x)^n e^{-i \pi n} (x-\omega_b)^{-n} e^{-x t} + e^{-i \pi} \int_{\omega_b}^{\omega_h} (\omega_h-x)^n e^{i \pi n} (x-\omega_b)^{-n} e^{-x t} \\ + e^{-i \pi} \int_{\omega_b}^{\infty} dx \, (x-\omega_h)^n (x-\omega_b)^{-n} e^{-x t} $$
where $c \gt -\omega_b$.  Note that the factor of $-1$ was replaced with $e^{i \pi}$ above the branch cut $\operatorname{Re}{z} < -\omega_b$ and $e^{-i \pi}$ below the branch cut.  The phase factors in the second and fifth integrals vanish, so that the second and fifth integrals cancel.
By Cauchy's theorem, the contour integral is zero.  We may then write down the inverse Laplace transform in terms of a real integral as follows:
$$\begin{eqnarray} \frac1{i 2 \pi} \int_{c-i \infty}^{c+i \infty} ds \, (s+\omega_h)^n (s+\omega_b)^{-n} e^{s t} &= \frac{\sin{\pi n}}{\pi} \int_{\omega_b}^{\omega_h} (\omega_h-x)^n (x-\omega_b)^{-n} e^{-x t} \\ &= \frac{\sin{\pi n}}{\pi} e^{-\omega_b t} \int_0^{\omega_h-\omega_b} dx \, x^n (\omega_h-\omega_b-x)^{-n} e^{-x t} \\ &= \frac{\sin{\pi n}}{\pi} e^{-\omega_b t} (\omega_h-\omega_b) \int_0^1 du \, u^n (1-u)^{-n} e^{-(\omega_h-\omega_b) t u} \end{eqnarray} $$
The latter integral is a confluent hypergeometric function:
$$\int_0^1 du \, u^n (1-u)^{-n} e^{- \beta u} = \frac{\pi n}{\sin{\pi n}} {}_1F_1(1+n,2,-\beta) $$
Thus, the ILT of the specified function is, for $t \gt 0$,

$$\frac1{i 2 \pi} \int_{c-i \infty}^{c+i \infty} ds \, \left (\frac{1+s/\omega_h}{1+s/\omega_b} \right )^n e^{s t} = n (\omega_h-\omega_b) \left (\frac{\omega_b}{\omega_h} \right )^n  e^{-\omega_b t} {}_1F_1(1+n,2,-(\omega_h-\omega_b) t)$$

For $t \lt 0$, it should be clear that the contour $C$ is closed to the right.  Because the contour encloses no singularities, the ILT is zero in this case.
