function ps = children(adj_mat, i)
%
% c = children(adj_mat, i)
% Return the list of children of node i

ps = find(adj_mat(i,:));
