Given a set of data points @math{(x_1, y_1) \dots (x_n, y_n)} the routines described in this section compute a continuous interpolating function @math{y(x)} such that @math{y_i = y(x_i)}. The interpolation is piecewise smooth, and its behavior at the points is determined by the type of interpolation used.