import numpy as np from scipy.differentiate import derivative def f(x): return np.sin(x) x = np.pi / 4