D.3.2.22 jordan

Usage:

jordan(M); matrix M

Assume:

eigenvalues of M in basefield

Return:

list l; Jordan data of M ideal l[1]; number l[1][i]; eigenvalue of i-th Jordan block of M intvec l[2]; int l[2][i]; size of i-th Jordan block of M intvec l[3]; int l[3][i]; multiplicity of i-th Jordan block of M

LIB "linalg.lib";
ring R=0,x,dp;
matrix M[3][3]=3,2,1,0,2,1,0,0,3;
print(M);
==> 3,2,1,
==> 0,2,1,
==> 0,0,3 
jordan(M);
==> [1]:
==>    _[1]=2
==>    _[2]=3
==> [2]:
==>    1,2
==> [3]:
==>    1,1