#include<stdio.h>
void gaussianElinination (float a[][4],float b[],int n){
int i,j,k;
float factor;
for(i=0;i<n-1;i++){
for(j=i+1;j<n;j++){
factor=a[j][i]/a[i][i];
for (k=i;k<=n;k++){
a[j][k]-=factor*a[i][k];
}
b[j]-=factor*b[i];
}
}
for(i=n-1;i>=0;i--){
b[i]/=a[i][i];
a[i][i]/=a[i][i];
for (j=i-1;j>=0;j--){
b[j]-=a[j][i]*b[i];
a[j][i]/=a[i][i];
}
}
}
int main(){
float a[3][4]={{1,2,2,1},
{2,5,3,0},
{1,4,3,0}};
float b[3]={1,0,0};
int n=3;
int i;
gaussianElinination(a,b,n);
printf("Solution:\n");
for (i=0;i<n;i++){
printf("x%d=%f\n",i+1,b[i]);
}
printf("王佳康2023151420");
return 0;
}