Matlab中矩阵的乘法C++代码

void Matrix_Multip(double*A_matrix,double *B_matrix,double*AB_matrix,int A_Rows,int A_Colunms,int B_Rows,int B_Colunms)
{
 //矩阵乘法 double*AB_matrix最后的结果
 for (int i=0;i<A_Rows;i++)
 {
  for (int j=0;j<B_Colunms;j++)
  {
   AB_matrix[i*B_Colunms+j]=0;

   for (int k=0;k<A_Colunms;k++)
   {
    AB_matrix[i*B_Colunms+j]=AB_matrix[i*B_Colunms+j]+A_matrix[i*A_Colunms+k]*B_matrix[k*B_Colunms+j];
   }   
  }
 }
}