close?all;
clear;
clc

%計算可達矩陣%
A=xlsread（‘adjacency_matrix.xlsx‘）;?%讀入鄰接矩陣
n=size（A）;
E=eye（n）;
R=A+E;
k=0;
while?1
???Rnew=R*（A+E）>0;
???if?isequal（RRnew）;
???????k=k+1;
break;
???end
???R=Rnew;
???k=k+1;
end
disp（Rnew）;
A=Rnew;

%級別劃分
r=1;
M=zeros（n）;
while（~isequal（AM））
%區域劃分%
????for?i=1:n
????????P=find（A（i:））;
????????Q=f