[SciPy-User] Error estimation when using solve_ivp

[Lingyi Hu]

 Hi,

I would like to get an error estimate on my results from
`scipy.integrate.solve_ivp`. I am using it to solve a few differential
equations, and using the its root-finding capabilities to get the
coordinates when an event happens. This is all good, but I'm not sure how
to put an error estimate on my result.

I have an idea of how to estimate the global error of a Runge Kutta method
with a fixed step size, but it seems `solve_ivp` uses adaptive step sizes.
It would be great if someone who has experience with this could point me in
the right direction. I don't exactly need a highly rigorous estimation of
the error, just a rough upper bound will do.

Thanks!
Lingyi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-user/attachments/20180401/2d2bc36e/attachment.html>