The Use of Multi-Dimensional Cubic Spline Functions for Regression and Smoothing