One dimensional initial boundary delay pseudo-parabolic problem is being considered. To solve this problem numerically, we construct higher order difference method for approximation to the considered problem and obtain the error estimate for its solution. Based on the method of energy estimate the fully discrete scheme is shown to be convergent of order four in space and of order two in time. Numerical example is presented. (C) 2017 Elsevier B.V. All rights reserved.