g = np.linspace(-6, 4.5, 501)