This paper presents a theoretical and numerical study of a class of discontinuous Galerkin methods that shows the approximation of the gradient superconverges at the zeros of the Legendre polynomials on a model I D elliptic problem. Numerical experiments validate the theoretical results. (C) 2003 Elsevier B.V. All rights reserved.