Implementation of classical numerical analysis algorithms in Wolfram Mathematica for solving ordinary differential equations