Error estimation when using solve_ivp

 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

_______________________________________________
SciPy-User mailing list
[email protected]
https://mail.python.org/mailman/listinfo/scipy-user