linalg.lib
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 print(jordannf(M)); ==> 2,0,0, ==> 0,3,0, ==> 0,1,3