Comprobacion

> rank(augment(T,z)); z es un vector de T

> evalm(htranspose(u-z)&*T); u-z es ortogonal a T