Editing Euler–Maclaurin formula (section)

===The remainder term===
{{see also|Bernoulli polynomials}}

The remainder term arises because the integral is usually not exactly equal to the sum.  The formula may be derived by applying repeated [[integration by parts]] to successive intervals {{math|[''r'', ''r'' + 1]}} for {{math|''r'' {{=}} ''m'', ''m'' + 1, …, ''n'' − 1}}.  The boundary terms in these integrations lead to the main terms of the formula, and the leftover integrals form the remainder term.

The remainder term has an exact expression in terms of the periodized Bernoulli functions {{math|''P<sub>k</sub>''(''x'')}}.  The Bernoulli polynomials may be defined recursively by {{math|''B''<sub>0</sub>(''x'') {{=}} 1}} and, for {{math|''k'' ≥ 1}},
<math display=block>\begin{align}
  B_k'(x) &= kB_{k - 1}(x), \\
  \int_0^1 B_k(x)\,dx &= 0.
\end{align}</math>
The periodized Bernoulli functions are defined as
<math display=block>P_k(x) = B_k\bigl(x - \lfloor x\rfloor\bigr),</math>
where {{math|⌊''x''⌋}} denotes the largest integer less than or equal to {{mvar|x}}, so that {{math|''x'' − ⌊''x''⌋}} always lies in the interval {{math|[0,1)}}.

With this notation, the remainder term {{mvar|R<sub>p</sub>}} equals
<math display="block">R_{p} = (-1)^{p+1}\int_m^n f^{(p)}(x) \frac{P_p(x)}{p!}\,dx. </math>

When {{math|''k'' > 0}}, it can be shown that for {{math|0 ≤ ''x'' ≤ 1}},
<math display=block>\bigl|B_k(x)\bigr| \le \frac{2 \cdot k!}{(2\pi)^k}\zeta(k),</math>
where {{mvar|ζ}}  denotes the [[Riemann zeta function]]; one approach to prove this inequality is to obtain the Fourier series for the polynomials {{math|''B<sub>k</sub>''(''x'')}}. The bound is achieved for even {{mvar|k}} when {{mvar|x}} is zero.  The term {{math|''ζ''(''k'')}} may be omitted for odd {{mvar|k}} but the proof in this case is more complex (see Lehmer).<ref name="Lehmer">{{cite journal|last1=Lehmer|first1=D. H.|author-link=D. H. Lehmer |title=On the maxima and minima of Bernoulli polynomials | date=1940 | journal=[[The American Mathematical Monthly]]|volume=47|issue=8|pages=533–538 |doi=10.2307/2303833|jstor=2303833}}</ref> Using this inequality, the size of the remainder term can be estimated as
<math display=block>\left|R_p\right| \leq \frac{2 \zeta(p)}{(2\pi)^p}\int_m^n \left|f^{(p)}(x)\right|\,dx.</math>