Procedure from library linalg.lib (see linalg_lib).
linalg.lib
jordannf(M); matrix M
eigenvalues of M in basefield
matrix J; Jordan normal form of M
Example:
LIB "linalg.lib"; ring R=0,x,dp; matrix M[3][3]=3,2,1,0,2,1,0,0,3; print(M); print(jordannf(M));