Comprobacion
> evalm(transpose(U)&*U); U es ortogonal
> evalm(transpose(V)&*V); V es ortogonal
> equal(A,U&*Sigma&*V);