A cubic spline function with NB - 1 segments is a function consisting of
NB - 1 pieces, each of which is a cubic polynomial. At the abscissae,
called knots, at which adjacent segments meet, the function has C2
continuity, i.e. continuity in value, first derivative, and second
derivative.
This software is part of the MATH77/mathc90 collection here.