Abstract:The asymptotic stable region and long-time decay rate of solutions to linear homogeneous Caputo time fractional ordinary differential equations (F-ODEs) are known to be completely determined by the eigenvalues of the coefficient matrix. Very different from the exponential decay of solutions to classical ODEs, solutions of F-ODEs decay only polynomially, leading to the so-called Mittag-Leffler stability, which was already extended to semi-linear F-ODEs with small perturbations. This work is mainly devoted to the qualitative analysis of the long-time behavior of numerical solutions. By applying the singularity analysis of generating functions developed by Flajolet and Odlyzko (SIAM J. Disc. Math. 3 (1990), 216-240), we are able to prove that both $\mathcal{L}$1 scheme and strong $A$-stable fractional linear multistep methods (F-LMMs) can preserve the numerical Mittag-Leffler stability for linear homogeneous F-ODEs exactly as in the continuous case. Through an improved estimate of the discrete fractional resolvent operator, we show that strong $A$-stable F-LMMs are also Mittag-Leffler stable for semi-linear F-ODEs under small perturbations. For the numerical schemes based on $\alpha$-difference approximation to Caputo derivative, we establish the Mittag-Leffler stability for semi-linear problems by making use of properties of the Poisson transformation and the decay rate of the continuous fractional resolvent operator. Numerical experiments are presented for several typical time fractional evolutional equations, including time fractional sub-diffusion equations, fractional linear system and semi-linear F-ODEs. All the numerical results exhibit the typical long-time polynomial decay rate, which is fully consistent with our theoretical predictions.

Numerical oscillations of the θ- method for advanced delay differential equations with piecewise continuous arguments

Numerical Simulation for the System of Ordinary Differential Equations

Numerical Analysis of a Fractional Nonlinear Oscillator with Coordinate-Dependent Mass

Oscillation and Asymptotic Behavior of a Special Delay Third Order Nonlinear Neutral Functional Differential Equation

A Schrödinger Equation with Time-Oscillating Critical Nonlinearity

Oscillation Criteria for Delay and Advanced Differential Equations with Nonmonotone Arguments

Functional Differential Equations with Distributed Deviating Arguments: Oscillatory Features of Solutions

The asymptotic behavior and oscillation for a class of third-order nonlinear delay dynamic equations

Numerical solution of highly oscillatory ordinary differential equations

Oscillation Theorems for Second-order Damped Differential Equations with Distributed Deviating Arguments

Stability and Numerical Analysis of a Coupled System of Piecewise Atangana–Baleanu Fractional Differential Equations with Delays

New oscillation theorems for a class of even-order neutral delay differential equations

Applications of the Integral Equation Method to Delay Differential Equations

Analytical solution of fractional oscillation equation with two Caputo fractional derivatives

Oscillation Criteria for Nonlinear Third-Order Delay Dynamic Equations on Time Scales Involving a Super-Linear Neutral Term

Theta-type convolution quadrature OSC method for nonlocal evolution equations arising in heat conduction with memory

Comments on "A perturbation-iterative method for determining limit cycles of strongly non-linear oscillators" - Reply

A Perturbation–incremental Method for Strongly Non-Linear Non-Autonomous Oscillators

Oscillation for Systems of a Class of Second Order Partial Eifferential Equations with Deviating Arguments

Oscillation theorems for fourth-order quasi-linear delay differential equations

Mittag--Leffler stability of numerical solutions to time fractional ODEs