and validated numerics to rigorously enclose solutions of initial-value problems. Our
approach uses a combination of contraction mapping arguments based on a Newton–
Kantorovich type theorem and piecewise polynomial interpolation. Completing multiple
steps of integration is challenging, and we resolve it by smooth interpolation of the previous
solution, resulting in an interval-valued polynomial initial condition for the subsequent step …