function C=matrixmul(A,B)
%此程序用于實(shí)現(xiàn)任意兩個(gè)規(guī)模為m*n和n*m的矩陣相乘
%格式為C=matrixmul(A,B)
%A和B分別為兩個(gè)任意階矩陣，但要確保階數(shù)符合矩陣相乘的規(guī)則
[m,n]=size(A);
[c,d]=size(B);
if (c==n)==1
    for i=1:m
        for j=1:m
            temp=0;
            for k=1:n
                temp=temp+A(i,k)*B(k,j);
            end
            C(i,j)=temp;
        end
    end
    if m>=n
        C; 
    else
        C=C(1:m,1:d);
    end
else
    error('Inner matrix dimensions must agree.')
end