The development of theory and algorithms relating to interpolation to data by functions which preserve the monotonicity and/or convexity of the data is presented. The functions used for interpolation are polynomials, piecewise polynomials, polynomial splines and exponential splines. The techniques emphasized are the shape-preserving splines of the author, although some discussion of alternate techniques is given.