(* Linear Algebra: Theory and Applications *)
(* Ward Cheney & David Kincaid *)
(* JBPub.com (c) 2008 *)
(* File: exp8310.mat *)
(* Example 8-3.10 *)

A = {{1 3 5},{7 -8 1},{2 4 -2}}
x = {1, 1, 1}
For[k=1, k<<=50, k++
  y = A.x
  r = y(1)/x(1)
  x = y/norm(y)
]

c=Eigenavalues[A]
det(A - c(2)*eye(3))