Some ordinary differential equations do not have exact solutions. Their solutions can be approximated by numerical methods. This thesis presents several numerical methods for solving IVPs. Moreover, the methods are compared in terms of their accuracy, convergence, and consistency. Among all the selected methods, the four-stage Runge-Kutta method of fourth-order is verified to be the most effective method as it could be expected.