[c]代码库#include
#include
#define MAX_N 20
main()
{
int i,j,k;
int n;
int i_k; //列主元所在行的指标
double m; //列主元的值
double a,b; //换行时做终结制值
double D=1; //行列式的值
double A[MAX_N][MAX_N];
double B[MAX_N];
printf("请输入系数矩阵的阶数n:\n");
scanf("%d",&n);
printf("请输入系数矩阵A:\n");
for(i=0;i
for(j=0;j
scanf("%lf",&A[i][j]);
printf("请输入右端项b:\n");
for(i=0;i
scanf("%lf",&B[i]);
for(k=0;k
{
m=0; //开始按列选主元
i_k=k;
for(i=k;i
{
if(fabs(A[i][k])>m)
{
m=fabs(A[i][k]);
i_k=i; //记录主元所在行的指标
}
}
if(m==0)
{
printf("系数矩阵奇异!\ndetA=0\n");
break;
}
else if(i_k!=k)
{
for(j=k;j
{
a=A[i_k][j];
A[i_k][j]=A[k][j];
A[k][j]=a;
}
b=B[i_k];
B[i_k]=B[k];
B[k]=b;
D=(-1)*D; //换行,行列式变号
}
D=D*A[k][k];
for(i=k+1;i
{
A[i][k]=A[i][k]/A[k][k];
}
for(i=k+1;i
{
for(j=k+1;j
{
A[i][j]=A[i][j]-A[i][k]*A[k][j];
}
B[i]=B[i]-A[i][k]*B[k];
}
}
for(i=n-1;i>=0;i--) //回代;临时用b来存储这里的和
{
b=0;
for(j=i+1;j
{b+=A[i][j]*B[j];}
B[i]=(B[i]-b)/A[i][i];
}
D=D*A[n-1][n-1]; //在k进行循环的时候只循环到n-2(也可以在上面让k循环到n-1)
printf("系数矩阵行列式的值为%lf。\n",D);
printf("解为x_i=\n");
for(i=0;i
printf("%lf\n",B[i]);
return 0;
}