Powers of Matrices and Eigenvectors

> with(linalg):

Warning, new definition for norm

Warning, new definition for trace

> A:= matrix(2,2,[[-1,3],[-2,4]]);

> id := diag(1,1);

> vects := 'vects';

> evalf(Eigenvals(A,vects));

> print(vects);

> n := 10;

> v := vector([1,0]);

The power method for finding eigenvalues and eigenvectors.

> for i from 1 to 10 do
v := evalm(A &* v);
v_len := evalf(sqrt(dotprod(v,v))):
v := map(evalf,evalm(1/v_len * v)):
print(v):

od:

>