source('utils.R')

n = 21 ## official problem size
A = build_tridiagonal_matrix(n, 0, 1, 1)

pm = power_method(A, tol=1e-6, init_choice=2)

exact = exact_tridiagonal_eigensystem(n, 0, 1)
print('exact eigenvalues=')
print(exact$lambda)