I'd like to calculate $$ \int_0^R r^2 \left[ I_0(r)K_1(r) + I_1(r)K_0(r) \right] \mathrm{d}r $$
Earlier, I have used the property $$\int x I_0(x) \,\mathrm{d}x = x I_1(x)$$ but I can't find a similar standard integral in books and tables for my problem. From another source I am hoping the answer to be like $\frac{R^2}{2}$, but Mathematica gives a complicated answer involving Meijer G-functions. Integrating by parts doesn't seem to help.