Numerical Methods for Ordinary Differential Systems: The Initial Value Problem