High Order Approximations of Solutions to Initial Value Problems for Linear Fractional Integro-Differential Equations

Abstract

We consider a general class of linear integro-differential equations with Caputo fractional derivatives and weakly singular kernels. First, the underlying initial value problem is reformulated as an integral equation and the possible singular behavior of its exact solution is determined. After that, using a suitable smoothing transformation and spline collocation techniques, the numerical solution of the problem is discussed. Optimal convergence estimates are derived and a superconvergence result of the proposed method is established. The obtained theoretical results are supported by numerical experiments.