compute an vector best orthogonal to a list of directions
N = perp(v) N = perp(v, 'robust') % ignore outliers