2011软考程序员算法实例解析：矩阵求逆算法

 　/**

　　* 求矩阵A的逆矩阵Ai

　　*@param A 源矩阵

　　*@param Ai 逆矩阵

　　*@param size 矩阵的大小

　　*@return 求解成功返回非零值，失败返回零

　　*/

　　int InverseMatrix(double **Ai, double **A, int size)

　　{

　　int i, j;

　　double *b, *x;

　　b = (double*)malloc(sizeof(double)*size);

　　x = (double*)malloc(sizeof(double)*size);

　　for (i = 0; i < size; i++)

　　{

　　memset(b, 0, sizeof(double)*size);

　　b[i] = 1;

　　if (!LinearEquation(A, x, b, size))

　　{

　　free(b);

　　free(x);

　　return 0;

　　}

　　for (j = 0; j < size; j++)

　　Ai[j][i] = x[j];

　　}

　　free(b);

　　free(x);

　　return 1;

　　}