function [the_genes] = show_clusters(all_genes,U,show_me) [nclust,N] = size(U); j = 1; for i = 1:N, if U(show_me,i) == 1, the_genes(j,:) = all_genes(i,:); j = j + 1; end end