Also, see Best Solver for Bessel-like Differential Equations with Singularity; unexpected Interpolation Behavior - #2 by ChrisRackauckas