linalg.lib
matrix J; Jordan matrix with list(e,s,m)==jordan(J)
LIB "linalg.lib"; ring R=0,x,dp; ideal e=ideal(2,3); intvec s=1,2; intvec m=1,1; print(jordanmatrix(list(e,s,m))); ==> 2,0,0, ==> 0,3,0, ==> 0,1,3