Abstract: In this paper, we present a numerical integration method to solve singularly perturbed
delay differential equations. In this method, we first convert the second order singularly
perturbed delay differential equation to first order neutral type delay differential equation and
employ the numerical integration. Then, linear interpolation is used to get three term recurrence
relation which is solved easily by discrete invariant imbedding algorithm. The method is
demonstrated by implementing several model examples by taking various values for the delay
parameter and perturbation parameter.