Matlab求矩阵的逆矩阵的C++代码

void Inverse_Matrix(double * Ora_Matrix,int m_Image_Rows,int m_Image_Colunms)
{
 //求矩阵的逆矩阵

 int rows=m_Image_Rows;
 int colunms=m_Image_Colunms;
 double *Re_Matrix=new double[2*rows*colunms];

 for (int i=0;i<rows;i++)
 {
  for (int j=0;j<colunms;j++)
  {
   Re_Matrix[i*2*colunms+j]=Ora_Matrix[i*colunms+j];
  }
 }

 for (int i=0;i<rows;i++)
 {
  //在原始矩阵右侧加入同大小的单位矩阵
  for (int j=colunms;j<colunms*2;j++)
  {
   Re_Matrix[i*2*colunms+j]=(j-i==colunms)?1:0;
  }
 }

 //进行初等变换

 for (int i=0;i<rows;i++)
 {
  if (Re_Matrix[i*2*colunms+i]!=1.0)
  {
   double tem=Re_Matrix[i*2*colunms+i];

   Re_Matrix[i*2*colunms+i]=1.0;