Higher order schemes and Richardson extrapolation for singular perturbation problems

Semilinear singular perturbation problems are solved numerically by using finite–difference schemes on non-equidistant meshes which are dense in the layers. The fourth order uniform accuracy of the Hermitian approximation is improved by the Richardson extrapolation.