Note that the GAUSS Run-Time Library now contains 2- and 3-dimensional cubic spline functions called spline2d() and spline3d().

Cubic Spline Smoothing Routines