from scipy.integrate import trapz, simps, quad import numpy as np def f(t): return np.sin(t) x = np.linspace(0, np.pi, num=100, endpoint=True) y = f(x) I3 = quad(f, x[0], x[-1]) print("QUADPACK: ", I3)
