Numerical solution of Volterra integral-algebraic equations using block pulse functions

This paper presents a method for computing numerical solutions for linear Volterra integral-algebraic equations using block pulse functions. The problem is transformed to a linear lower triangular system of algebraic equations using the operational matrix associated with block pulse functions. Convergence result and numerical examples are presented to illustrate the efficiency and applicability of the method.