📄 jacobi_fun.h
字号:
/*头文件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);/*正常返回*/
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -