Abstract
In this paper, an implementation of an efficient numerical method of linear fractional integro-differential equations (LFIDEs) by least squares method with aid of shifted Chebyshev polynomials of the third kind method. The fractional derivative is described in the Caputo sense. The method is based upon shifted Chebyshev polynomials of the third kind approximations is introduced. Some numerical examples are presented to illustrate the theoretical results and compared with the results obtained by other numerical methods. We have computed the numerical results using Mathematica 9 programming.