matlab运行时提示Error: Missing variable or function.
请各位高手看看以下程序的问题出在哪里？谢谢！

function [gtt,ctt]=ccut(g,i,j,ctt,cs)
gtt=g;
l1=find(gtt(i,:)==j);
for k=1:length(l1)
ctt(i,l1(k))=cs(i,l1(k));
ctt(l1(k),i)=ctt(i,l1(k));
gtt(i,l1(k))=10000;
gtt(l1(k),i)=gtt(i,l1(k));
end

l1=find(gtt(j,:)==i);
for k=1:length(l1)
ctt(j,l1(k))=cs(j,l1(k));
ctt(l1(k),j)=ctt(j,l1(k));
gtt(j,l1(k))=10000;
gtt(l1(k),j)=gtt(j,l1(k));
end

function [gtt,ctt]=cupdate(g,i,j,ctt,cs)
gtt=g;
l1=find(gtt(i,:)==10000);
for k=1:length(l1)
if l1(k)~=j && cs(i,j)*cs(j,l1(k))>cs(i,l1(k))
ctt(i,l1(k))=cs(i,j)*cs(j,l1(k));
ctt(l1(k),i)=ctt(i,l1(k));
gtt(i,l1(k))=j;
gtt(l1(k),i)=gtt(i,l1(k));
end
end

l1=find(gtt(j,:)==10000);
for k=1:length(l1)
if l1(k)~=i && cs(j,i)*cs(i,l1(k))>cs(j,l1(k))
ctt(j,l1(k))=cs(j,i)*cs(i,l1(k));
ctt(l1(k),j)=ctt(j,l1(k));
gtt(j,l1(k))=i;
gtt(l1(k),j)=gtt(j,l1(k));
end
end

function [gt,umean,uvar]=groupdemo(cs)
N=size(cs,1);
%----------------------------------------------
%更新量，初始化
ct=zeros(N,N);
ct=cs;
gt=zeros(N,N);%0 代表不连接，10000 代表连接,其他数即为间接连接的中间结
点号

%----------------------------------------------
count=0;
cut=0;
while 1
    % % 所有 ij 不属于 g 都不满足条件 2 时，收敛，退出循环；否则，将 ij加入网络 gt
    % % end
    cflag=1;%判断收敛否
    count=count+1;
    [a,b]=find(gt==0);

    for i=1:length(a)
        gtt=gt;
        ctt=ct;
    if a(i)~=b(i)
           gtt(a(i),b(i))=10000;
           gtt(b(i),a(i))=10000;
           [gtt,ctt]=cupdate(gtt,a(i),b(i),ctt,cs);%更新
           ui_gt=U(gt,a(i),ct);
           ui_gtt=U(gtt,a(i),ctt);
           uj_gt=U(gt,b(i),ct);
           uj_gtt=U(gtt,b(i),ctt);
    if ui_gtt>=ui_gt && uj_gtt>=uj_gt
    gt=gtt;cflag=0;ct=ctt;%更新
    end
end
end
if cflag==1
break;
end
if count==20000
break;
end

% %所有 ij 属于 g,若 ij 属于 g 满足条件 1，连接保持；否则，连接删除
[a,b]=find(gt>0);

    for i=1:length(a)
     gtt=gt;
     ctt=ct;
    if a(i)~=b(i)
     ui_gt=U(gt,a(i),ct);
     uj_gt=U(gt,b(i),ct);
     gtt(a(i),b(i))=0;
     gtt(b(i),a(i))=0;
     [gtt,ctt]=ccut(gtt,a(i),b(i),ctt,cs);%更新
     ui_gtt=U(gtt,a(i),ctt);
     uj_gtt=U(gtt,b(i),ctt);
     if ui_gt>=ui_gtt && uj_gt>=uj_gtt
     gt=gt;
     else
     cut=cut+1;gt=gtt;ct=ctt;%更新
     end
    end
    end
