for k in range(n - 1):  # k = subdiagonal