书本74面方程组求解

#include <stdio.h>
#include <stdlib.h>
#include<math.h>

int main() {
    double x[3]= {0},y[3];
    double fc[3][3]= {{8,-3,2},{4,11,-1},{6,3,12}};
    double g[3]= {20,33,36};
    double eps=0.000001;
    double temp;
    int flag[3]= {1};
    int i;
    do {
        temp=0;
        for(i=0; i<3; i++) {
            flag[i]=0;
            y[i]=x[i];
            x[i]=(g[i]-fc[i][0]*x[0]*flag[0]-fc[i][1]*x[1]*flag[1]-fc[i][2]*x[2]*flag[2])/fc[i][i];
            flag[i]=1;
        }
        for(i=0; i<3; i++) {
            printf("x[%d]=%lf,y[%d]=%lf\n",i,x[i],i,y[i]);
        }
        for(i=0; i<3; i++) {
            if(temp<fabs(y[i]-x[i])) {
                temp=fabs(y[i]-x[i]);
            }
        }
    } while(temp>eps);

return 0;
}