function [P,v]=dettra(A,B)
P=A*B;
v=[det(P),trace(P)];