Computational method for singularly perturbed delay differential equations with twin layers or oscillatory behaviour

In this paper, we have presented a computational method for solving singularly perturbed delay differential equations with twin layers or oscillatory behaviour. In this method, the original second order singularly perturbed delay differential equation is replaced by an asymptotically equivalent first order neutral type delay differential equation. Then, we have employed numerical integration and linear interpolation to get tridiagonal system. This tridiagonal system is solved efficiently by using discrete invariant imbedding algorithm. Several model examples are solved, and computational results are presented by taking various values of the delay parameter and perturbation parameter. We have also discussed the convergence of the method.