Look like we have bessel functions of the first kind (pt.j0 and pt.j1), which you could use to build the second kind (check the article on Bessel functions).
Some of the math heavy hitters might have better insight into whether this might be more complex (har har) than it seems – @maresb @aseyboldt ?