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