cosmo sample unique hdr

function samples = cosmo_sample_unique(k, n, count, varargin)