Numerical methods for evolutionary equations with delay and software package PDDE The paper gives a survey of the author’s results on the grid-based numerical algorithms for solving the evolutionary equations (parabolic and hyperbolic) with the effect of heredity on a time variable. From uniform positions we construct analogs of schemes with weights for the one-dimensional heat conduction equation with delay of general form, analog of a method of variable directions for the equation of parabolic type with time delay and two spatial variables, analog of the scheme with weights for the equation of hyperbolic type with delay. For the one-dimensional heat conduction equation and the wave equation we obtained conditions on the weight coefficients that ensure stability on the prehistory of the initial function. Numerical algorithms are implemented in the form of software package partial delay differential equations (PDDE) toolbox.