function cluster = induced_cluster( node, graph )
% add each of node and its neighbours to cluster
nset = find(graph(:,node));
cluster = sort([node,nset']);