Understanding abstol and reltol in DifferentialEquations.jl

Are you modifying u in f? What is plotted doesn’t look natural, it looks like it’s violating something. And the instabilities around low tolerances is a hallmark that you’re doing something that violates step rejection.

See