#include<stdio.h>
#include<math.h>
 /**********迭代求根*********/
 double jacobi(double *x)
   {
        double x_new[3];

        double error;
        int i;

          error=0.0;

            x_new[0]=0.125*(20+3*x[1]-2*x[2]);
            error=error+fabs(x[0]-x_new[0]);

            x_new[1]=(33-4*x[0]+x[2])/11;
            error=error+fabs(x[1]-x_new[1]);

            x_new[2]=(12-2*x[0]-x[1])/4;
            error=error+fabs(x[2]-x_new[2]);

       for(i=0;i<3;i++)
       *(x+i)=*(x_new+i);

      return error;
   }

/**************main**************/

   main()
  {
   double x[3];
   double tol;
   int i;
   int k,tmax;
      tmax=1000;
      k=0;
      tol=1e-6;
   for(i=0;i<3;i++)
       x[i]=0.0;
   while(jacobi(x)>=tol&&k<tmax)
     {
            k++;

            printf("changing :  x1=%f  x2=%f  x3= %f",x[0],x[1],x[2]);
            printf("   %d  times\n",k);
     }
      printf("\n Roots are:  x1=%f  x2=%f  x3=%f\n\n",x[0],x[1],x[2]);

     return 0;
  }