dt1 = t_vec[1] - t_vec[0]