linalg.lib
LIB "linalg.lib"; ring r=0,(t,x),lp; matrix A[2][2]=1,x2,x,x2+3x; print(A); ==> 1,x2, ==> x,x2+3x matrix Adj[2][2]=adjoint(A); print(Adj); //Adj*A=det(A)*E ==> x2+3x,-x2, ==> -x, 1 print(Adj*A); ==> -x3+x2+3x,0, ==> 0, -x3+x2+3x