Approximate solutions of delay parabolic equations with the Neumann condition

An approximate solution of the initial-boundary value problem for the delay parabolic partial differential equation is considered. Stable difference schemes of first and second orders of accuracy for this problem are investigated. Convergence estimates for the solution of these difference schemes in Hölder norms are established. Theoretical statements are supported by numerical examples.