Numerical methods for solving a hereditary equation of hyperbolic type

A family of grid methods is constructed for the numerical solution of a wave equation with delay of general form; the methods are based on the idea of separating the current state and the history function. A theorem on the order of convergence of the methods is obtained by means of embedding into a general difference scheme with aftereffect. Results of calculating test examples with constant and variable delays are presented.