只要矩阵可逆,用程序算的方法与该矩阵具体性质没关系.
a[]={3,4 ,2,1,5 ,
0 ,0 ,3 ,1 ,3 ,
0 ,0 ,2 ,5 ,7 ,
0,0,0,1,3 ,
0 ,0 ,0 ,2 ,2 };
int n = 5;
int MatrixInvert(float a[],int n)//实矩阵求逆
{
int *is,*js,i,j,k,l,u,v;
float d,p;
is=(int*)malloc(n*sizeof(int));
js=(int*)malloc(n*sizeof(int));
for (k=0; k