Inverse Fourier transform of modified Bessel function I'm relatively new to Fourier transforms so apologize in advance if this problem seems trivial.
In order to solve a second order PDE I have defined the following sine Fourier transform
$$V(r,\lambda)=\sqrt{\frac{2}{\pi}}\int_{0}^{\infty}v(r,z)\sin(\lambda z)\,\textrm{d}z.$$
By doing so I arrive at the following solution for $V$
$$V(r,\lambda)=\sqrt{\frac{2}{\pi}}\frac{1}{\lambda}\frac{I_{1}(\lambda r)}{I_{1}(\lambda)},$$
where $I_{\alpha}$ is the modified Bessel function of the first kind.
My goal is to invert the Fourier transform and obtain the solution for $v$. So far I have that
\begin{align*}
v(r,z)&=\sqrt{\frac{2}{\pi}}\int_{0}^{\infty}V(r,\lambda)\sin(\lambda z)\,\textrm{d}\lambda
\\
&=\frac{2}{\pi}\int_{0}^{\infty}\frac{\sin(\lambda z)}{\lambda}\frac{I_{1}(\lambda r)}{I_{1}(\lambda)}\,\textrm{d}\lambda.
\end{align*}
I know that the definition of the modified Bessel function of the first kind gives
$$I_{1}(\lambda r)=\sum_{m=0}^{\infty}\frac{1}{m!\Gamma(m+2)}\left(\frac{\lambda r}{2}\right)^{2m+1}.$$
Therefore
\begin{align*}
v(r,z)&=\sqrt{\frac{2}{\pi}}\int_{0}^{\infty}V(r,\lambda)\sin(\lambda z)\,\textrm{d}\lambda
\\
&=\frac{2}{\pi}\int_{0}^{\infty}\frac{\sin(\lambda z)}{\lambda}\frac{\sum_{m=0}^{\infty}\frac{1}{m!\Gamma(m+2)}\left(\frac{\lambda r}{2}\right)^{2m+1}}{\sum_{m=0}^{\infty}\frac{1}{m!\Gamma(m+2)}\left(\frac{\lambda}{2}\right)^{2m+1}}\,\textrm{d}\lambda.
\end{align*}
This looks horrible! Can anyone help me from here on out?
If it's any help I know that the solution should be 
$$v(r,z)=2\sum_{\beta_{n}}\frac{1}{\beta_{n}}\frac{J_{1}(r\beta_{n})}{J_{0}(\beta_{n})}e^{-\beta_{n}z},$$
where $J_{\alpha}$ is the Bessel function of the first kind and $\beta_{n}$ are the zeros of $J_{1}(\beta)$.
Thanks!

A Fourier-Bessel representation for the ratio 
\begin{equation}
\frac{I_1(\lambda r)}{I_1(\lambda)}=\frac{J_1(i\lambda r)}{J_1(i\lambda)}
\end{equation} 
can be found in Ederlyi, Higher Transcendental functions II, 7.10.4 (56):
\begin{equation}
\frac{I_1(\lambda r)}{I_1(\lambda)}=2\sum_{m=1}^\infty\frac{\beta_mJ_1(\beta_m r)}{\left( \beta_m^2+\lambda^2 \right)J_2(\beta_m)}
\end{equation} 
Then,
\begin{equation}
v(r,z)=\frac{4}{\pi}\sum_{m=1}^\infty\frac{\beta_mJ_1(\beta_m r)}{J_2(\beta_m)}\int_{0}^{\infty}\frac{\sin(\lambda z)}{\lambda\left( \beta_m^2+\lambda^2 \right)}\,\textrm{d}\lambda
\end{equation} 
This integral can be evaluated easily
\begin{equation}
\int_{0}^{\infty}\frac{\sin(\lambda z)}{\lambda\left( \beta_m^2+\lambda^2 \right)}\,\textrm{d}\lambda=\frac{\pi}{2\beta_m^2}\left( 1-e^{-\beta_mz} \right)
\end{equation} 
thus
\begin{align}
v(r,z)&=-2\sum_{m=1}^\infty\frac{J_1(\beta_m r)\left( 1-e^{-\beta_mz} \right)}{\beta_mJ_0(\beta_m)}\\
&=2\sum_{m=1}^\infty\frac{J_1(\beta_m r)e^{-\beta_mz}}{\beta_mJ_0(\beta_m)}+r
\end{align} 
where we have used $J_0(\beta_m)=-J_2(\beta_m)$ and the Fourier Bessel representation of $r$, which can be found in Ederlyi7.10.4 (54), for example. Obtained result is slightly different from the one you quote, though both forms are very similar. Notice however that the present expression is zero when $z=0$, as expected from the integral expression.
