Is there an algorithm or implementation? I have searched the monotone cubic interpolation or monotone piecewise cubic interpolation. It seems that both the two methods cannot preserve the derivatives at the two given points. I am not really sure about this, because I am not a mathematician. Can you provide some information or ideas? Anything will be greatly appreciated.