In recent years much attention has been given to the numerical solution of ODEs. Of particular interest has been the solution of singularly perturbed and stiff problems. These types of problems arise in various fields of science and engineering such as fluid mechanics, physics, chemistry, mechanics, chemical reactor theory, convection diffusion processes, optimal control and other branches of applied mathematics. Singular perturbation problems depend on the presence of a small, positive parameter which provides a multi-scale character to the solution. That is the solution varies very rapidly...