In this paper, the numerical technique based on hybrid Bernoulli and Block-Pulse functions has been developed to approximate the solution of system of linear Volterra integral equations. System of Volterra integral equations arose in many physical problems such as elastodynamic, quasi-static visco-elasticity and magneto-electro-elastic dynamic problems. These functions are formed by the hybridi...