Abstract
The aim of this paper is to develop numerical methods for solving systems of integral equations with cardinal splines. The unknown functions are expressed as a linear combination of horizontal translations of certain cardinal spline functions with small compact supports. Then a simple system of equations on the coefficients is obtained for the system of integral equations. It is relatively straight forward to solve the system of unknowns and an approximation of the original solution with high accuracy is achieved. Several cardinal splines are used in the paper to enhance the accuracy. The sufficient condition for the existence of the inverse matrix is examined and the convergence rate is investigated. Examples are given to demonstrate the benefits of the methods.