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...