Comprobacion
> evalm(htranspose(Q)&*Q): Q es unitaria
> map(simplify,%);
> evalm(inverse(Q)&*B&*Q-T):