function Auc=multi1(Aeec)               % 函數名與輸入、輸出變量名multil.m、Aeec、uc

%
      [ ps,pq] = size(Aeec);            % 求出輸入的行和列數
       pn = ps/2;                       % 求出兩輸入變量的個數
       pe = [ eye(pn) zeros(pn)]*Aeec;  % 制作兩兩相乘的矩陣
       pec = [ zeros(pn) eye(pn)]*Aeec;
       pec = pec';
 pa = pe(:,1)*pec(1,:);                 % 計算矩陣的乘積
 Auc = pa(1,:);                         % 制作Auc矩陣
 for i = 2:pn
     Auc = [Auc pa(i,:)];
 end
 
 for j = 2:pq                           %  制作Aucj矩陣
        pa = pe(:,j)*pec(j,:);
        Aucj = pa(1,:);
      for i = 2:pn
          Aucj = [Aucj pa(i,:)];
      end
    Auc =[Auc;Aucj];
 end    
 
 Auc = Auc';                            % 最后的結果