A numerical method for singularly perturbed weakly coupled system of two second order ordinary differential equations with discontinuous source term (Q875362)

The authors present a numerical method for a singularly perturbed weakly coupled system of two reaction-diffusion Dirichlet boundary-value problems with a discontinuous source term. Here, the classical finite difference schemes are used to solve the system on a layer-adapted Shishkin mesh. Because of the discontinuity in the source term there is an interior layer in addition to the boundary layers. The Shishkin mesh is constructed following the one-dimensional case. Singular perturbation parameter-uniform error estimates are derived for the proposed scheme. Two numerical examples are carried out to verify the theoretical results. Although the theoretical error estimate shows the first order uniform convergence up to a logarithmic factor, the first component of the solution \(y_1\) has a convergence rate even less than half in both the examples. The second component \(y_2\) reflects the theoretical rate of convergence.