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

**URL:** <https://discourse.julialang.org/t/best-solver-for-bessel-like-differential-equations-with-singularity-unexpected-interpolation-behavior/41482>\
**Category:** Numerics\
**Created:** [June 15, 2020, 11:39pm UTC](https://discourse.julialang.org/t/best-solver-for-bessel-like-differential-equations-with-singularity-unexpected-interpolation-behavior/41482 "2020-06-15T23:39:00Z")\
**Posts on this page:** 1\
**Showing post:** 2

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [June 15, 2020, 11:43pm UTC](https://discourse.julialang.org/t/best-solver-for-bessel-like-differential-equations-with-singularity-unexpected-interpolation-behavior/41482/2 "2020-06-15T23:43:22Z")

</div>

> [@CosmoProf](#):
>
> oscillatory in general, with a singularity in the function and its derivatives at `t=0` .

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/](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.

---

_[View the full topic](https://discourse.julialang.org/t/best-solver-for-bessel-like-differential-equations-with-singularity-unexpected-interpolation-behavior/41482)._