Before revision:  We adopt this methodology to derive finite-difference schemes for heat equation with Dirichlet and Neumann boundary conditions, whose convergence rates are higher than their truncation errors.
After revision:  We adopt this methodology to design BFD schemes for the heat equation with Dirichlet and Neumann boundary conditions, whose convergence rates are higher than their truncation errors.