Double-shifted Chebyshev polynomials are developed in this study to approximate the solutions of the convolution integral, Volterra integral equation, and Fredholm integral equation. This method simplifies the computations of integral equations to the successive solutions of a linear algebraic equation in matrix form. In addition, the computational complexity can be reduced remarkably. Three examples are illustrated. It is seen that the proposed approach is straightforward and convenient, and converges faster in finding approximations than other existing orthogonal function methods.