But does your patch to create_set() + + if (nargin == 1) + y = unique(x)'; ^ maybe require to use non-hermitian transpose operator? I think at present your patch could conjugate a set of complex numbers for example. Best, -Muthu