A numerical approach for solving systems of nonautonomous ordinary differential equations (ODEs) is proposed under suitable assumptions. This approach is based on expansion of the solutions of ODEs by Volterra series and allows to estimate the distance between the obtained approximation and the true trajectory.