ABSTRACT The aim of this report is to show the method of numerical integration for solving second order linear singularly perturbed delay differential equations with negative shift, whose solution has boundary layer of the left and right ends of the solution domain. First, the second order singularly perturbed delay differential equation is replaced by an asymptotically equivalent first order delay differential equation. Then, trapezoidal integration formula along with linear interpolation are employed to get three term recurrence relation which is solved easily by Discrete Invariant Imbedding Algorithm. The efficiency of the method is demonstrated by implementing it on several model examples by taking various values for the delay parameter, 𝛿 and the perturbation parameter, 𝜀.