  In this paper we consider the numerical solution of fractional differential
equations. In particular, we study a step-by-step graded mesh procedure based
on an expansion of the vector field using orthonormal Jacobi polynomials. Under
mild hypotheses, the proposed procedure is capable of getting spectral
accuracy. A few numerical examples are reported to confirm the theoretical
findings.

Numerical Analysis;Numerical Analysis