As of 1 April 2026, the PSJD database will become an archive and will no longer accept new data. Current publications from Polish scientific journals are available through the Library of Science: https://bibliotekanauki.pl
We 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.