A numerical method for solution of ordinary differential equations of fractional order

In this paper we propose an algorithm for the numerical solution of arbitrary differential equations of fractional order. The algorithm is obtained by using the following decomposition of the differential equation into a system of differential equation of integer order connected with inverse forms of Abel-integral equations. The algorithm is used for solution of the linear and non-linear equations.