EXPONENTIALLY FITTED NUMERICAL SCHEME FOR SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS INVOLVING SMALL DELAYS

This paper deals with numerical treatment of singularly perturbed differential equations involving small delays. The highest order derivative in the equation is multiplied by a perturbation parameter epsilon taking arbitrary values in the interval (0, 1]. For small epsilon, the problem involves a boundary layer of width O(epsilon), where the solution changes by a finite value, while its derivative grows unboundedly as epsilon tends to zero. The considered problem contains delay on the convection and reaction terms. The terms with the delays are approximated using Taylor series approximations resulting to asymptotically equivalent singularly perturbed BVPs. Inducing exponential fitting factor for the term containing the singular perturbation parameter and using central finite difference for the derivative terms, numerical scheme is developed. The stability and uniform convergence of difference schemes are studied. Using a priori estimates we show the convergence of the scheme in maximum norm. The scheme converges with second order of convergence for the case epsilon = O(N-1) and for the case epsilon << N-1, the scheme converge uniformly with first order of convergence, where N is number of mesh intervals in the domain discretization. We compare the accuracy of the developed scheme with the results in the literature. It is found that the proposed scheme gives accurate result than the one in the literatures.