newman快速算法源代码

function [Z] = newman（A）
n=size（A1）;
Z=[];
e=A/sum（sum（A==1））;
clusters=[1:1:n;1:1:n;zeros（1n）]; 
k=1;
while（sum（clusters（3:）==0）~=1）
    a=sum（e）;
    DeltaQ=[];
    for l=1:size（e1）
        for s=l+1:size（e1）
            deltaQ=2*（e（ls）-a（l）*a（s））;
            DeltaQ=[DeltaQ [l;s;deltaQ]];
        end
    end
    [maxDeltaQid]=max（DeltaQ（3:））;
    c1=DeltaQ（1id）;
    c2=DeltaQ（2id）;     
    f