from __future__ import print_function
​
import numpy as np
from scipy.integrate import simps
from numpy import trapz
​
​
# The y values. A numpy array is used here,
# but a python list could also be used.
y = np.array([5, 20, 4, 18, 19, 18, 7, 4])
​
# Compute the area using the composite trapezoidal rule.
area = trapz(y, dx=5)
print("area =", area)
​
# Compute the area using the composite Simpson's rule.
area = simps(y, dx=5)
print("area =", area)