A Mathematica package for the solution of IVPs