In this project report, numerical integration method with exponential integrating factor is presented to solve singularly perturbed delay differential equations with negative shift, whose solution has boundary layer of the left and the right ends of the solution domain. First, the given second order singularly perturbed delay differential equation is replaced by an asymptotically equivalent first order neutral type delay differential equation. An exponential integrating factor is introduced into the first order delay differential equation. Then, trapezoidal rule along with linear interpolation has been employed to get a three term recurrence relation which is solved by Thomas algorithm. Convergence of the proposed method has also discussed. The efficiency of this method is demonstrated by implementing it on two modal examples by taking different values for delay parameters ii  and the perturbation parameters  .Different value of perturbation parameters (  ) and mesh sizes ( h ) are considered and shown in table 1, table 2 and table 3.The effect of delay parameters on the boundary layer solutions has been investigated and presented in (figures 1-10). Maximum absolute errors are computed, tabulated and compared with the result of Kadalbajoo and Sharma (2004).