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

with(LinearAlgebra):
A := Matrix([1 3 5],[7 -8 1],[2 4 -2]);
x := Vector(1, 1, 1);
for k from 1 to 50 do
  y := A.x;
  r := y(1)/x(1);
  x := y/norm(y);
  z := [k,x']
do 

%checking
c:=eig(A)
det(A - c(2)*eye(3))