Apply trapezium rule (#590)

Apply trapezium rule (#590)

Approximate ∫abf(x)dx≈h2[f(x0)+2∑i=1n−1f(xi)+f(xn)]\int_a^b f(x) dx \approx \frac{h}{2} \left[ f(x_0) + 2\sum_{i=1}^{n-1} f(x_i) + f(x_n) \right] with h=b−anh = \frac{b-a}{n}.