wolffd@0: function ps = parents(adj_mat, i)
wolffd@0: % PARENTS Return the list of parents of node i
wolffd@0: % ps = parents(adj_mat, i)
wolffd@0: 
wolffd@0: ps = find(adj_mat(:,i))';