c语言——求逆矩阵，伴随矩阵，行列式

算法思想：利用完全展开式计算行列式|A|;

                递归计算n阶行列式；终止条件n=2：|A|元素交叉相乘相减

               公式： A-1=A*/|A|

要点：要动态开辟二维空间来储存降阶矩阵

float determinant(float **c,int n) {//按第一行完全展开式计算|A|
	float det,t,**temp1;
	temp1=(float **)malloc(sizeof(float *)*(n-1));//为储存降阶矩阵开辟空间
	for(int p=0; p=i)?k+1:k];//若列数小于i则不变;若列数大于等于i则向后移动一列;从而记录余子式的矩阵
				}
			}
			t=determinant(temp1,n-1);//递归计算
			if(i%2==0) {//判断余子式的正负;与第一行元素相乘;相加得行列式
				det+=c[0][i]*t;
			} else {
				det-=c[0][i]*t;
			}
		}
	}
	return det;
}
float **adjoint(float **c,int n) {//计算每一行每一列的每个元素所对应的余子式，组成A*
	float **temp2,**adj;
	temp2=(float **)malloc(sizeof(float *)*(n-1));         //为n-1阶矩阵开辟空间
	for(int p=0; p=i)?k+1:k][(t>=j)?t+1:t];//剔除元素所在行与列之后的矩阵
				}
			}
			adj[j][i]=determinant(temp2,n-1);//计算代数余子式Aji为  转置  后的
			if((i+j)%2==1) {//判断符号(-1)^(i+j)
				adj[j][i]=-adj[j][i];
			}
		}
	}
	return adj;
}

在输出时加一句判断|A|==0；是否可逆。

bt=adjoint(a,n);
z=determinant(a,n);//返回行列式的值
if(z==0) { 
    printf("n矩阵不可逆！！");
} else {
    printf("逆矩阵为：n");
    for(int i=0; i