Heat equation solver by Thomas algorithm (tridiagonal matrix)