What about Stineman interpolation ( https://pages.uoregon.edu/dgavin/software/stineman.pdf )?Could someone suggest interpolation algorithm alternative to cubic splines? Interpolated function has to be continuous, with continuous first differential and perfect fit with input data.