Estimate error for trapezium rule (#594)

Estimate error for trapezium rule (#594)

Use error bound ET=−(b−a)312n2f′′(ξ)E_T = -\frac{(b-a)^3}{12n^2} f''(\xi) for some ξ∈(a,b)\xi \in (a,b).