jacobi_fun.h

来自「包含冒泡排序算法函数模板和数据结构中的字符串查找例子及解线性方程组的3个典型算法」· C头文件 代码 · 共 30 行

/*头文件jacobi_fun.h*/
#include<math.h>
#include<stdio.h>
#define N 10
int jacobi(int n,double a[],double b[],double x[],double eps)
{
	int i,j;
	int k=0;
	double p,s,q;
	double t[N];
	for(i=0;i<N;i++) t[i]=0.0;
	p=eps+1.0;
	while(p>=eps) /*不满足精度要求，继续迭代*/
	{   k++;
		p=0.0;
		for(i=0;i<=n-1;i++)
		{
			s=0.0;
			for(j=0;j<=n-1;j++)
				if(j!=i) s=s+a[i*n+j]*t[j];
			x[i]=(b[i]-s)/a[i*n+i];
			q=fabs(x[i]-t[i])/(1.0+fabs(x[i]));
			if(q>p) p=q;
		}
		printf("第%d次迭代:\n",k);
		for(i=0;i<=n-1;i++)
		{t[i]=x[i]; printf("x[%d]=%3.3e\n",i,x[i]);}
	}
	return(1);/*正常返回*/
}