Estimate error for Simpson’s rule (#595)

Estimate error for Simpson’s rule (#595)

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