In this ‎article‎‎, ‎an ‎ap‎plied matrix method, which is based on Bernouli Polynomials, has been presented to find approximate solutions of ‎high order ‎Volterra ‎integro-differential‎ equations. Through utilizing this approach, the proposed equations reduce to a system of algebric equations with unknown Bernouli coefficients. A number of numerical ‎illustrations‎ have been ‎solved‎ to ‎assert‎ the ‎credibility‎ and ‎practically‎ of ‎this‎ ‎method