利用克拉默法则求X=B／A

#include 
#include 
#include 
using  namespace std;
#define N 50
//得到给定矩阵src的逆矩阵保存到des中。
double getA(double arcs[N][N],int n);
//计算每一行每一列的每个元素所对应的余子式，组成A*
void  getAStart(double arcs[N][N],int n,double ans[N][N])
{
	if(n==1)
	{
		ans[0][0] = 1;
		return;
	}
	int i,j,k,t;
	double temp[N][N];
	for(i=0;i=i?k+1:k][t>=j?t+1:t];
				}
			}


			ans[j][i]  =  getA(temp,n-1);
			if((i+j)%2 == 1)
			{
				ans[j][i] = - ans[j][i];
			}
		}
	}
}
bool GetMatrixInverse(double src[N][N],int n,double des[N][N])
{
	double flag=getA(src,n);
	double t[N][N];
	if(flag==0)
	{
		return false;
	}
	else
	{
		getAStart(src,n,t);
		for(int i=0;i=i)?k+1:k];

			}
		}
		double t = getA(temp,n-1);
		if(i%2==0)
		{
			ans += arcs[0][i]*t;
		}
		else
		{
			ans -=  arcs[0][i]*t;
		}
	}
	return ans;
}
int main(){
	double a[N][N];
	double b[N][N];
	double c[N][N];
	double d[N][N];
	int n,m,q;
	int i,j,k,kk=0;
	cout<<"陈扬送给给位没有matlab又深陷线性代数的道友^-^"<>m>>n;
		 for(i=0;i>d[i][j];
		GetMatrixInverse(d, n, a);
		   cout<<"输入矩阵b的大小n，q："<>n>>q;
		   for(i=0;i>b[i][j];
		for(i=0;irun：
 
  
陈扬送给给位没有matlab又深陷线性代数的道友^-^
 
  
求X=B／A的程序，利用克拉默法则：
 
  
输入行列式a的大小m,n：
 
  
5 5
 
  
    0     1     1     1     1
 
  
     1     0     1     1     1
 
  
     1     1     0     1     1
 
  
     1     1     1     0     1
 
  
     1     1     1     1     0
 
  
输入矩阵b的大小n，q：
 
  
5 1
 
  
1  2 3 4 5
 
  
矩阵c为矩阵b/a的乘积, 相当于c为克拉默法则运算后的X的解的矩阵
 
  
   2.75
 
  
   1.75
 
  
   0.75
 
  
  -0.25
 
  
  -1.25
 
  
|A|=4