Comprobacion
>
evala(evalm(htranspose(Q)&*Q));
Las columnas de Q son ortonormales
>
evalm(Q&*R-B);