define:
Martix(A)=|a11|a12|a13|
-------------- |a21|a22|a23|
---------------|a31|a32|a33|
C1=(a11,a21,a31)^t
…Cn=(a11,a21,…,an1)^t
Martix(B)=|b11|b12|b13|
-------------- |b21|b22|b23|
---------------|b31|b32|b33|
R1=(b11,b12,b13)
…
solve:
CnRn=|a11 b11| a11b12|…a11b1n|
----------|a21b11|…a21b1n|
…
----------|an1b11|…an1b1n|
matAmatB=sum(matA.CnmatB.Rn)