Numerical solutions of system of two-dimensional Volterra integral equations via operational matrices of hybrid functions

The main target of this paper is to solve a system of two-dimensional volterra integral equations (2-dvies). operational matrices of two-dimensional hybrid of block-pulse functions and legendre polynomials are applied to reduce these systems of integral equations to a system of algebraic equations. the main benefit of these basic functions is their efficiency in dealing with non-sufficiently smooth functions. an error bound is provided and some examples are prepared to verify the applicability of the offered numerical technique.