  Functions in a Sobolev space are approximated directly by piecewise affine
interpolation in the norm of the space. The proof is based on estimates for
interpolations and does not rely on the density of smooth functions.
