  We study in this paper the function approximation error of linear
interpolation and extrapolation. Several upper bounds are presented along with
the conditions under which they are sharp. All results are under the
assumptions that the function has Lipschitz continuous gradient and is
interpolated on an affinely independent sample set. Errors for quadratic
functions and errors of bivariate linear extrapolation are analyzed in depth.
