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(hessenberg(M)); ==> 3,2,1, ==> 0,2,1, ==> 0,0,3