A Second Order Approximation for the Caputo Fractional Derivative

From MaRDI portal
Publication:3384994

zbMATH Open1488.26013arXiv1502.00719MaRDI QIDQ3384994FDOQ3384994

Yuri Dimitrov

Publication date: 17 December 2021

Abstract: When 0<alpha<1, the approximation for the Caputo derivative y^{(alpha)}(x) = frac{1}{Gamma(2-alpha)h^alpha}sum_{k=0}^n sigma_k^{(alpha)} y(x-kh)+O�igl(h^{2-alpha}�igr), where sigma0(alpha)=1,sigman(alpha)=(n1)1an1a and sigma_k^{(alpha)} = (k-1)^{1-alpha}-2k^{1-a}+(k+1)^{1-alpha},quad (k=1...,n-1), has accuracy . We use the expansion of sumk=0nkalpha to determine an approximation for the fractional integral of order 2alpha and the second order approximation for the Caputo derivative y^{(alpha)}(x) = frac{1}{Gamma(2-alpha)h^alpha}sum_{k=0}^n delta_k^{(alpha)} y(x-kh)+O�igl(h^{2}�igr), where deltak(alpha)=sigmak(alpha) for 2leqkleqn, delta_0^{(alpha)} = sigma_0^{(alpha)}-zeta(alpha-1), delta_1^{(alpha)} = sigma_1^{(alpha)}+2zeta(alpha-1),delta_2^{(alpha)} = sigma_2^{(alpha)}-zeta(alpha-1), and zeta(s) is the Riemann zeta function. The numerical solutions of the fractional relaxation and subdiffusion equations are computed.


Full work available at URL: https://arxiv.org/abs/1502.00719




Recommendations





Cited In (11)





This page was built for publication: A Second Order Approximation for the Caputo Fractional Derivative

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