sor_fun.h

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

/*头文件sor_fun.h*/

#include<math.h>
#include<stdio.h>
int SOR(int n,float w,double a[],double b[],double x[],double eps)
{
	int i,j;
	int k=0;
	double p,t,s,q;
	p=eps+1.0;
	while(p>=eps)
	{
		k=k+1;
		p=0.0;
		for(i=0;i<=n-1;i++)
		{
			t=x[i];  s=0.0;
			for(j=0;j<=n-1;j++)
				s=s+a[i*n+j]*x[j];
			x[i]=t+w*(b[i]-s)/a[i*n+i];  /*用松弛因子w校正*/
			q=fabs(x[i]-t)/(1.0+fabs(x[i]));
			if(q>p)  p=q;
		}
		printf("第%d次迭代结果:\n",k);
		for(i=0;i<=n-1;i++)
			printf("x[%d]=%2.2f\n",i,x[i]);
	}
	return(1);
}