S = np.linspace(1e6, 1000e6)