Comprobacion

> rank(augment(S,w)); w es un vector de S

> evalm(transpose(v-w)&*S); v-w es ortogonal a S