LAUSR

ABSTRACT In this study, a numerical method is presented in order to approximately solve singularly perturbed second-order differential equations given with boundary conditions. The method uses the set of monomials whose degrees do not exceed a prescribed as the set of base functions, resulting from the supposition that the approximate solution is a polynomial of degree whose coefficients are to be determined. Then, following Galerkin’s approach, inner product with the base functions are applied to the residual of the approximate solution polynomial. This process, with a suitable incorporation of the boundary conditions, gives rise to an algebraic linear system of size . The approximate polynomial solution is then obtained from the solution of this resulting system. Additionally, a technique, called residual correction, which exploits the linearity of the problem to estimate the error of any computed approximate solution is discussed briefly. The numerical scheme and the residual correction technique are illustrated with two examples.