Paul Dierckx: Algorithm/algorithmus 42 an algorithm for cubic spline fitting with convexity constraints. Computing 24(4): 349-371 (1980)