Abstract
Although there have existed some numerical algorithms for the fractional differential equations, developing high-order methods (i.e., with convergence order greater than or equal to 2) is just the beginning. Lubich has ever proposed the high-order schemes when he studied the fractional linear multistep methods, where he constructed the pth order schemes (p = 2, 3, 4, 5, 6) for the αth order Riemann-Liouville integral and αth order Riemann-Liouville derivative. In this paper, we study such a problem and develop recursion formulas to compute these coefficients in the higher-order schemes. The coefficients of higher-order schemes (p = 7,8, 9,10) are also obtained. We first find that these coefficients are oscillatory, which is similar to Runge's phenomenon. So, they are not suitable for numerical calculations. Finally, several numerical examples are implemented to testify the efficiency of the numerical schemes for p = 3,…, 6.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
