我c语言学的不好，然后都不大会调试，这个作业马上要交了，我已经自己很努力了，但是依然不能得到结果
方程组是：
8.6x+1.5+3.4z=1.9;
2.1x+9.8y+3.4z=6.7;
2.7x+1.8y+7.2z=2.4;
精确解是x=0.0498801697;y=0.617376037;z=0.160284261
我写的代码如下：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
double dmax(double a,double b,double c)
    {
        double m,max;
        m=a>b?a:b;
        max=m>c?m:c;
        return max;
    }
void main()
{
    double a[3][3]={8.6,1.5,3.4,2.1,9.8,3.4,2.7,1.8,7.2};
    double b[3]={1.9,6.7,2.4};
    double x[3]={0},d[3]={0};
    double temp;
    double M=0;
    int i,j,k,h;
   
    for(k=0;M<0.00005;k++)
    {
        for(i=0;i<3;i++)
        {
            for(j=0;j<3;j++)
            {
                if(j!=i)
                    temp+=a[i][j]*d[j];
            }
            x[i]=(b[i]-temp)/a[i][i];
        }
        for(h=0;h<3;h++)
        {
            b[h]=x[h];
        }
        M=dmax(fabs(x[0]-0.0498801697),fabs(x[1]-0.617376037),fabs(x[2]-0.160284261));
        printf("%lf",M);
    }
    printf("迭代次数：  %d\n",k);
    printf("x,y,z的值分别是%.5lf   %.5lf   %.5lf",x[0],x[1],x[2]);
}
求帮助 谢谢

用线性代数克莱姆法则会好做一些

/*我c语言学的不好，然后都不大会调试，这个作业马上要交了，我已经自己很努力了，但是依然不能得到结果
方程组是：
8.6x+1.5y+3.4z=1.9;
2.1x+9.8y+3.4z=6.7;
2.7x+1.8y+7.2z=2.4;
精确解是x=0.0498801697;y=0.617376037;z=0.160284261
我写的代码如下：*/

#include <stdio.h>
#define E(a,b,c,d,e,f,g,h,i) a*e*i+b*f*g+c*d*h-c*e*g-b*d*i-a*f*h;
int main()
{
    double x,y,z,d,d1,d2,d3;
    d = E(8.6,1.5,3.4,2.1,9.8,3.4,2.7,1.8,7.2);
    d1 = E(1.9,1.5,3.4,6.7,9.8,3.4,2.4,1.8,7.2);
    d2 = E(8.6,1.9,3.4,2.1,6.7,3.4,2.7,2.4,7.2);
    d3 = E(8.6,1.5,1.9,2.1,9.8,6.7,2.7,1.8,2.4);
    x=d1/d;
    y=d2/d;
    z=d3/d;
    printf("x=%.10f\ny=%.10f\nz=%.10f",x,y,z);
}

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

int main()
{
   
    double a[3][3]={8.6,1.5,3.4,2.1,9.8,3.4,2.7,1.8,7.2};
    double b[3]={1.9,6.7,2.4};
    double x[3]={0},d[3]={0};
    double m=1,temp;
    int i,j;
    while(m>0.00005)
    {
        for(i=0;i<3;++i)
            d[i]=x[i];
        for(i=0;i<3;++i)
        {
            x[i]=0;
            for(j=0;j<3;++j)
            {
                if(j!=i)
                {
                    x[i]+=a[i][j]*d[j];
                }
            }
            x[i]=(b[i]-x[i])/a[i][i];
        }
        m=fabs(x[0]-d[0]);
        for(i=1;i<3;++i)
        {
            temp=fabs(x[i]-d[i]);
            if(temp>m)
                m=temp;
        }
    }
    for(i=0;i<3;++i)
        printf("%.5lf\n",x[i]);
    getch();
    return 0;
}

论坛真是高人多啊，第二楼的真的好厉害 但是我看不懂 呵呵 是c吗？
还有 你们都没给出求迭代次数的k值