First, we investigate the existence and uniqueness of solution for a k-dimensional system of multi-term fractional integro-differential equations. Also, we apply shifted Chebyshev and shifted Legendre polynomials to obtain an approximation solution for the k-dimensional system. Finally, we provide some examples illustrating the presented methods.