Abstract-A new and effective direct method to determine the numerical solution of specific nonlinear Volterra-Fredholm integral and integro-differential equations is proposed. The method is based on vector forms of block-pulse functions (BPFs). By using BPFs and its operational matrix of integration, an integral or integro-differential equation can be transformed to a nonlinear system of algebraic equations.Some numerical examples are provided to illustrate accuracy and computational efficiency of the method. Finally, the error evaluation of this method is presented. The benefits of this method are low cost of setting up the equations without applying any projection method such as Galerkin, collocation, . . . . Also, the nonlinear system of algebraic equations is sparse.