Integro-differential equations with time-varying delay can provide us with realistic models of many real world phenomena. Delayed Lotka-Volterra predator-prey systems arise in ecology. We investigate the numerical solution of a system of two integro-differential equations with time-varying delay and the given initial function. We will present an approach based on $q$-step methods using quadrature formulas.