Abstract
We propose a numerical verification method for initial value problems(IVP) of ODEs in order to enclose the solutions at the end time by considerably small bounds. The method is based on Nakao's theory which is established for numerical verification methods for PDEs. We construct boundary value problems(BVP) of ODEs whose solutions can be used to define the bounds of the solutions to the IVP. Then Nakao's theory is applied to the BVPs. Numerical examples are shown to compare our method with existing validated computation methods for ODEs.
