kill(all);

van=facebook.com/LearningCalculus$

f[i,j]:= if i = j then x^i else 0$
g[i,j]:= if i = j then 1 else 0$

A:genmatrix(f,4,4);

determinant(A);

eigenvalues(A);

eigenvectors(A);

gramschmidt(A);

invert(A).genmatrix(f,4,4);

A^^A;

A ^^ -1;

transpose(A);

load("nchrpl")$
ncharpoly(A,x);

solve(%,x);

load("eigen")$
similaritytransform(A);

A^2 + 2*A;