end
    uvar=0;
    umean=0;
    for i=1:N
        u(i)=U(gt,i,ct);
    end
    umean=sum(u)/N;
   uvar=var(u);

   function [u_gt]=U(g,i,c)
   sita=0.8;
   transc=0.8;

   lsum=0;
   costd=0;
   costid=0;
   %与 i 直接连接的 j
   l1=find(g(i,:)==10000);
   for j=1:length(l1)
   lsum=lsum+transc*c(i,l1(j));
   costd=costd+0.01^c(i,l1(j));
   end
   %与 i 间接连接的 j
   l1=find(g(i,:)<10000&g(i,:)>0);
   for j=1:length(l1)
   lsum=lsum+(1-transc)*c(i,l1(j));

   costid=costid+0.01^c(i,l1(j));
   end

   u_gt=1+sita*lsum-1-1.5*costd-1*costid;

   function draw(g,rx,ry)
   N=size(g,1);
   figure;
   for i=1:N
       x=rx(i);
       y=ry(i);
       plot(x,y,'o');
       text(x,y,[' ' num2str(i)]);
       hold on;
   end

   for i=1:N
       count=0;
       for j=1:i
       if g(i,j)==10000
       x=[rx(i),rx(j)];
       y=[ry(i),ry(j)];
       plot(x,y,'-');
       hold on;
   end
   end
end

function test
close all;

N=size(cs,1);
N1=10;
N2=30;

