矩阵求逆

矩阵求逆过程中对零判断条件很重要啊！！！！！！！！！！！！！！

Matrix matrixInvese(Matrix mat)
{
if(mat.n==0 || mat.n!=mat.m)
{
printf("不能求逆\n");
exit(0);
}
Matrix temp;
temp.m=mat.m;
temp.n=mat.n;
temp.ary=(double *)malloc(sizeof(double)*temp.m*temp.n);
int i;
int j;
for(i=0; i<mat.n; i++)
{
for(j=0; j<mat.n; j++)
{
if(i==j)
{
temp.ary[i*temp.n+j]=1;
}
else
{
temp.ary[i*temp.n+j]=0;
}
}
}
for(i=0; i<mat.n; i++)
{
for(j=0; j<mat.n; j++)
{
if(j!=i)
{
int k;
double temp1;
temp1=mat.ary[i*mat.n+i]/mat.ary[j*mat.n+i];
for(k=0; k<mat.n; k++)
{
mat.ary[j*mat.n+k]=mat.ary[j*mat.n+k]*temp1-mat.ary[i*mat.n+k];
temp.ary[j*temp.n+k]=temp.ary[j*temp.n+k]*temp1-temp.ary[i*temp.n+k];
}
}
}
}
for(i=0; i<mat.n; i++)
{
for(j=0; j<mat.n; j++)
{
temp.ary[i*temp.n+j]/=mat.ary[i*mat.n+i];
}
}
return temp;
}