Interpolation by a local taut cubic piecewise polynomial
The interpolation of a discrete set of data, on the interval [a,b], representing the function g(x) is obtained using a local piecewise polynomial of order 3. This piecewise cubic interpolant has the following properties: monotonicity and/or convexity or turning points that are present in the data are preserved and no extraneous turning points are created. This approximation method is also presented in the form of FORTRAN 77 subroutine.
scientific article; zbMATH DE number 3986513