“…If we take h ≥ ε, the existing classical numerical methods produce oscillatory solution and pollute the solution in the entire interval, because of boundary layer behavior. In connection to this, there are some numerical methods suggested by various authors for solving self-adjoint singular perturbation problems, namely, initial value technique [5], quintic spline method [6], nonpolynomial spline functions method [7], difference scheme using cubic spline [8], finite difference method with variable mesh [9], fitted mesh B-spline collocation method [10], higher order numerical methods [11].…”