Superconvergence of Centered Finite Difference Approximations

Abstract

The finite difference (FD) method is commonly used to approximate derivatives of smooth functions, with accuracy typically determined by stencil size and derivative order. However, certain centered stencils exhibit unexpectedly higher accuracy, a phenomenon known as superconvergence, which has been observed in practice but lacks rigorous explanation. We present a mathematical framework for superconvergence in centered FD approximations based on Taylor expansions of the truncation error and the resulting linear system for the FD coefficients. By analyzing symmetry properties of these coefficients and their interaction with the parity of the derivative order, we identify conditions under which higher-order error terms cancel. We show that superconvergence occurs for odd-order derivatives with even centered stencils and for even-order derivatives with odd centered stencils, while no superconvergence occurs for even derivatives with even centered stencils. Numerical experiments in MATLAB confirm the predicted convergence rates.

To be published in PUMP Math Journal