add function to calculate thermal gradient for velocity gradient

add function of Runge-Kutta arithmatic, but not finished.