Abstract: This study describes a three-step method for the numerical solution of ordinary differential equations with singularities. The scheme is based on rational functions approximation technique and its development and analysis is based on power series expansions (Taylor and Binomial) and Dalhquist stability test method. The scheme is convergent and A-stable. Numerical results show that the scheme is accurate, effective and efficient.