Abstract
In this paper, we implement Chebyshev pseudo-spectral method for solving numerically system of linear and non-linear fractional integro-differential equations of Volterra type. The proposed technique is based on the new derived formula of the Caputo fractional derivative. The suggested method reduces this type of systems to the solution of system of linear or non-linear algebraic equations. We give the convergence analysis and derive an upper bound of the error for the derived formula. To demonstrate the validity and applicability of the suggested method, some test examples are given. Also, we present a comparison with the previous work using the homotopy perturbation method. (C) 2013 Elsevier Inc. All rights reserved.