Reinsch's method for solving the smoothing problem in several dimensions

Reinsch's method for constructing cubic smoothing splines is generalized to the multidimensional case. Solution of the smoothing problem reduces to repeated solution of systems of equations with positive definite matrices. The algorithm is described in detail.