\[\int_{a}^{b}f(x)\,\mathrm{d}x=\tfrac{1}{3}h(f_{0}+4f_{1}+2f_{2}+4f_{3}+2f_{4}+%
\cdots+4f_{n-1}+f_{n})+E_{n}(f),\]