Comprobacion
>
evalm(transpose(Q)&*Q);
Las columnas de Q son ortonormales
>
evalm(Q&*R-A);