indAlpha = np.argwhere(np.diff(np.sign(f - z))).flatten()