Computation of Smoothing and Interpolating Natural Splines via Local Bases

It is shown how smoothing splines can be represented in terms of a local basis, and that the coefficients can be obtained by solution of a banded linear system. Recursion relations are developed which permit rapid and accurate calculation of the necessary basis elements.