A numerical method for a system of singularly perturbed reaction–diffusion equations

A Dirichlet problem for a system of two coupled singularly perturbed reaction–diffusion ordinary differential equations is examined. A numerical method whose solutions converge pointwise at all points of the domain independently of the singular perturbation parameters is constructed and analysed. Numerical results are presented, which illustrate the theoretical results.