I am trying to prove that the integral \begin{align} \int_{0}^{\infty } e^{-\frac{r^2}{2B}} r^{l-n} L_n^{l-n}\left(\frac{r^2}{C}\right) I_{l-n}\left(\rho r \right) r dr \end{align}

has the form \begin{align} B^{l+1} e^{\frac{B}{2}\rho^2} \rho^{l-n} L_n^{l-n}\left(\frac{\rho^2}{C} \right), \end{align} where $L_n^{l-n}$ is the generalised Laguerre function, and $I_{l-n}$ is the modified Bessel function of the first kind.

Expanding the modified Bessel function into an infinite sum, and then using Eq. 7, section 7.414 (pg. 809) from Tables of Integrals, Series & Products (Ed. 7) (by I.A. Gradshteyn & I.M.Ryzhik), which is \begin{align} \int_{0}^{\infty} e^{-st} t^{\beta} L_n^{\alpha}(t) dt = \frac{\Gamma(\beta+1) \Gamma(\alpha+n+1)} {n! \Gamma(\alpha+1) s^{\beta+1}} {}_1F_2\left(-n, \beta+1; \: \alpha+1; \: \frac{1}{s} \right), \end{align} I can get close, but not close enough!

Does anyone know how to do this?

Please let me know if you need any further information. I haven't done this very often, but the few times I have I always seem to forget something pertinent!