Multi-order fractional differential equations and their numerical solution. (English) Zbl 1060.65070

The paper presents the numerical solution of multi-order fractional differential equations of the general (possibly nonlinear) form \[ y^{(\alpha )}(t)=f(t, y(t), y^{(\beta_1)}(t), y^{(\beta_2)}(t),\dots, y^{(\beta_n)}(t)) \] with \(\alpha>\beta_n>\beta_{n-1}>...>\beta_1\) and \(\alpha -\beta_n\leq1\), \(\beta_j -\beta_{j-1}\leq 1\) for all \(j\) and \(0<\beta_1\leq 1\). Its linear case is \(y^{(\alpha )}(t)=\lambda_oy(t)+ \sum_{j=1}^n\lambda_jy^{(\beta_j)}(t)+f(t)\). The initial conditions have the form \(y^{k}(t)=y_o^{(k)}\), \(k=0,1,\dots,\lceil\alpha\rceil -1\). The derivatives are understood in the Caputo sense.
A generalization of an approach employed in the solution of ordinary differential equations of order two or higher converting such equation to a system of equations of order one is used. It uses the fact that any real number can be approximated arbitrarily closely by a rational number. Thereby, the assumption on the commensuracy for fractional order equations can be ensured by an appropriate order approximation. A simple generalization of the theorem on the equivalence of a nonlinear system and the linear systems theory is presented first. Then, the nonlinear problem includes two Gronwall-type results for a two-term equation, the general existence-uniqueness as well as the structural stability results. A convergent and stable Adams-type numerical method is proposed including a specific numerical example.


65L05 Numerical methods for initial value problems involving ordinary differential equations
65L20 Stability and convergence of numerical methods for ordinary differential equations
34A34 Nonlinear ordinary differential equations and systems
26A33 Fractional derivatives and integrals
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations


Full Text: DOI


[3] Diethelm, K., An algorithm for the numerical solution of differential equations of fractional order, Elect. Transact. Numer. Anal., 5, 1-6 (1997) · Zbl 0890.65071
[4] Lubich, C., Fractional linear multistep methods for Abel-Volterra integral equations of the second kind, Math. Comp., 45, 463-469 (1985) · Zbl 0584.65090
[5] Brunner, H.; van der Houwen, P. J., The Numerical Solution of Volterra Equations (1986), North Holland: North Holland Amsterdam · Zbl 0611.65092
[6] Diethelm, K.; Ford, N. J., Numerical solution of the Bagley-Torvik equation, BIT, 42, 490-507 (2002) · Zbl 1035.65067
[9] Podlubny, I., Fractional Differential Equations (1999), Academic Press: Academic Press San Diego · Zbl 0918.34010
[10] Diethelm, K.; Freed, A. D., On the solution of nonlinear fractional differential equations used in the modeling of viscoplasticity, (Keil, F.; Mackens, W.; Voß, H.; Werther, J., Scientific Computing in Chemical Engineering II-Computational Fluid Dynamics, Reaction Engineering, and Molecular Properties (1999), Springer: Springer Heidelberg), 217-224
[11] Oldham, K. B.; Spanier, J., The Fractional Calculus, Vol. 111 of Mathematics in Science and Engineering (1974), Academic Press: Academic Press New York, London
[12] Samko, S. G.; Kilbas, A. A.; Marichev, O. I., Fractional Integrals and Derivatives: Theory and Applications (1993), Gordon and Breach: Gordon and Breach Yverdon · Zbl 0818.26003
[13] Gaul, L.; Klein, P.; Kempfle, S., Damping description involving fractional operators, Mech. Syst. Signal Process., 5, 81-88 (1991)
[14] Glöckle, W. G.; Nonnenmacher, T. F., A fractional calculus approach to self-similar protein dynamics, Biophys. J., 68, 1, 46-53 (1995)
[15] Mainardi, F., Fractional calculus: Some basic problems in continuum and statistical mechanics, (Carpinteri, A.; Mainardi, F., Fractals and Fractional Calculus in Continuum Mechanics (1997), Springer: Springer Wien), 291-348 · Zbl 0917.73004
[16] Metzler, R.; Schick, W.; Kilian, H.-G; Nonnenmacher, T. F., Relaxation in filled polymers: A fractional calculus approach, J. Chem. Phys., 103, 7180-7186 (1995)
[17] Caputo, M., Linear models of dissipation whose \(Q\) is almost frequency independent, II, Geophys. J. Royal Astronom. Soc., 13, 529-539 (1967)
[18] Edwards, J. T.; Ford, N. J.; Simpson, A. C., The numerical solution of linear multi-term fractional differential equations: Systems of equations, J. Comput. Appl. Math., 148, 401-418 (2002) · Zbl 1019.65048
[19] Diethelm, K.; Ford, N. J., Analysis of fractional differential equations, J. Math. Anal. Appl., 265, 229-248 (2002) · Zbl 1014.34003
[21] Diethelm, K.; Freed, A. D.; Ford, N. J., A predictor-corrector approach for the numerical solution of fractional differential equations, Nonlinear Dyn., 29, 3-22 (2002) · Zbl 1009.65049
[23] Ford, N. J.; Simpson, A. C., The numerical solution of fractional differential equations: Speed versus accuracy, Numer. Algorithms, 26, 333-346 (2001) · Zbl 0976.65062
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.