On the convergence of spline collocation methods for solving fractional differential equations

In the first part of this paper we study the regularity properties of solutions of initial value problems of linear multi-term fractional differential equations. We then use these results in the convergence analysis of a polynomial spline collocation method for solving such problems numerically. Using an integral equation reformulation and special non-uniform grids, global convergence estimates are derived. From these estimates it follows that the method has a rapid convergence if we use suitable nonuniform grids and the nodes of the composite Gaussian quadrature formulas as collocation points. Theoretical results are verified by some numerical examples.