Abstract: This paper investigated extrapolation methods for systems of ordinary differential equations and presented a high accuracy algorithm,with the following characteristic features:(1)extrapolated by continued fractions;(2)based on a cheaper smoothing scheme;(3)adopted an efficient error criteria,It is demonstrated that the new algorithm is more efficient than the polynomial extrapolation algorithm and GBS algorithm.