Numerical Solutions of Initial Value Problems Using Mathematica