t = t_vec[-1] - t_vec[0]