Abstract
A class of non-local boundary value problems for linear fractional differential equations with Caputo-type differential operators is considered. By using integral equation reformulation of the boundary value problem, we study the existence and smoothness of the exact solution. Using the obtained regularity properties and spline collocation techniques, we construct two numerical methods (Method 1 and Method 2) for finding approximate solutions. By choosing suitable graded grids, we derive optimal global convergence estimates and obtain some super-convergence results for Method 2 by requiring additional assumptions on equation and collocation parameters. Some numerical illustrations for verification of theoretical results is also presented.