Exponential quadrature rules for linear fractional differential equations

From MaRDI portal
Publication:2018704



Abstract: This paper focuses on the numerical solution of initial value problems for fractional differential equations of linear type. The approach we propose grounds on expressing the solution in terms of some integral weighted by a generalized Mittag-Leffler function. Then suitable quadrature rules are devised and order conditions of algebraic type are derived. Theoretical findings are validated by means of numerical experiments and the effectiveness of the proposed approach is illustrated by means of comparisons with other standard methods.


The authors investigate numerical solutions of linear fractional differential equations of the following form \[ D^{\alpha}_{t_{0}}y(t)+\lambda y(t)=f(t), \] where \( \alpha >0\) is a real number, \( \lambda \in \mathbb R \), \( y:[t_{0},T]\rightarrow \mathbb R \) and the forcing term \( f:[t_{0},T]\rightarrow \mathbb R \) is sufficiently smooth. The symbol \( D^{\alpha}_{t_{0}} \) denotes the derivative operator of non-integer order \( \alpha \), with respect to the origin \( t_{0}\), according to the Caputo definition \[ D^{\alpha}_{t_{0}}y(t)= \frac{1}{\Gamma(m-\alpha)}\int^{t}_{t_{0}}\frac{y^{(m)}(u)}{(t-u)^{\alpha+1-m}}du,\tag{\(*\)} \] with the gamma function \( \Gamma(.) \), \( m=\lceil\alpha\rceil \) the smallest integer such that \( \alpha < m \), \( y^{(m)} \) is the derivative of \( y \) of integer order \( m \) and the initial conditions \( y^{(k)}(t_{0})=y_{0,k}\), \(k=0,1,\dots,m-1 \). In order to obtain high accuracy under suitable and reasonable assumptions of smoothness on \( f(t)\), Eq.(*) is reformulated in terms of some special weighted integrals involving \( f(t) \) alone as the integrand and then specific quadrature rules are formulated for their numerical approximation. Numerical examples are given to illustrate the results obtained. The effectiveness of the present approach is demonstrated by means of comparisons with other standard methods.



Cites work



Describes a project that uses

Uses Software






This page was built for publication: Exponential quadrature rules for linear fractional differential equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2018704)