function ns = neighbors(adj_mat, i)
%
% ns = neighbors(adj_mat, i)

ns = union(children(adj_mat, i), parents(adj_mat, i));

