Implicit midpoint rule and extrapolation to singularly perturbed boundary alue problems

Linear singularly perturbed boundary value probleme εy″?py = f(x), y(0) = y(l) = 0 is solved numerically by reducing to the first order linear system and applying the implicit midpoint rule on equidistant meshes. Using the asymptotic expansion of the global error, the second order of convergence is improved by Richardson extrapolation when h ²≤ε. Some numerical examples are given in illustration of this theory.