In PDE, discretization of 1D heat leads to a tridiagonal system — Thomas algorithm O(n)