DIRERCKX (Bell Labs Fortran Routines for Spline Fit or Smoothing)
Netlib (Fortran Routines)
Numerical Recipes