Edit. Primera versión carecía de claridad y había algunos errores.
Primera nota de que
\begin{align}
I_n(t) =\frac{1} {i2} \int^\infty_0 x^{n-1}e^{(i-t) x} \,dx - \frac{1} {i2} \int^\infty_0 x^{n-1}e^{(-i-t) x} \,dx
\end{align}
Tenemos dos integrales. Ambos pueden estar relacionados con una integral de contorno. Uno también puede ver de alguna manera el $\Gamma$-función en ella, de modo que sería muy bueno expres estas integrales en términos de $\Gamma$-función. Que es lo que vamos a hacer y le muestra sólo para la primera integral. Para que considere la posibilidad de:
\begin{align}
\tag{1}
\oint_C z^{n-1}e^z\, dz
\end{align}
El contorno de $C$ es un círculo sector con ángulo de $\arctan(1/t)$ que va desde el origen de todo el camino hasta en el segundo cuadrante y regresar todo el camino hacia abajo a través de la negativa del eje real. Deje que el radio del "círculo" ser $R$. El contorno puede ser visto a continuación.
$\hspace{70pt} $
Al $R$ va al infinito la circular contribución va a cero por el ML-lema. Puesto que el integrando es holomorphic, por Cauchy teorema tenemos que la integral en $(1)$ es cero y por parametrización (y $R\to \infty $) :
\begin{align}
\int^0_{-\infty}x^{n-1}e^x\, dx+(i-t)^n\int^\infty_0 x^{n-1}e^{(i-t) x} \, dx=0
\end{align}
Sustituimos $u=-x$ para la primera integral:
\begin{align}
\int^\infty_{0}(-1)^{n-1}u^{n-1} e^{-u}\,du + (i-t)^n\int^\infty_0 x^{n-1}e^{(i-t) x} \, dx=0
\end{align}
Damas y caballeros, eso es lo que queríamos: la primera integral es el $\Gamma$-función y todo lo que ello implica:
\begin{align}
\frac{1} {i2} \int^\infty_0 x^{n-1}e^{(i-t) x} \,dx=\frac{(n-1)!(-1)^{n}} {i2(i-t)^n} =\frac{(n-1)!(t+i)^n}{i2(1+t^2)^n}
\end{align}
Del mismo modo, se obtiene:
\begin{align}
\frac{1} {i2} \int^\infty_0 x^{n-1}e^{(-i-t) x} \,dx=\frac{(n-1)!(t-i)^n}{i2(1+t^2)^n}
\end{align}
La adición de ellos se obtiene el resultado deseado.