雅克比迭代法求解线性方程组

#include<bits/stdc++.h>
using namespace std;
const int maxn=110;
double a[110][110],b[110];
double x[110],tmpx[110];
const double eps=1e-6;
int main(){
	double epsilon,sum;
	int k=0;
	int n;scanf("%d",&n);
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			scanf("%lf",&a[i][j]);
	for(int i=0;i<n;i++) 
		scanf("%lf",&b[i])
	while(1){
		epsilon=0;
		k++;
		for(int i=0;i<n;i++){
			sum=0;
			for(int j=0;j<n;j++){
				if(i==j) continue;
				sum+=a[i][j]*x[j];
			}
			tmpx[i]=(b[i]-sum)/a[i][i];
			epsilon=max(epsilon,fabs(tmpx[i]-x[i]));
		}
		for(int i=0;i<n;i++) x[i]=tmpx[i];
		if(epsilon<eps){
			printf("迭代次数为%d\n",k);
			printf("方程组的解为：\n");
			for(int i=0;i<n;i++)
				printf("x%d=%.6f\n",i,x[i]); 
			return 0;
		}
		if(k>=100){
			puts("-1");
			return 0;
		}
	}
	
	return 0;
}
/*
4
4 0.24 -0.08 
0.09 3 -0.15
0.04 -0.08 4
8 9 20
*/