A Chebyshev spectral method for solving Riemann–Liouville fractional boundary value problems

The authors derive a series of explicit formulas to approximate the Riemann–Liouville derivative and integral of arbitrary order by shifted Chebyshev polynomials. This is then applied to solve boundary value problems involving Riemann–Liouville derivatives.