Applications of mesh adaptation techniques could be found in the numerical solution of PDE's or in the optimal triangulation of surfaces for shape representation or graphic display. The scope of this work is to verify through numerical experiments the effectiveness of some algorithms for the control of the L^2 error norm for piecewise linear approximation on 2D unstructured triangular meshes. The analysis could be extended to parametric surfaces and to the 3D case.