用matlab怎么做DFS呢？

有一个有n个节点的有向图的邻接矩阵array，现在要用深度优先算法得到一片深度优先森林，并将每棵树用不同的编号标记起来，熬了一晚上都没做出来……囧……这个要怎么做呢？麻烦哪位大神帮忙说说算法……