linalg.lib
LIB "linalg.lib"; ring r=0,(x),dp; matrix A[4][4]=6,0,0,0,0,0,6,0,0,6,0,0,0,0,0,6; print(A); ==> 6,0,0,0, ==> 0,0,6,0, ==> 0,6,0,0, ==> 0,0,0,6 diag_test(A); ==> 1