A uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction-diffusion problems

A coupled system of two singularly perturbed linear reaction–diffusiontwo-point boundary value problems is examined. The leading termof each equation is multiplied by a small positive parameter,but these parameters may have different magnitudes. The solutionsto the system have boundary layers that overlap and interact.The structure of these layers is analysed, and this leads tothe construction of a piecewise-uniform mesh that is a variantof the usual Shishkin mesh. On this mesh central differencingis proved to be almost first-order accurate, uniformly in bothsmall parameters. Supporting numerical results are presentedfor a test problem.