In this paper we present and discuss an algorithm for the numerical solution of some kinds of fractional differential equations, i.e. of order n+v, n ≤ 3, 2 nd order 2v and nv, n ≥ 3. The algorithm for the numerical solution of these equations is based on a Newton approach. The stability and convergence of the fractional order numerical method are described. Finally, some numerical examples are provided to show that the numerical method for solving the fractional differential equation is an effective solution method.