In this paper, we considered singularly perturbed self–adjoint boundary-value problems and proposed a computational technique based on spline scheme, which is also suitable for parallel computing. The whole domain is divided into three non-overlapping subdomains and corresponding subproblems are obtained by using zeroth order approximations of the solution at boundaries of these subproblems. The subproblems corresponding to boundary layer regions are solved using adaptive spline scheme. Numerical example is provided to show the efficiency and accuracy.
Subject Classification: AMS 65L10 CR G1.7.