Fortran subroutines for computing smoothing and interpolating natural splines

Abstract FORTRAN subroutines for the stable and fast computation of smoothing and interpolating natural splines solving certain data fitting problems in one dimension are presented. The subroutines are based on algorithms discussed in Ref. 1.