#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;
}