Superconvergence of high order finite difference schemes based on variational formulation for elliptic equations

The classical continuous finite element method with Lagrangian $Q^k$ basis reduces to a finite difference scheme when all the integrals are replaced by the $(k+1)\times (k+1)$ Gauss-Lobatto quadrature. We prove that this finite difference scheme is $(k+2)$-th order accurate in the discrete 2-norm for an elliptic equation with Dirichlet boundary conditions, which is a superconvergence result of function values.