In this paper, the system of Fredholm integral equations of the second kind are investigated by using the direct computation method. A degenerate approxima- tion of a non-degenerate kernel is obtained by using the interpolation. Lagrange polynomials are adopted for the interpolation. The error and convergence of the algorithm are given strictly. The efficiency of the approach will be shown by applying the procedure on some ill-posed and well-posed prototype examples.