  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.
