AN ALGORITHM FOR THE NUMERICAL SOLUTION OF DIFFERENTIAL EQUATIONS OF FRACTIONAL ORDER

We present and discuss an algorithm for the numerical solution of initial value problems of the form $D_*^\alpha$y(t) = f(t, y(t)), y(0) = y0, where $D_*^\alpha$y is the derivative of y of order $\alpha$ in the sense of Caputo and 0<${\alpha}{\leq}1$. The algorithm is based on the fractional Euler's method which can be seen as a generalization of the classical Euler's method. Numerical examples are given and the results show that the present algorithm is very effective and convenient.