A quadrature method for numerical solutions of fractional differential equations

In this article, a numerical method is developed to obtain approximate solutions for a certain class of fractional differential equations. The method reduces the underlying differential equation to system of algebraic equations. An algorithm is presented to compute the coefficient matrix for the resulting algebraic system. Several examples with numerical simulations are provided to illustrate effectiveness of the method.