Objectives: The main objective of this work is to design an efficient numerical scheme is proposed for solving singularly perturbed time delayed parabolic problems with two parameters.
Results: The scheme is constructed via the implicit Euler and non-standard finite difference method to approximate the time and space derivatives, respectively. Besides, to enhance the accuracy and order of convergence of the method Richardson extrapolation technique is employed.