,…,n计算
3 如果则输出结果;否则执行4
4 如果则不收敛,终止程序;否则,转2
1.2.2 程 序
#include
#include
#define eps 1e-6
#define max 100
void jacobi(float *a,int n,float x[])
{
int i,j,k=0;
float epsilon,s;
float *y= new float [n];
for(i=0;i while(1)
{
epsilon=0;
k++;
for(i=0;i {
s=0;
for(j=0;j {
if(j==i)continue;
s+=*(a+i*(n+1)+j)*x[j];
}
y[i]=(*(a+i*(n+1)+n)-s)/(*(a+i*(n+1)+i));
epsilon+=fabs(y[i]-x[i]);
}
for(i=0;i if(epsilon {printf("die dai ci shu wei:%d\n",k);return;}
if(k>=max)
{printf("die dai fa san");return;}
}
delete y;
}
void main()
{s
int i;
float a[4][5]={10,-1,2,0,-11,0,8,-1,3,-11,2,-1,10,0,6,-1,3,-1,11,25};
float x[4];
jacobi(a[0],4,x);
for(i=0;i<4;i++)printf("x[%d]=%f\n",i,x[i]);
}
1.2.3 输出结果
1.2.4 结果分析
迭代次数增加时,精度越高。从输出结果可以看出此方程组的迭代次数为17,迭代结果越来越接近精确解了,于是
=-1.467391, =-2.358696, =0.657609,
上一页 [1] [2] [3] [4] [5] [6] 下一页