A high order schema for the numerical solution of the fractional ordinary differential equations

作者：Cao Junying; Xu Chuanju^*

来源：Journal of Computational Physics, 2013, 238: 154-168.

DOI：10.1016/j.jcp.2012.12.013

• Fractional ordinary differential equations
• High order schema
• Stability and convergence analysis

In this paper we present a general technique to construct high order schemes for the numerical solution of the fractional ordinary differential equations (FODEs). This technique is based on the so-called block-by-block approach, which is a common method for the integral equations. In our approach, the classical block-by-block approach is improved in order to avoiding the coupling of the unknown solutions at each block step with an exception in the first two steps, while preserving the good stability property of the block-by-block schemes. By using this new approach, we are able to construct a high order schema for FODEs of the order alpha, alpha > 0. The stability and convergence of the schema is rigorously established. We prove that the numerical solution converges to the exact solution with order 3 + alpha for 0 < alpha <= 1, and order 4 for alpha > 1. A series of numerical examples are provided to support the theoretical claims.

• 出版日期2013-4-1
• 单位厦门大学