-In this paper, we study the convergence of the finite difference discretization of a second order elliptic equation with variable coefficients subject to general boundary conditions. We prove that the scheme exhibits the phenomenon of supraconvergence on nonuniform grids, i.e., although the truncation error is in general of the first order alone, one has second order convergence. More precisely, for s ∈ (1/2, 2] the optimal order O(h s )-convergence of the finite difference solution and its gradient appears if the exact solution is in the Sobolev -Slobodetskij space H 1+s (Ω). All error estimates are strictly local. Another result of the paper is a close relationship between finite difference scheme and linear finite element methods combined with a special kind of quadrature. As a consequence, the results of the paper can be viewed as the introduction of a fully discrete finite element method for which the gradient is superclose, i.e., the error of the approximate gradient with respect to the linear interpolation of the solution u is of the second order if u ∈ H 3 (Ω). A numerical example is given.2000 Mathematics Subject Classification: 65N06, 65N12, 65N15, 65N30.