Daniel Wilczak
Schedae Informaticae, Volume 20, 2011, pp. 9-42
https://doi.org/10.4467/20838476SI.11.001.0287We present a Lohner type algorithm for the computation of rigorous bounds for the solutions of ordinary differential equations and its derivatives with respect to the initial conditions up to an arbitrary order.