Implements a cubic monotonic Hermite interpolator as described in
Dougherty, R. L., Edelman, A. S., & Hyman, J. M. (1989). Nonnegativity-, monotonicity-, or convexity-preserving cubic and quintic Hermite interpolation. Mathematics of Computation, 52(186), 471-494.