This article describes a general-purpose method for computing interpol
ating polynomial splines with arbitrary constraints on their shape and
satisfying separable or nonseparable boundary conditions. Examples of
applications of the related Fortran code are periodic shape-preservin
g spline interpolation and the construction of visually pleasing close
d curves.