function s_set = CS_s_set(n,k,num_s)
%
%   n (int): H_n
%   k (int): singleton point in an H_k
%   num_s (int): number of singleton point H_k's desired
%
%   s_set (num_sxn array): s_set descriptors


s_set = zeros(num_s,n);

s = 0;
while s<num_s
    s = s + 1;
    ss = round(rand(1,n));
    indexes = randperm(n);
    ss(indexes(1:k)) = -1;
    s_set(s,:) = ss;
    if not new
        lower s and loop
    end
end
