In this paper, we present a new numerical technique based on Block-pulse functions to solve two-dimensional Volterra-Fredholm integral equations of the second kind. To produce Block-pulse functions, the orthogonal Legendre polynomials is used. Furthermore, operational matrix is applied to convert two-dimensional Volterra-Fredholm integral equations to a linear algebraic system. The convergence analysis of the new method is discussed. Finally, some numerical examples are given to confirm the applicability and efficiency of the new method for solving two-dimensional Volterra-Fredholm integral equations of the second kind. Manuscript profile