Unhelpful interpolation for solution from differential equation

[cwittens]

change

return solve(prob, Tsit5(), callback=_cb_extinction)

to

return solve(prob, Tsit5(), callback=_cb_extinction, reltol=1e-6, abstol=1e-6)

or even lower/stricter tolerances such as 1e-9 or 1e-12.