Abstract 
Consider a triangulation of the xy plane, and a general surface z = ƒ (xy). The points of the triangle, when lifted to the surface, form a linear spline approximation to the surface. We are interested in the error between the surface and the linear approximant. In fact, we are interested in building triangulations in the plane such that the induced linear approximant is nearoptimal with respect to a given error.
Here we describe a new method, which iteratively adds points to a "Delaunaylike" triangulation of the plane. We locally approximate ƒ by a quadratic surface and utilize this quadrilateral that enables us to minimize the error between the surface and the triangulation.
