WolframAlpha tells me that $I_0(x)K_1(x)+I_1(x)K_0(x)=\frac{1}{x}$ for any $x$, but I have no idea how to prove it.
Should I use the recurrence relations for Bessel functions? Or should I use the series expansion around 0?
Is this identity widely known? If so, I would like a reference. It greatly simplifies the expression that appears in my physics application.