A numerical integration method for solving singular perturbation problems (Q917232)

This paper deals with linear singularly perturbed two-point boundary value problems with a boundary layer on the left end of the interval. The original problem is approximated by a first order differential equation with a small deviating argument. Then classical numerical methods (Simpson's rule, Taylor's expansions and finite differences) are used, together with the given boundary conditions, to obtain a tridiagonal algebraic system of equations. The solution of this system is an approximation of the solution of the original problem. The proposed method is iterative on the deviating argument. Error analysis and the solvability of the tridiagonal system are not considered here.