Update calculate_r_t.m