f = np.linspace(0, 0.49999, 1000)