We propose a local method for the construction of differentiable funct
ions which interpolate a set of gridded data and are monotonicity pres
erving, that is increasing or decreasing exactly where the data are. T
he scheme is based upon the boolean sum of cubic interpolating operato
rs, which blend together shape-preserving, one dimensional, interpolat
ing polynomials of adaptive degree. Error estimates are also provided.