Abstract
This work presents a projection method based on Vieta-Lucas polynomials and an effective approach to solve a Cauchy-type fractional integro-differential equation system. The suggested established model overcomes two linear equation systems. We prove the existence of the problem's approximate solution and conduct an error analysis in a weighted space. The theoretical results are numerically supported.