Searching for C/C++ implementation of Hermite Interpolation. Please suggest simple/readable/academic program.
Inputs: x and y arrays of doubles with 4 values each. Representing 4 2d points. Please note I do not have tangents with me.
Output: Now given any new value of x (within range of given x values), need to find new y. Please note I do not have parameter 't' with me.
Please suggest formulation and show C/C++ implementation.