Comprobacion

> evalm(transpose(Q)&*Q): Q es ortogonal

> map(simplify,%);

> evalm(inverse(Q)&*A&*Q-T):

> map(simplify,%);