> A:=Matrix([[a,b],[c,d]]);

> B:=Matrix([[p,q],[r,s]]);

> with(linalg);

Warning, premature end of input

Warning, the protected names norm and trace have been redefined and unprotected

> trace(A.B);

> trace(B.A);

> trace(A.B)-trace(B.A);

>