function niantongxi(z,m)
z=[6 10;10 13;13 8;10 1;1 5;5 11;5 9;9 2;3 12;12 7;12 4]   %p12圖2-2一個(gè)具有兩個(gè)連通系的信息矩陣
m=[0 1 1 0 0 1 0 1 0 0 1 0 0]
d=qiu_d(z)
[m1,m2]=size(m)
vs=zeros(m1,m2)
qc=[]
n=0
[m3,mm]=find(m==1)
[m5,m6]=size(mm)
for i=1:m6
    qc(1)=mm(i)
    n=n+1
    [e,f]=find(vs==0)
    [g,h]=size(f)
    if (h==0)                  %當(dāng)vs的每個(gè)元素均不為0時(shí)結(jié)束循環(huán)
        break
    end
    while(length(qc)~=0)
        a=qc(1)
        qc(1)=[]
        vs(a)=n
         for j=1:m2
             if(d(a,j)==1&vs(j)~=n)
                    qc(length(qc)+1)=j
             end
         end
     end
end
aa=max(vs)
bb=length(vs)
for i=1:aa
    gg=1
    for j=1:bb
        if (vs(j)==i)
            s{i}(gg)=j
            gg=gg+1
        end
    end
end
fprintf('此配電網(wǎng)的連通系為:\n')
for i=1:aa
    fprintf('s[%d]',i)
    disp(s{i})
end