T = np.logspace(1, 5, 10)