A novel high-order algorithm for the numerical estimation of fractional differential equations

This paper uses polynomial interpolation to design a novel high-order algorithm for the numerical estimation of fractional differential equations. The Riemann-Liouville fractional derivative is expressed by using the Hadamard finite-part integral and the piecewise cubic interpolation polynomial is utilized to approximate the integral. The detailed error analysis is presented and it is established that the convergence order of the algorithm is O(h4−α). Asymptotic expansion of the error for the presented algorithm is also investigated. Some numerical examples are provided and compared with the exact solution to show that the numerical results are in well agreement with the theoretical ones and also to illustrate the accuracy and efficiency of the proposed algorithm.