The two-point boundary layer problem, described by the second order differential equation with coefficients depending on the small perturbation parameter, is considered. The approximate solution is constructed as a sum of the layer function and a truncated orthogonal series, using the asymptotic expansion of the first order. The theoretical results are illustrated by two numerical examples.