“…Collocation approximate solutions of nonlinear systems of VIDEs with delay arguments of the general type, N 1 (t) = f 1 (N 1 (t), N 2 (t)) + t t−τ F 1 (t − s)G 1 (N 1 (t), N 2 (s))ds, N 2 (t) = f 2 (N 1 (t), N 2 (t)) + t t−τ F 2 (t − s)G 2 (N 2 (t), N 1 (s))ds, (3) which includes system (1) as a special case were analyzed in [14] based on a constant stepsize strategy. Discrete RK methods and their numerical stability for VIDEs with constant delay have been investigated in [15,16,17] assuming a fixed stepsize strategy. These methods are mainly extensions of the classical Pouzet and Beltyukov RK methods where the discretization is a part of the method and no continuous approximation is produced.…”