In this paper we present a triangular mesh simplification algorithm which produces accurate approximations of the original models. The simplification is realized using iterative edge contractions. The accuracy is obtained using a symmetric error metric and generating sample points over the simplified mesh. The sample points are generated using iterative 1 : 4 subdivisions of each triangle. The number of subdivisions for each triangle depends on the area of the triangle.