(求助)简单迭代格式编程
我要编的是用简单迭代法德迭代格式1解方程组(三元一次),运行不出结果。。。请各位前辈帮忙看看迭代公式X(k+1)=CX(k)+B
#include <stdio.h>
#include<math.h>
void main()
{
float a[3][3]={{20,2,3},{1,8,1},{2,-3,15}},b[3]={24,12,30};
int I[3][3]={{1,0,0},{0,1,0},{0,0,1}};
float c[3][3],x[3]={0 }, x1[3];
int i,j,k;
float m;
for(i=0;i<3;i++)
for(j=0;j<3;j++)
{
c[i][j]=I[i][j]-a[i][j];
}
printf("c=\n");
for(i=0;i<3;i++){
for(j=0;j<3;j++)
printf("%10f",c[i][j]);
printf("\n");}
getch();
for(k=0;k<500;k++)
{
for(i=0;i<3;i++)
{m=0;
for(j=0;j<3;j++)
{
m=m+c[i][j]*x[j];
}
x1[i]=m+b[i];
}
while(fabs(x1[1]-x[1])<0.00001&&fabs(x1[2]-x[2])<0.00001&&fabs(x1[3]-x[3])<0.00001) break;
for(i=0;i<3;i++) x[i]=x1[i];
}
for(i=0;i<3;i++)
printf("x[%d]=%10f\n",i,x[i]);
getch();
}