高斯消元线性方程组

#include <cmath>
#include <stdio.h>

int main()
{
	int n,i,j,k;
	double a[100][100],b[100],x;
	printf("请输入未知数个数：");
	scanf("%d",&n);
	printf("请输入%d个方程组(未知数系数及等号右边常数):\n",n);
	for (i = 1;i <= n;i++)
		for (j = 1;j <= n+1;j++)
			scanf("%lf",&a[i][j]);

	for (i=1;i<=n;i++)
		for (j=i+1;j<=n;j++)
			if (fabs(a[j][i])>1e-7)
			{
				x=a[i][i]/a[j][i];
				for (k=i;k<=n+1;k++)
				a[j][k]=a[j][k]*x-a[i][k];
			}
	for (i = n;i > 0;i--)
	{
		b[i] = a[i][n+1]/a[i][i];
		for (j = i-1;j > 0;j--)
			a[j][n+1]=a[j][n+1]-b[i]*a[j][i];
	}
	printf("解得:\n");
	for (i=1;i<=n;i++)
    {
         printf("x%d = %lf\n",i,b[i]);
    }
	return 0;
}