Abstract
In this paper, the predictor-corrector approach is used to propose two algorithms for the numerical solution of linear and non-linear fractional differential equations (FDE). The fractional order derivative is taken to be in the sense of Caputo and its properties are used to transform FDE into a Volterra-type integral equation. Simpson's 3/8 rule is used to develop new numerical schemes to obtain the approximate solution of the integral equation associated with the given FDE. The error and stability analysis for the two methods are presented. The proposed methods are compared with the ones available in the literature. Numerical simulation is performed to demonstrate the validity and applicability of both the proposed techniques. As an application, the problem of dynamics of the new fractional order non-linear chaotic system introduced by Bhalekar and Daftardar-Gejji is investigated by means of the obtained numerical algorithms.