Use error bound ES=−(b−a)5180n4f(4)(ξ)E_S = -\frac{(b-a)^5}{180n^4} f^{(4)}(\xi)ES=−180n4(b−a)5f(4)(ξ) for some ξ∈(a,b)\xi \in (a,b)ξ∈(a,b).