# Integrals of two Bessel functions of the first kind and a modified bessel function of the second kind

I'm searching for a suitable (hopefully simple enough) solution to the following form of integral:

$$\int_0^\infty \mathrm{d}x~x^n J_\nu(a x) J_\nu(b x) K_\mu(c x)$$

Where $n$, $\nu$, and $\mu$ are all integers, and $a$, $b$, and $c$ are all real and positive.

If not generally, a specific case would be quite helpful:

$$\int_0^\infty \mathrm{d}x~x^2 J_\nu(a x) J_\nu(b x) K_1(c x)$$

I am aware of the following:

• Gradshteyn & Ryzhik eq. (6.522.3), which calculates the following integral as a relatively simple function:

$$\int_0^\infty \mathrm{d}x~x K_0(ax) J_\nu(b x) J_\nu(c x)$$

$$\pi \mathrm{i} J_\nu(z) = \mathrm{e}^{-\nu\pi\mathrm{i}/2} K_\nu(-\mathrm{i} z) - \mathrm{e}^{\nu\pi\mathrm{i}/2} K_\nu(\mathrm{i} z), ~~~~|~\mathrm{arg}~z~| \leq \pi/2$$

I simply fail to see how this formula results in a simple change in his equations, especially because of the complex nature of the arguments introduced by the above formula, and the involved definition of $K_\mu$ for integer order.

Should I instead be looking for specific application of G&R eq. 6.522.17-18, which could provide the required formulas, or are there better approaches to this problem? It seems that Fabrikant in the above-linked article says the validity of these formulas is more strictly bounded than shown in G&R.

Any help here would be much appreciated.

• If this better suits math.stackexchange.com, feel free to migrate. I couldn't decide, so I took after this question and put it here. – rubenvb Apr 25 '14 at 8:29
• For the specific case: apply differentiation $\frac{\partial}{\partial a}$ on the formula that you cited from G&R and use $\frac{d}{dz}K_{0}(z)=-K_{1}(z)$ . – Johannes Trost Apr 28 '14 at 15:56
• @Johannes: Yes, indeed, that hit me after writing up this question. The property $\partial_A A^\mu K_\mu(A r) = - r A^\mu K_{\mu-1}(A r)$ combined with $K_{-\mu}(x) = K_\mu(x)$ works wonders for the specific integrals I'm coming across. I'm still interested in how the relation between $J_\nu$ and $K_\nu$ transforms the integral in the Fabrikant paper. The result seems almost trivial. – rubenvb Apr 29 '14 at 8:33

Starting from G&R eq. (6.522.3), $$\int_0^\infty \mathrm{d} r~ r K_0(a r) J_\nu(b r) J_\nu(c r) = \left( \frac{r_2 - r_1}{r_2 + r_1} \right)^{|\nu|} \frac{1}{r_1 r_2}$$ where $r_1^2 = a^2 + (b-c)^2$ and $r_2^2 = a^2 + (b+c)^2$, one can use the derivative and symmetry properties of the bessel-$K$ functions to change the order from 0, incrementing the power of $r$ in the integrand at the same time. Simply calculating the derivative $c~\partial_c$ of the above righthand side gives my example from the question.