Algorithm 547

Algorithm 547 Fortran routines for discrete cubic spline interpolation and smoothing.