clc,clear;
disp('-----層次分析系統-----');
xx=input('請輸入層數:');
A(1,1)=1;
shang=1;
for j=1:(xx-1)  
    disp('請輸入第'),disp(j+1),disp('層的元素個數:');
    z=input('');
     CCR=1;
     while (CCR>0.1)
     n=1;
     for i=1:shang
     CR(i)=1;
      while (CR(i)>0.1)
         w=input('請輸入合適的判斷比較矩陣 :')
         [V,D]=eig(w);
         x=max(max(D));
         [y yy]=find(D==x); 
         [k,kk]=size(w);
         CI=(x-k)/(k-1);
         switch k
         case 3
          RI=0.58;
         case 4
          RI=0.90;
         case 5
          RI=1.12;
         case 6
          RI=1.24;
         case 7
          RI=1.32;
         case 8
          RI=1.41;
         case 9
          RI=1.45;
         end
       if k==1|k==2 CR(i)=0;
           RI=0;
         else CR(i)=CI/RI;
       end
     end
      x
      CI
      RI
      CR(i)
      ss=sum(V(:,y));
      V(:,y)=V(:,y)/ss;
      M=V(:,y)
      tzxl(n:(n+k-1))=V(:,y).*A(j,i);
      n=n+k;
  end
   CR
   CCR=sum(CR.*A(j,1:shang))
end
    A(j+1,1:z)=tzxl./sum(tzxl);
    shang=z;
end
A