Computational schemes for the collocation method and the mechanical quadrature method for the approximate solution of systems of singular integro-differential equations with a Cauchy kernel are elaborated. The case where the systems of equations are defined on an arbitrary smooth closed contour of a complex plane is examined. The methods researched are based on Fej茅r points. Estimates of the rate of convergence in Lebesgue space are obtained.