Given a table of values of a function and its first derivative, this
subroutine does table look-up and Hermite cubic interpolation to compute an
interpolated value of the function for an arbitrary argument.
Optionally the subroutine can compute values of the first, second, or third
derivative of the interpolation function. All features here are also
available, with more generality, in Chapter 12.1 here
.
This software is part of the MATH77/mathc90 collection here.