Differential equations of fractional order are believed to be more challenging to compute compared to the integer-order differential equations due to its arbitrary properties. This study proposes a multistep method to solve fractional differential equations. The method is derived based on the concept of a third-order Adam–Bashforth numerical scheme by implementing Lagrange interpolation for fractional case, where the fractional derivatives are defined in the Caputo sense. Furthermore, the study includes a discussion on stability and convergence analysis of the method. Several numerical examples are also provided in order to validate the reliability and efficiency of the proposed method. The examples in this study cover solving linear and nonlinear fractional differential equations for the case of both single order as α∈(0,1) and higher order, α∈1,2, where α denotes the order of fractional derivatives of Dαy(t). The comparison in terms of accuracy between the proposed method and other existing methods demonstrate that the proposed method gives competitive performance as the existing methods.

Fractional differential equations have recently demonstrated their importance in a variety of fields, including medicine, applied sciences, and engineering. The main objective of this study is to propose an Adams-type multistep method for solving differential equations of fractional order. The method is developed by implementing the Lagrange interpolation and taking into account the idea of the Adams–Moulton method for fractional case. The fractional derivative applied in this study is in the Caputo derivative operator. The analysis of the proposed method is presented in terms of order of the method, order of accuracy, and convergence analysis, with the proposed method being proved to converge. The stability of the method is also examined, where the stability regions appear to be symmetric to the real axis for various values of α. In order to validate the competency of the proposed method, several numerical examples for solving linear and nonlinear fractional differential equations are included. The method will be presented in the numerical predict–correct technique for the condition where $\alpha \in (0,1)$ α ∈ ( 0 , 1 ) , in which α represents the order of fractional derivatives of $D^{\alpha }y(t)$ D α y ( t ) .