A numerical study of fractional relaxation–oscillation equations involving ψ-Caputo fractional derivative

We provide a numerical method to solve a certain class of fractional differential equations involving ψ -Caputo fractional derivative. The considered class includes as particular case fractional relaxation–oscillation equations. Our approach is based on operational matrix of fractional integration of a new type of orthogonal polynomials. More precisely, we introduce ψ -shifted Legendre polynomial basis, and we derive an explicit formula for the ψ -fractional integral of ψ -shifted Legendre polynomials. Next, via an orthogonal projection on this polynomial basis, the problem is reduced to an algebraic equation that can be easily solved. The convergence of the method is justified rigorously and confirmed by some numerical experiments.publishe

Uniform approximation of fractional derivatives and integrals with application to fractional differential equations

It is well known that for every f ε Cm there exists a polynomial pn such that pn (k) → f(k), k = 0,..m,Here we prove such a result for fractional (non-integer) derivatives. Moreover, a numerical method is proposed for fractional differential equations. The convergence rate and stability of the proposed method are obtained. Illustrative examples are discussed