x = np.linspace(1.5e-2, 5, 55)