In this study, linear Fredholm fractional integro-differential equations (FIDEs) are numerically solved, where the fractional derivative is considered in the Caputo sense. In this work, the least squares method (LSM) using a compact combination of shifted Chebyshev polynomials (SCP) of the first Kind is applied to solving a class of FIDEs. Our aim is to write the unknown function as a series of a linear combination of SCP, and then to reduce the problem to a system of linear algebraic equations, which will be solved for the unknown constants associated with the approximate solution, using MATLAB R2020a. Finally, numerical examples are presented to confirm the reliability, applicability, and efficiency of this method, in addition, various comparisons are also shown.