This paper deals with spline collocation methods for fractional differential equations, introduced by Pedas and Tamme (2014). Some practical formulas are derived, for the computation of fractional integrals involved in the method, useful for implementation. Linear stability analysis is carried out and stability regions of several methods are provided. Numerical experiments on linear and nonlinear test problems confirm theoretical expectations