This paper deals with the analytical approximate solution of Abel differential equations (ADEs) of the type du/dt = Σmk=0 gk(t)uk by proposing a new modified version of the parametric iteration method (PIM). The modified algorithm analytically approximates the solution of ADEs in a sequence of subintervals which is continuous everywhere. The local convergence and the stability of the algorithm are discussed in details. Also we show how to characterize the stability function and the region on which the algorithm is presented. Some examples are given to demonstrate the efficiency and accuracy of this algorithm. Comparison with numerical Runge-Kutta methods (RK) shows that the modified algorithm presented in this paper has the advantage of giving an analytical form of the solution, which is not possible in the purely numerical RK techniques. Moreover, the approximations obtained by the new algorithm converge faster than the numerical RK4 methods, as will be shown. The obtained results reveal that the present algorithm is a promising iterative analytical method for solving ADEs. Furthermore, the proposed algorithm provides us with an easy way to modify the rate and region convergence of the solution. Most promising, however, it seems that the newly developed technique can be further implemented easily to solve other nonlinear ordinary differential equations (ODEs) of physical interests.

Keywords