Comprobando

> op(Q);

> H:=evalm(inverse(Q)&*B&*Q);