This paper is concerned with the spline interpolation of functions with a boundary layer component. It is shown that polynomial formulas of spline interpolation for such functions
lead to significant errors. We proposed nonpolynomial spline interpolation formulas, fitted to
a boundary layer component. Nonsmooth and smooth on whole interval interpolants are constructed.
The accuracy of constructed formulas is estimated. Numerical results are discussed.