Bien, estamos mirando a la inversa de la transformada de Laplace de:
$$\text{y}_{\space\text{n}}\left(t\right):=\mathscr{L}_\text{s}^{-1}\left[\exp\left(-\text{n}\cdot\text{s}\right)\cdot\mathcal{I}_{\space0}\left(\text{n}\cdot\text{s}\right)\right]_{\left(t\right)}\tag1$$
Donde $\mathcal{I}_{\space\text{p}}\left(\text{z}\right)$ es la función modificada de Bessel de primera especie.
Ahora, usando el teorema de convolución, de la transformada de Laplace, podemos escribir:
$$\text{y}_{\space\text{n}}\left(t\right)=\int_0^t\mathscr{L}_\text{s}^{-1}\left[\mathcal{I}_{\space0}\left(\text{n}\cdot\text{s}\right)\right]_{\left(\tau\right)}\cdot\mathscr{L}_\text{s}^{-1}\left[\exp\left(-\text{n}\cdot\text{s}\right)\right]_{\left(t-\tau\right)}\space\text{d}\tau\tag2$$
El uso de la tabla seleccionada de las transformadas de Laplace, podemos escribir:
$$\mathscr{L}_\text{s}^{-1}\left[\exp\left(-\text{n}\cdot\text{s}\right)\right]_{\left(t-\tau\right)}=\delta\left(t-\tau-\text{n}\right)\tag3$$
Donde $\delta\left(x\right)$ es la función Delta de Dirac.
Así, podemos reescribir la ecuación de $\left(2\right)$ como sigue:
$$\text{y}_{\space\text{n}}\left(t\right)=\int_0^t\mathscr{L}_\text{s}^{-1}\left[\mathcal{I}_{\space0}\left(\text{n}\cdot\text{s}\right)\right]_{\left(\tau\right)}\cdot\delta\left(t-\tau-\text{n}\right)\space\text{d}\tau\tag4$$
Utilizando la definición de la función modificada de Bessel de primera especie, podemos escribir:
$$\mathcal{I}_{\space0}\left(\text{n}\cdot\text{s}\right)=\sum_{\text{k}=0}^\infty\frac{1}{\text{k}!}\cdot\frac{1}{\Gamma\left(1+\text{k}\right)}\cdot\left(\frac{\text{n}\cdot\text{s}}{2}\right)^{2\text{k}}\tag5$$
Donde $\Gamma\left(\text{s}\right)$ es la función Gamma.
Así que el uso de la tabla seleccionada de las transformadas de Laplace, podemos escribir:
$$\mathscr{L}_\text{s}^{-1}\left[\mathcal{I}_{\space0}\left(\text{n}\cdot\text{s}\right)\right]_{\left(\tau\right)}=\sum_{\text{k}=0}^\infty\frac{1}{\text{k}!}\cdot\frac{1}{\Gamma\left(1+\text{k}\right)}\cdot\mathscr{L}_\text{s}^{-1}\left[\left(\frac{\text{n}\cdot\text{s}}{2}\right)^{2\text{k}}\right]_{\left(\tau\right)}=$$
$$\sum_{\text{k}=0}^\infty\frac{1}{\text{k}!}\cdot\frac{1}{\Gamma\left(1+\text{k}\right)}\cdot\frac{1}{\Gamma\left(-2\text{k}\right)}\cdot\left(\frac{\text{n}}{2}\right)^{2\text{k}}\cdot\frac{1}{\tau^{1+2\text{k}}}\tag6$$
Ahora, podemos reescribir la ecuación de $\left(4\right)$ como sigue:
$$\text{y}_{\space\text{n}}\left(t\right)=\int_0^t\left\{\sum_{\text{k}=0}^\infty\frac{1}{\text{k}!}\cdot\frac{1}{\Gamma\left(1+\text{k}\right)}\cdot\frac{1}{\Gamma\left(-2\text{k}\right)}\cdot\left(\frac{\text{n}}{2}\right)^{2\text{k}}\cdot\frac{1}{\tau^{1+2\text{k}}}\right\}\cdot\delta\left(t-\tau-\text{n}\right)\space\text{d}\tau=$$
$$\sum_{\text{k}=0}^\infty\frac{1}{\text{k}!}\cdot\frac{1}{\Gamma\left(1+\text{k}\right)}\cdot\frac{1}{\Gamma\left(-2\text{k}\right)}\cdot\left(\frac{\text{n}}{2}\right)^{2\text{k}}\int_0^t\frac{\delta\left(t-\tau-\text{n}\right)}{\tau^{1+2\text{k}}}\space\text{d}\tau\tag7$$
Al $\text{n}>0\space\wedge\space\text{k}\ge0\space\wedge\space t\ge0$, podemos escribir:
$$\int_0^t\frac{\delta\left(t-\tau-\text{n}\right)}{\tau^{1+2\text{k}}}\space\text{d}\tau=\frac{2\cdot\theta\left(t\right)-1}{\left(t-\text{n}\right)^{1+2\text{k}}}\cdot\theta\left(t-\text{n}-t\cdot\theta\left(-t\right)\right)\cdot\theta\left(\text{n}-t+t\cdot\theta\left(t\right)\right)\tag8$$
Donde $\theta\left(\text{a}\right)$ es la función escalón unitario.