pub fn power_iteration(a: &Matrix, max_iter: usize, tol: f64) -> (f64, Vec<f64>)
Power iteration for the dominant eigenvalue/eigenvector.