In this article, a numerical method based on improvement of block-pulse functions (IBPFs) is discussed for solving the system of linear Volterra and Fredholm integral equations. By using IBPFs and their operational matrix of integration, such systems can be reduced to a linear system of algebraic equations. An efficient error estimation and associated theorems for the proposed method are also presented. Some examples are given to clarify the efficiency and accuracy of the method.
