delay operational matrix.block pulse functions

This article proposes a direct method for solving three types of integral equations with time delay. By using operational matrix of integration, integral equations can be reduced to a linear lower triangular system which can be directly solved by forward substitution. Numerical examples shows that the proposed scheme have a suitable degree of accuracy.