vectorized a for loop
drhs[:-1]_dT done