矩阵求逆 代码

* 矩阵求逆函数(Metrix = Metrix(-1))，参数说明:
*   Metrix - 原始矩阵，也是求逆之后的矩阵，以一维数组形式存储，n行n列

void Inverse(double *Metrix, int n)
{
int i, j, k;
    for(k=0; k<n; k++)
{
   for(i=0; i<n; i++)
   {
    if(i != k)
     Metrix[i*n+k] = - Metrix[i*n+k] / Metrix[k*n+k];
   }
   Metrix[k*n+k] = 1/Metrix[k*n+k];
  
   for(i=0; i<n; i++)
   {
    if(i != k)
    {
     for(j=0; j<n; j++)
     {
      if(j != k)
       Metrix[i*n+j] += Metrix[k*n+j] * Metrix[i*n+k];
     }
    }
   }
  
   for(j=0; j<n; j++)
   {
    if(j != k)
     Metrix[k*n+j] *= Metrix[k*n+k];
   }
}
}