function f = family(A,i)

f = union(parents(A,i), i);
