This paper presents a numerical matrix method based on Bernstein polynomials (BPs) for approximate the solution of a system of m-th order nonlinear Volterra integro-differential equations under initial conditions. The approach is based on operational matrices of BPs. Using the collocation points,this approach reduces the systems of Volterra integro-differential equations associated with the given conditions, to a system of nonlinear algebraic equations. By solving such arising non linear system, the Bernstein coefficients can be determined to obtain the finite Bernstein series approach. Numerical examples are tested and the resultes are incorporated to demonstrate the validity and applicability of the approach. Comparisons with a number of conventional methods are made in order to verify the nature of accuracy and the applicability of the proposed approach. Keywords: Systems of nonlinear Volterra integro-differential equations; The Bernstein polyno- mials and series; Collocation points. 2010 AMS Subject Classication: 34A12, 34A34, 45D05, 45G15, 45J05, 65R20. پرونده مقاله