Multi-order fractional differential equations and their numerical solution

摘要

We consider the numerical solution of (possibly nonlinear) fractional differential equations of the form y(α)(t)=f(t,y(t),y(β1)(t),y(β2)(t),…,y(βn)(t)) with α>βn>βn−1>⋯>β1 and α−βn⩽1, βj−βj−1⩽1, 0<β1⩽1, combined with suitable initial conditions. The derivatives are understood in the Caputo sense. We begin by discussing the analytical questions of existence and uniqueness of solutions, and we investigate how the solutions depend on the given data. Moreover we propose convergent and stable numerical methods for such initial value problems.