  In this paper, we use a numerical method that involves hybrid and block-pulse
functions to approximate solutions of systems of a class of Fredholm and
Volterra integro-differential equations. The key point is to derive a new
approximation for the derivatives of the solutions and then reduce the
integro-differential equation to a system of algebraic equations that can be
solved using classical methods. Some numerical examples are dedicated for
showing efficiency and validity of the method that we introduce.
