假設鄰接矩陣A,I鏈路度矩陣經n步達矩陣M
則：
B=(A+I)^n=I+A+A^貳+..A^n
再B非零元素改依零元素變即An步達矩陣M；
代碼：
function
M=reach(A,I,n)
%A,In要給定
[row,cow]=size(A)
%處rowcow應該相等
E=ones(row,cow)
%單位矩陣E
B=zeros(row,cow)
%初始全0矩陣B
B=B+I;
%循環+A^n
for
i=依:n
B=B+A^i;
end
%B非零元素改依零元素變,即讓BE進行交運算
M=B&E;
en