wolffd@0: % Do the example in fig 23.4 p479 of Cormen, Leiserson and Rivest (1994)
wolffd@0: 
wolffd@0: u = 1; v = 2; w = 3; x = 4; y = 5; z = 6;
wolffd@0: n = 6;
wolffd@0: dag=zeros(n,n);
wolffd@0: dag(u,[v x])=1;
wolffd@0: dag(v,y)=1;
wolffd@0: dag(w,[y z])=1;
wolffd@0: dag(x,v)=1;
wolffd@0: dag(y,x)=1;
wolffd@0: dag(z,z)=1;
wolffd@0: 
wolffd@0: [d, pre, post, cycle, f, pred] = dfs(dag, [], 1);
wolffd@0: assert(isequal(d, [1 2 9 4 3 10]))
wolffd@0: assert(isequal(f, [8 7 12 5 6 11])
wolffd@0: assert(cycle)
wolffd@0: 
wolffd@0: % Now give it an undirected cyclic graph
wolffd@0: G = mk_2D_lattice(2,2,0);
wolffd@0: % 1 - 3
wolffd@0: % |   |
wolffd@0: % 2 - 4
wolffd@0: [d, pre, post, cycle, f, pred] = dfs(G, [], 0);
wolffd@0: % d = [1 2 4 3]
wolffd@0: assert(cycle)
wolffd@0: 
wolffd@0: % Now break the cycle
wolffd@0: G(1,2)=0; G(2,1)=0;
wolffd@0: [d, pre, post, cycle, f, pred] = dfs(G, [], 0);
wolffd@0: assert(~cycle)