Tp = np.linspace(2, 16)