An algorithm for computing shape-preserving interpolating splines of arbitrary degree

摘要

The aim of this paper is to describe an algorithm for computing co-monotone and/or co-convex splines of degree m and deficiency m − k at the knots (0 < k < m − k), which are interpolant or osculatory to a given set of data. The method is based upon some existence properties recently developed. Graphical examples and a listing of the FORTRAN code SPISP1 are given.