#include<stdio.h》
#include<math.h>
void main(){
FILE *f;
double a[20][20],b[20],x[20];
double s,m;
int i,j,k,n;
f=fopen("D:\\实验3高斯数据.txt","r");
fscanf(f,"%d",&n);
for(i=1;i<=10;i++){
for(j=1;j<=10;j++)
fscanf(f,"%lf",&a[i][j]);
fscanf(f,"%lf",&b[i]);
}
fclose(f);
for(k=1;k<=n;k++)
for(i=k+1;i<=n;i++){
for(j=k+1;j<=n;j++){
m=a[i][k]/a[k][k];
a[i][j]=a[i][j]-m*a[k][j];
}
b[i]=b[i]-m*b[k];
}
//迭代
x[n]=b[n]/a[n][n];
for(i=n-1;i>0;i--){
s=0;
for(j=i+1;j<=n;j++)
s=s+a[i][j]*x[j];
x[i]=(b[i]-s)/a[i][i];
}
for(i=1;i<=10;i++){
printf("x[%d]=%0.f\n",i,x[i]);
}
}
#include<math.h>
void main(){
FILE *f;
double a[20][20],b[20],x[20];
double s,m;
int i,j,k,n;
f=fopen("D:\\实验3高斯数据.txt","r");
fscanf(f,"%d",&n);
for(i=1;i<=10;i++){
for(j=1;j<=10;j++)
fscanf(f,"%lf",&a[i][j]);
fscanf(f,"%lf",&b[i]);
}
fclose(f);
for(k=1;k<=n;k++)
for(i=k+1;i<=n;i++){
for(j=k+1;j<=n;j++){
m=a[i][k]/a[k][k];
a[i][j]=a[i][j]-m*a[k][j];
}
b[i]=b[i]-m*b[k];
}
//迭代
x[n]=b[n]/a[n][n];
for(i=n-1;i>0;i--){
s=0;
for(j=i+1;j<=n;j++)
s=s+a[i][j]*x[j];
x[i]=(b[i]-s)/a[i][i];
}
for(i=1;i<=10;i++){
printf("x[%d]=%0.f\n",i,x[i]);
}
}