Numerical treatment of fractional differential models

Several of fracture and fatigue problems are modelled by fractional differential equations. This paper deals with the numerical solution by a class of one and two step spline collocation methods. These methods have higher order of convergence with respect to most numerical methods for fractional differential equations. The methods are illustrated, and their convergence properties are discussed. Several numerical experiments confirm theoretical results and compare one and two step spline collocation methods.