Euler–MacLaurin formulas via differential operators

Recently there has been a renewed interest in asymptotic Euler–MacLaurin formulas, because of their applications to spectral theory of differential operators. Using elementary means, we recover such formulas for compactly supported smooth functions on intervals, polygons, and three-dimensional polytopes, where the coefficients in the asymptotic expansion are sums of differential operators involving only derivatives of the function in directions normal to the faces of the polytope. Our formulas apply to wedges of any dimension.