This article introduces a new projection method via shifted Legendre polynomials and an efficient procedure for solving a system of integro-differential equations of the Cauchy type. The proposed computational process solves two systems of linear equations. We demonstrate the existence of the solution to the approximate problem and conduct an error analysis. Numerical tests provide theoretical results.