计算图的可达矩阵MATLAB程序

% 计算图的可达矩阵
function P=dgraf(A)
n=size(A,1) ;
P=A;
% 计算矩阵Bn
for i=2:n
P=P+A^i;
end
P(P~=0) =1; % 将不为0 的元素变为1
P;
end