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)