A Robust Exponentially Fitted Scheme for Singularly Perturbed Fractional‐Order Delay Parabolic Equation

This study introduces a fitted numerical approach for solving singularly perturbed time‐fractional parabolic differential equations incorporating a delay term. The stability of the method is rigorously examined using the comparison principle and solution bounds, while its convergence is analyzed through the barrier function approach and the Peano kernel. To validate the theoretical framework, two numerical test examples are considered for different values of the perturbation parameter. The proposed scheme demonstrates more accurate, stable, and uniform convergence compared to existing methods in the literature.