数值分析里面用高斯列主元方法解方程的，后8个解求得正确，最后一个解出错（也就是回代的最后一步出错了），我实在找不出原因，大家帮帮忙啊

#include<iostream.h>
#include<math.h>
#include<iomanip.h>
void main(void)
{
    double a[9][10]={{31,-13,0,0,0,-10,0,0,0,-15},{-13,35,-9,0,-11,0,0,0,0,27},{0,-9,
        31,-10,0,0,0,0,0,-23},{0,0,-10,79,-30,0,0,0,-9,0},{0,0,0,-30,57,-7,0,-5,0,-20},
    {0,0,0,0,-7,47,-30,0,0,12},{0,0,0,0,0,-30,41,0,0,-7},{0,0,0,0,-5,0,0,27,-2,7},{
        0,0,0,-9,0,0,0,-2,29,10}};
    double x[8]={0,0,0};
    double z=0;//交换两行时用的中间值//
    int k=0;//列主元时记录最大主元的所在行//
    for(int i=0;i<=7;i++)//找最大主元//
    {
        for(int j=1;j+i<=8;j++)
        {
            if(fabs(a[i+j][i])>fabs(a[i][i]))
                k=i+j;
            else
                k=i;
        }
        for(int v=i;v<=9;v++)//交换两行//
        {
            z=a[i][v];
            a[i][v]=a[k][v];
            a[k][v]=z;
        }
        for(int w=i+1;w<=8;w++)//列主元后消元，w为计算时的行空竹量//
        {
            double H=-a[w][i]/a[i][i];
            for(int x=0;i+x<=9;x++)//i+x为计算的列控制量//
            {
                a[w][i+x]=a[w][i+x]+a[i][i+x]*H;
            }
        }
    }
        for(int m=0;m<=8;m++)//输出处理后的矩阵//
        {
            for(int n=0;n<=9;n++)
            {
                if(n!=0&&n%9==0)
                cout<<setw(14)<<a[m][n]<<endl;
            else
                cout<<setw(14)<<a[m][n];
            
            }
        }
        cout<<a[0][0]<<endl;//检测a[0][0]//
        x[8]=a[8][9]/a[8][8];
        for(int n=7;n>=0;--n)//回代过程开始//
        {
            double g=0;
            for(int m=n+1;m<=8;m++)
            {
                g+=a[n][m]*x[m];
            }
            x[n]=(a[n][9]-g)/a[n][n];
            cout<<setw(18)<<a[n][n]<<setw(18)<<g<<setw(18)<<x[n]<<endl;
        }//回代结束//
        cout<<a[0][0]<<endl;//检测a[0]//
        for(int t=0;t<=8;t++)
            cout<<setw(14)<<x[t];
        cout<<endl;
}
运行结果那个-30.8937是错的（应该是-0.28923），后面几个是对的。