function test
close all;
cs=[1 0.8 0.6 0.7 0.5 0.9 0.8 0.4 0.7 0.3 0.2 0.6 0.2 0.3 0.2 0 0 0
0 0.1 0 0.5 0.2 0.3 0 0 0.2 0 0 0 0.4 0.7 0.1 0 0.1 0.9
0.3 0.3 0.2 0.2 0.3 0.1 0.7 0.3 0.6 0.7 0.6 0.8 0.4 0.5 0.3 0.6 0.7 0.7
0.6 0.4 0.6 0.3 0.9 0.7
0.8 1 0.4 0.6 0.3 0.7 0.8 0.2 0.8 0.1 0.2 0.3 0 0.2 0 0.4 0.5 0.3
0.2 0.3 0 0.8 0.4 0.3 0 0 0 0 0 0.2 0.8 0.8 0.4 0.9 0.2 0.5
0.5 0.4 0.4 0.7 0.6 0.6 0.6 0.2 0.3 0.8 0 0 0.8 0.2 0 0.4 0.8 0.3
0.2 0.3 0.7 0.9 0.4 0.7
0.6 0.4 1 0.4 0.5 0.6 0.8 0.3 0.7 0.2 0.5 0.4 0.2 0.2 0 0 0 0
0.2 0.2 0 0.8 0 0 0 0 0 0 0 0 0.8 0.2 0.5 0.5 0.3 0.9
0.9 0.3 0.2 0.2 0.3 0.7 0.2 0.5 0.2 0.4 0.5 0.2 0.6 0.6 0.6 0.5 0.7 0.3
0.2 0.8 0.2 0.8 0.5 0.9
0.7 0.6 0.4 1 0.2 0.5 0.8 0.2 0.6 0.2 0.4 0.5 0.2 0.2 0.3 0.1 0 0 0
0.4 0.2 0.8 0.2 0.2 0 0 0.3 0.4 0.2 0.3 0 0.7 0.8 0.3 0.1 0.1 0.7 0
0.4 0.4 0.2 0.5 0.5 0.7 0.2 0.9 0.1 0.1 0.6 0.7 0.6 0.1 0.4 0.4 0.7 0.1
0.7 0.8 0.4 0.7
0.5 0.3 0.5 0.2 1 0.2 0.6 0 0.7 0 0.5 0.6 0 0.2 0.3 0.6 0.2 0 0
0 0.1 0.8 0.4 0.3 0.4 0.2 0.3 0 0 0 0.2 0.1 0.5 0.8 0.7 0.9 0.6
0.9 0.6 0.5 0.2 0.9 0.7 0.9 0.3 0.4 0.5 0.5 0.2 0.9 0.2 0.7 0.6 0.7 0.1
0.9 0.4 0.1 0.5 0.1
0.9 0.7 0.6 0.5 0.2 1 0.8 0.2 0.7 0.2 0.4 0.5 0.2 0.6 0.5 0.4 0.3 0.2
0.5 0 0 0 0.6 0.2 0.3 0.3 0 0 0.2 0.4 0 0.5 0.4 0.3 0.1 0.4
0.5 0.8 0.4 0.8 0 0.5 0.5 0.6 0.5 0.1 0.5 0.1 0.9 0.8 0.2 0.8 0.9 0.3
0.4 0.2 0.3 0.1 0.1 0.8
0.8 0.8 0.8 0.8 0.6 0.8 1 0.2 0.8 0.2 0.5 0.6 0 0 0 0.3 0 0
0.2 0 0 0.8 0.6 0 0 0 0 0 0 0 0.5 0.4 0.5 0.2 0.1 0.4
0.2 0.7 0.7 0.1 0.1 0.4 0.7 0.5 0.8 0.4 0.9 0.7 0.4 0.8 0.1 0.8 0.7 0.5
0.1 0.8 0.9 0 0.2 0.8
0.4 0.2 0.3 0.2 0 0.2 0.2 1 0.2 0 0 0.4 0.2 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0.8 0 0 0.7 0 0.7 0.6 0.5 0.1 0.3
0.8 0.4 0.3 0.4 0.6 0.3 0.9 0.4 0.9 0.9 0.3 0.2 0.5 0.7 0.9 0.3 0.2 0.1
0.7 0.7 0.9 0.8 0.3
0.7 0.8 0.7 0.6 0.7 0.7 0.8 0.2 1 0.2 0.5 0.6 0.3 0.2 0.2 0.6 0.4 0.2
0.4 0.5 0 0.7 0.8 0 0 0 0 0 0.2 0.4 0.8 0.1 0.1 0 0.7 0.2
0.7 0.7 0.8 0.3 0.7 0.7 0.9 0.2 0.2 0.8 0.6 0.2 0.2 0.8 0.5 0.8 0.6 0.7
0.5 0.5 0.2 0.1 0.6 0.1
0.3 0.1 0.2 0.2 0 0.2 0.2 0 0.2 1 0 0 0 0 0 0 0.4 0.6
0.7 0.8 0.3 0.2 0 0.5 0.3 0.5 0 0 0.3 0.8 0.1 0.9 0.3 0.9 0.6 0.1
0.3 0.5 0 0.4 0.5 0.6 0.5 0 0.4 0 0.5 0.8 0.7 0.3 0.5 0.4 0.5 0.5 0
0.3 0.9 0.9 0.3 0.3
0.2 0.2 0.5 0.4 0.5 0.4 0.5 0 0.5 0 1 0.5 0.2 0 0 0.4 0 0.2 0
0 0 0.6 0.4 0 0 0 0 0 0 0 0.9 0.7 0.2 0.5 0.2 0.8 0.5
0.6 0.6 0.1 0.3 0 0.5 0.8 0.6 0.2 0.8 0.8 0.2 0.5 0 0.1 0.1 0.4 0.6
0.9 0.3 0.5 0.2 0.8
0.6 0.3 0.4 0.5 0.6 0.5 0.6 0.4 0.6 0 0.5 1 0 0 0 0.4 0.2 0 0
0 0 0.6 0.3 0 0 0 0 0 0 0 0.2 0.1 0.8 0.3 0.5 0.1 0.2
0.2 0.9 0.5 0.3 0.2 0.9 0.4 0.9 0.7 0.2 0.1 0.4 0.8 0.9 0.5 0.8 0.7 0.2
0.1 0.4 0.4 0.1 0.8
0.2 0 0.2 0.2 0 0.2 0 0.2 0.3 0 0.2 0 1 0 0.4 0.2 0.6 0.2 0
0 0 0 0 0 0.3 0.4 0 0 0.3 0 0.4 0.1 0.5 0.5 0.7 0.6 0.5
0.1 0.3 0.4 0.7 0.6 0.8 0.2 0.5 0.9 0.1 0.2 0.1 0.9 0.1 0.5 0 0.9 0.8
0.3 0.2 0.4 0.1 0.5
0.3 0.2 0.2 0.2 0.2 0.6 0 0 0.2 0 0 0 0 1 0.4 0 0 0.3 0
0 0 0 0 0 0 0 0.6 0 0 0.4 0.7 0.7 0.1 0 0.9 0.2 0.5
0.8 0.6 0.4 0.1 0.2 0.7 0.1 0.7 0.2 0.3 0.8 0.1 0.2 0.1 0.2 0.9 0.9 0.8
0.5 0.5 0.5 0.8 0.3
0.2 0 0 0.3 0.3 0.5 0 0 0.2 0 0 0 0.4 0.4 1 0 0.2 0.3 0
0 0 0 0 0 0.4 0.3 0.2 0 0.4 0.6 0.2 0.8 0.8 0 0.2 0.1 0.6 0
0.6 0.8 0.9 0.6 0.4 0.7 0.7 0.2 0.9 0.2 0.9 0.6 0.2 0.4 0 0.4 0.9 0.1
0.1 0.6 0.1 0.5
0 0.4 0 0.1 0.6 0.4 0.3 0 0.6 0 0.4 0.4 0.2 0 0 1 0 0 0
0 0 0.6 0.5 0 0 0 0 0 0 0 0.8 0.1 0.3 0.1 0.1 0.8 0.6
0.2 0.5 0.2 0.4 0.4 0.9 0.5 0.1 0.9 0.7 0.6 0.9 0.6 0.6 0.6 0.7 0.8 0
0.8 0.6 0.7 0.8 0.2
0 0.5 0 0 0.2 0.3 0 0 0.4 0.4 0 0.2 0.6 0 0.2 0 1 0.5
0.2 0.4 0.3 0 0 0.3 0 0 0 0 0.2 0.4 0.2 0.6 0.2 0.4 0.8 0.7
0.7 0.7 0.1 0.1 0.6 0.7 0.1 0.6 0.7 0.4 0.2 0.2 0.9 0.3 0.5 0.9 0.7 0.8
0.8 0.2 0.4 0.2 0.7 0
0 0.3 0 0 0 0.2 0 0 0.2 0.6 0.2 0 0.2 0.3 0.3 0 0.5 1
0.3 0.4 0.3 0.2 0 0 0 0 0 0 0.2 0.3 0.8 0.9 0.1 0.5 0.4 0.7
0.9 0.1 0.3 0.5 0.4 0.4 0.6 0.6 0.2 0.4 0.8 0.9 0.6 0.7 0.3 0.8 0.3 0.9
0.2 0.6 0.1 0.6 0.5 0.6
0 0.2 0.2 0 0 0.5 0.2 0 0.4 0.7 0 0 0 0 0 0 0.2 0.3 1
0.6 0.2 0.4 0 0 0 0.3 0 0 0.2 0.6 0.9 0.7 0.5 0.1 0.6 0.7 0.5
0.9 0.2 0.7 0.6 0.9 0.3 0.2 0.7 0.1 0.9 0.4 0.1 0.5 0.9 0.8 0.3 0.9 0.8
0.1 0.3 0.2 0.1 0.6
0.1 0.3 0.2 0.4 0 0 0 0 0.5 0.8 0 0 0 0 0 0 0.4 0.4
0.6 1 0 0 0 0.6 0 0 0 0 0.2 0.5 0.4 0.9 0.9 0.9 0.3 0.2
0.8 0.7 0.5 0.9 0.3 0.5 0.9 0.9 0.5 0.6 0.1 0.4 0.8 0.9 0.2 0.9 0.8 0.7 0
0.8 0.7 0.4 0.4 0.9
0 0 0 0.2 0.1 0 0 0 0 0.3 0 0 0 0 0 0 0.3 0.3
0.2 0 1 0.5 0.6 0.2 0.4 0 0 0 0.2 0 0.3 0.6 0.6 0.3 0.5 0.8
0.1 0.3 0.4 0.6 0.4 0.6 0.4 0.1 0.9 0.8 0.9 0.1 0.8 0.3 0.4 0.7 0.2 0.7
0.2 0.3 0.1 0.9 0.1 0.1
0.5 0.8 0.8 0.8 0.8 0 0.8 0 0.7 0.2 0.6 0.6 0 0 0 0.6 0 0.2
0.4 0 0.5 1 0.6 0 0.4 0.3 0 0 0 0 0.7 0.9 0.1 0.6 0.9 0.1
0.5 0.4 0.8 0.9 0.9 0.2 0.3 0.4 0.8 0.6 0.9 0.8 0.4 0.9 0.4 0.6 0.8 0.9
0.3 0.1 0.2 0.1 0.1 0.8
0.2 0.4 0 0.2 0.4 0.6 0.6 0 0.8 0 0.4 0.3 0 0 0 0.5 0 0 0
0 0.6 0.6 1 0 0 0 0 0 0 0 0.9 0.2 0 0.4 0.4 0.9 0
0.4 0.3 0.5 0.4 0.8 0.4 0.6 0.4 0.4 0.8 0.3 0.9 0.5 0.8 0.2 0.4 0.6 0.6
0.1 0.6 0.9 0.6 0.3
0.3 0.3 0 0.2 0.3 0.2 0 0 0 0.5 0 0 0 0 0 0 0.3 0 0
0.6 0.2 0 0 1 0.3 0.2 0 0 0.5 0.8 0.4 0 0.9 0 0.2 0.4 0.9
0.2 0.3 0.1 0.8 0.2 0.7 0.5 0.4 0.9 0.1 0.5 0 0.3 0 0.6 0.9 0.8 0.9
0.9 0.2 0.1 0.9 0.7 0.7
0.6 0 0.5 0.1 0.5 0.5 0.9 0.9 0.6 0.5 0.8 0.2 0.1 0.3 0.9 0.7 0.2 0.8
0.9 0.1 0.9 0.9 0.8 0.1 0.9 0.5 0.9 0.1 0.9 0.8 0.5 0 0.2 0 0.6 0.2
0.2 0.9 0.2 0.5 0.8 0.2 0.2 0.4 0.8 0.3 1 0.3 0.5 0.5 0.6 0.4 0.8 0.2
0.1 0.2 0.7 0.9 0.4 0.1
0.8 0 0.2 0.1 0.5 0.1 0.7 0.3 0.2 0.8 0.8 0.1 0.2 0.8 0.2 0.6 0.2 0.9
0.4 0.4 0.1 0.8 0.3 0.5 0.6 0.8 0.4 0.6 0.9 0.6 0.9 0.9 0 0.1 0.1 0.9 0
0.2 0.9 0.5 0.2 0.1 0.5 0.4 0.2 0.3 0.3 1 0.6 0.6 0 0.9 0.7 0.2 0.1
0.5 0.3 0.7 0.4 0.9
0.4 0.8 0.6 0.6 0.2 0.9 0.4 0.2 0.2 0.7 0.2 0.4 0.1 0.1 0.9 0.9 0.9 0.6
0.1 0.8 0.8 0.4 0.9 0 0.6 0.3 0.7 0.7 0.6 0.1 0.5 0.4 0.9 0.6 0.2 0.5
0.1 0.5 0.9 0.7 0.5 0.6 0.5 0.3 0.2 0.4 0.5 0.6 1 0.2 0.9 0.1 0.1 0.5
0.2 0.1 0.1 0.9 0.3 0.7
0.5 0.2 0.6 0.7 0.9 0.8 0.8 0.5 0.8 0.3 0.5 0.8 0.9 0.2 0.6 0.6 0.3 0.7
0.5 0.9 0.3 0.9 0.5 0.3 0 0.5 0.2 0.9 0.4 0.5 0.2 0.8 0.6 0.6 0.7 0.3
0.5 0.2 0.3 0.1 0.8 0.7 0.8 0.8 0.3 0.9 0.5 0.6 0.2 1 0.4 0.8 0.7 0.5
0.8 0 0.5 0.5 0 0.6
0.3 0 0.6 0.6 0.2 0.2 0.1 0.7 0.5 0.5 0 0.9 0.1 0.1 0.2 0.6 0.5 0.3
0.9 0.2 0.4 0.4 0.8 0 0.8 0.6 0.6 0.6 0 0.2 0.3 0.7 0.5 0.3 0.8 0.8
0.7 0.6 0.5 0.2 0.7 0.7 0.9 0.1 0.2 0.9 0.6 0 0.9 0.4 1 0.3 0.5 0.2
0.9 0 0.1 0.2 0.2 0.2
0.6 0.4 0.5 0.1 0.7 0.8 0.8 0.9 0.8 0.4 0.1 0.5 0.5 0.2 0.4 0.6 0.9 0.8
0.8 0.9 0.7 0.6 0.2 0.6 0.5 0.1 0.5 0.4 0.7 0.2 0.8 0.3 0.9 0.7 0.1 0.1
0.1 0.1 0.7 0.9 0.2 0.4 0.9 0.9 0.4 0.9 0.4 0.9 0.1 0.8 0.3 1 0.1 0.7
0.5 0.7 0.9 0.6 0.5 0.9
0.7 0.8 0.7 0.4 0.6 0.9 0.7 0.3 0.6 0.5 0.1 0.8 0 0.9 0 0.7 0.7 0.3
0.3 0.8 0.2 0.8 0.4 0.9 0.7 0.9 0.9 0.4 0.9 0.2 0.8 0.4 0 0.3 0.6 0.9
0.3 0.7 0.2 0.1 0 0.7 0.9 0.9 0.7 0.9 0.8 0.7 0.1 0.7 0.5 0.1 1 0.3
0.4 0.9 0.3 0.4 0.8 0.3
0.7 0.3 0.3 0.4 0.7 0.3 0.5 0.2 0.7 0.5 0.4 0.7 0.9 0.9 0.4 0.8 0.8 0.9
0.9 0.7 0.7 0.9 0.6 0.8 0.1 0.1 0 0.3 0.2 0.4 0.4 0.9 0.2 0.4 0.9 0.9
0.8 0.7 0.2 0.9 0 0.7 0.8 0.3 0.7 0.1 0.2 0.2 0.5 0.5 0.2 0.7 0.3 1
0.8 0.1 0.3 0.1 0.6 0.9
0.6 0.2 0.2 0.7 0.1 0.4 0.1 0.1 0.5 0 0.6 0.2 0.8 0.8 0.9 0 0.8 0.2
0.8 0 0.2 0.3 0.6 0.9 0.7 0.3 0.7 0.8 0.9 0.1 0.1 0.9 0.6 0.5 0.2 0.5
0.9 0.2 0.5 0.8 0.5 0.2 0.6 0.8 0.8 0.9 0.1 0.1 0.2 0.8 0.9 0.5 0.4 0.8 1
0.4 0.6 0.8 0.2 0.9
0.4 0.3 0.8 0.1 0.9 0.2 0.8 0.7 0.5 0.3 0.9 0.1 0.3 0.5 0.1 0.8 0.2 0.6
0.1 0.8 0.3 0.1 0.1 0.5 0.4 0.7 0.1 0.5 0.9 0.5 0.6 0 0.4 0.2 0.4 0.7
0.8 0.2 0.6 0.2 0.5 0.2 0.2 0 0.2 0.2 0.2 0.5 0.1 0 0 0.7 0.9 0.1
0.4 1 0.7 0.9 0.6 0.1
0.6 0.7 0.2 0.7 0.4 0.3 0.9 0.7 0.2 0.9 0.3 0.4 0.2 0.5 0.1 0.6 0.4 0.1
0.3 0.7 0.1 0.2 0.6 0.2 0.1 0.9 0.8 0.3 0.6 0.5 0.3 0.7 0.8 0.4 0.1 0.8
0.3 0.4 0.6 0.3 0.8 0.5 0.6 0.2 0.7 0.1 0.7 0.3 0.1 0.5 0.1 0.9 0.3 0.3
81
0.6 0.7 1 0.5 0.1 0.9
0.3 0.9 0.8 0.8 0.1 0.1 0 0.9 0.1 0.9 0.5 0.4 0.4 0.5 0.6 0.7 0.2 0.6
0.2 0.4 0.9 0.1 0.9 0.2 0.4 0.1 0.4 0.4 0.8 0.6 0.2 0.3 0.7 0.5 0.4 0.1
0.9 0.2 0.4 0.2 0.3 0.2 0.3 0.2 0.1 0.9 0.9 0.7 0.9 0.5 0.2 0.6 0.4 0.1
0.8 0.9 0.5 1 0.5 0.8
0.9 0.4 0.5 0.4 0.5 0.1 0.2 0.8 0.6 0.3 0.2 0.1 0.1 0.8 0.1 0.8 0.7 0.5
0.1 0.4 0.1 0.1 0.6 0.1 0.4 0.3 0.2 0.7 0.4 0.4 0.3 0.1 0.3 0.6 0.4 0.3
0.8 0.6 0.2 0.2 0.1 0.6 0.2 0.6 0.9 0.7 0.4 0.4 0.3 0 0.2 0.5 0.8 0.6
0.2 0.6 0.1 0.5 1 0.7
0.7 0.7 0.9 0.7 0.1 0.8 0.8 0.3 0.1 0.3 0.8 0.8 0.5 0.3 0.5 0.2 0 0.6
0.6 0.9 0.1 0.8 0.3 0.9 0.9 0.9 0.9 0.9 0.2 0.2 0.7 0.2 0.7 0.5 0.9 0.1
0.1 0.8 0.3 0.7 0.6 0.2 0.9 0 0.1 0.7 0.1 0.9 0.7 0.6 0.2 0.9 0.3 0.9
0.9 0.1 0.9 0.8 0.7 1
];

N=size(cs,1);
N1=10;
N2=30;

    rx=[37,45,78,3,1,80,58,74,89,41,55,46,38,13,84,36,44,65,30,20,5,4
8,75,64,90,22,27,79,56,86,33,16,53,52,87,68,28,18,50,54,40,15,71,31,6
3,73,23,69,14,19,42,81,49,60,76,34,72,10,8,25];
    ry=[46,26,48,33,12,38,7,89,62,8,1,77,20,14,56,37,66,73,10,81,67,2
1,6,36,55,61,58,35,27,64,43,28,25,45,34,87,78,75,13,16,63,9,39,19,47,
5,30,42,52,68,70,86,71,41,79,76,32,51,31,15];
    cs10=cs(1:N1,1:N1);
    [g1,u]=groupdemo(cs10);
    draw(g1,rx,ry);
    cs30=cs(1:N2,1:N2);
    [g2,u]=groupdemo(cs30);
    draw(g2,rx,ry);
    [g3,u]=groupdemo(cs);
    draw(g3,rx,ry);

    for i=10:30
        css=cs(1:i,1:i);
        [g,umean(i),uvar(i)]=groupdemo(css);
    end

    figure;
    plot(10:30,umean(10:30),'o-');
    title('均值 10~30');