Here we present an approximation method for a rather broad class of first ordervariational problems in spaces of piece-wise constant functions overtriangulations of the base domain. The convergence of the method is based on aninequality involving $W^{-1, p}$ norms obtained by Nečas and on the generalframework of Γ-convergence theory.