Hm0 = np.linspace(1, 20)