smooth functions and in their derivatives. The method is constructed from a local stencil of grid point values and is based on a polynomial annihilation technique. By varying the order of the method ...