wolffd@0: function sep = graph_separated(G, X, Y, S)
wolffd@0: 
wolffd@0: G2 = G;
wolffd@0: G2(S,:) = 0;
wolffd@0: G2(:,S) = 0;
wolffd@0: conn = reachability_graph(G2);
wolffd@0: conn2 = conn(X,Y);
wolffd@0: sep = all(conn2(:)==0);
wolffd@0: