Higher order global solution and normalized flux for singularly perturbed reaction–diffusion problems

In this paper a computational technique is proposed for obtaining a higher order global solution and global normalized flux of singularly perturbed reaction–diffusion two-point boundary-value problems. The HOC (higher order compact) finite difference scheme developed in Gracia et al. (2001) [4] and which is constructed on an appropriate piecewise uniform Shishkin mesh, has been considered to find an almost fourth order convergent solution at mesh points. Using these values, piecewise cubic interpolants based approximations for solution and normalized flux in whole domain have been defined. It has been shown that the global solution and the global normalized flux are also uniformly convergent. Moreover, for the global solution, the order of uniform convergence in the whole domain is optimal, i.e., it is the same as this one obtained at mesh points, whereas, for the global normalized flux, the uniform convergence is almost third order, except at midpoints of the mesh, where it is also almost fourth order. Theoretical error bounds have been provided along with some numerical examples, which corroborate the efficiency of the proposed technique to find good approximations to the global solution and the global normalized flux.