The spline collocation method for solving boundary value problems is adapted for problems with singular perturbation. Quadratic spline is used as an approximating function. The exponential features of the exact solution are transfered to spline coefficients by the fitting factor. In this way a uniform stability of the system of matrix is achieved, while the convergence for fixed $\varepsilon$ is preserved.