Best Solver for Bessel-like Differential Equations with Singularity; unexpected Interpolation Behavior

Yeah, I’m trying out everything, looking for the best results.

Thanks for all the advice!