An efficient numerical method for solving systems of fractional ordinary differential equations

Abstract

In this work, we apply an efficient method based on hybrid functions for solving linear and non-linear systems of fractional ordinary differential equations (SFODEs). Here, we consider the fractional derivatives in the Caputo sense. By using the present method, a system of FODEs is reduced to a system of algebraic equations which can be solved by a proper numerical method. In convergence discussion of the method, an upper bound of the error is obtained. To show the efficiency and the accuracy of this method, some examples are simulated and then some comparisons between the outputs with those of several other methods are carried out.