我要编的是用简单迭代法德迭代格式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();
 }