% Linear Algebra: Theory and Applications % Ward Cheney & David Kincaid % JBPub.com (c) 2008 % File: exp8310.m % Example 8.3.10 A = [1 3 5; 7 -8 1; 2 4 -2] x = [1; 1; 1] for k=1:50 y = A*x; r = y(1)/x(1); x = y/norm(y); z = [k,x'] end %checking pause format long c=eig(A) det(A - c(2)*eye(3))