AbstractThis paper presents an efficient numerical approach for first-order delay differential equations containing a piece-wise constant delay. The strategy is based on a five-point hybrid block method that has been developed for ordinary differential equations. We will use the interpolation technique for the evaluation of delay terms that are not defined at the grid points. The main characteristics are discussed, including zero stability, local truncation errors, convergence and stability region. The method is easy to implement, and the numerical experiments show the efficiency and accuracy of the proposed method, compared to other methods appeared in the literature.