Search

In this paper, fractional differential equations in the sense of Caputo-Fabrizio derivative are transformed into integral equations. Then a high order numerical method for the integral equation is investigated by approximating the integrand with a piece- wise quadratic interpolant. The scheme is capable of handling both linear and nonlinear fractional differential equations. A detailed error analysis and stability region of the numerical scheme is rigorously established. [ABSTRACT FROM AUTHOR]

In the present work, first of all, a new numerical fractional differentiation formula (called the CF2 formula) to approximate the Caputo-Fabrizio fractional derivative of order a; (0 <α < 1) is developed. It is established by means of the quadratic interpolation approximation using three points (t j-2;y(t j-2)); (t j-1;y(t j-1)), and (t j; y(t j)) on each interval [t j-1; t j] for (j ≥ 2), while the linear interpolation approximation are applied on the first interval [t0; t1]. As a result, the new formula can be formally viewed as a modification of the classical CF1 formula, which is obtained by the piecewise linear approximation for y(t). Both the computational efficiency and numerical accuracy of the new formula is superior to that of the CF1 formula. The coefficients and truncation errors of this formula are discussed in detail. Two test examples show the numerical accuracy of the CF2 formula. The CF1 formula demonstrates that the new CF2 is much more effective and more accurate than the CF1 when solving fractional differential equations. Detailed stability analysis and region stability of the CF2 are also carefully investigated. [ABSTRACT FROM AUTHOR]