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

Specific Domains Numerics
2

Try:

RKO65 - Olver’s 6 stage 5th order method. This method is robust on problems which have a singularity at t=0 .

https://diffeq.sciml.ai/latest/solvers/ode_solve/

If that doesn’t work for your problem, you might want to try ApproxFun. Most of the methods in DiffEq are not safe with a singularity at t=0: generally that might need a discretization that’s fully implicit over time.

2 Likes