Singularly perturbed integro-partial differential equations with reaction–diffusion behavior present significant challenges due to boundary layers arising from small perturbation parameters, which complicate the development of accurate and efficient numerical schemes for physical and engineering models. In this study, a uniformly convergent higher-order method is proposed to address these challenges. The approach applies the implicit Euler method for temporal discretization on a uniform mesh and central differences on a Shishkin mesh for spatial approximation, and utilizes the trapezoidal rule for evaluating integral terms; further, extrapolation techniques are incorporated in both time and space to increase accuracy. Numerical analysis demonstrates that the base scheme achieves first-order convergence, while extrapolation enhances convergence rates to second-order in time and fourth-order in space. Theoretical results confirm uniform convergence with respect to the perturbation parameter, and comprehensive numerical experiments validate these analytical claims. Findings indicate that the proposed scheme is reliable, efficient, and particularly effective in attaining fourth-order spatial accuracy when solving singularly perturbed integro-partial differential equations of reaction–diffusion type, thus providing a robust numerical tool for complex applications in science and engineering.
Govindarao et al. (Sun,) studied this question.