This paper presents a uniformly accurate difference approximation for a system of singularly perturbed reaction-diffusion equations with delay. The proposed method utilizes an appropriate combination of exponential and cubic spline difference schemes. It employs grid equidistribution to address the challenges posed by the multiscale nature of these systems, which often feature sharp gradients and boundary layers. The grid is generated based on the equidistribution of a positive monitor function,
