Method of interpolation for a boundary layer problem

A singularly perturbed boundary value problem for a second order ordinary differential equation is considered. It is assumed that the solution is found at the nodes of a uniform or nonuniform mesh. An interpolation method taking into account the boundary layer part of the solution is proposed. Using the constructed interpolation function, we find the derivative of the solution with an accuracy uniform with respect to a parameter at any point of the interval.