问题六十二：怎么求一元十次方程在区间内的所有不相等的实根（2）——修正“区间端点零值”问题

前续“问题六十二：怎么求一元十次方程在区间内的所有不相等的实根”和“问题五十九：怎么求一元六次方程在区间内的所有不相等的实根（3）——修正一个问题”

不管是求解一元六次方程还是一元十次方程，都存在这样一个区间端值问题：

当区间端值为0时，在求该区间上不等实根总个数时出错。

下方以一元六次方程为例说明问题：

如上图显示，该函数在区间[-2.5，2.5]上有六个不等实根，分别是-1/2, -1/3, -1/4, 1/4, 1/3, 1/2

所以，该函数在区间[-2.5，0]上应该有三个不等实根，分别是-1/2, -1/3, -1/4

但是，程序输出的结果是这样的：

显示不等实根个数还是6个，这明显是错误的。

我们将程序修改如下：（将端点零值用1e-6表示）

bool roots_num_equation_6th(float ee6[7], double a, double b, int &num) {

        double temp1, temp2 , aaa, bbb;

        double ee77[7][7] = {0.0};
        double fun_a[7] = {0.0};
        double fun_b[7] = {0.0};
        int change_num_a = 0;
        int change_num_b = 0;

        if (fabs(a)<1e-6) {
            aaa = 1e-6;
        }
        else {
            aaa = a;
        }
        if (fabs(b)<1e-6) {
            bbb = 1e-6;
        }
        else {
            bbb = b;
        }

        for (int i=0; i<6; i++) {
        //determine f0, f1
            ee77[0][i] = ee6[i];
            ee77[1][i] = (6-i)*ee6[i];
        }
        ee77[0][6] = ee6[6];

        for (int i=2; i<6; i++) {
        //determine f2, f3, f4, f5
            temp1 = ee77[i-2][0] / ee77[i-1][0];
            temp2 = (ee77[i-2][1] - temp1*ee77[i-1][1]) / ee77[i-1][0];
            for (int j=0; j<(6-i); j++) {
                ee77[i][j] = temp1*ee77[i-1][j+2] + temp2*ee77[i-1][j+1] - ee77[i-2][j+2];
                if (fabs(ee77[i][j])<1e-5) {
                    ee77[i][j] = 0.0;
                }
            }
            ee77[i][6-i] = temp2*ee77[i-1][6-i+1] - ee77[i-2][6-i+2];
            if (fabs(ee77[i][6-i])<1e-5) {
                ee77[i][6-i] = 0.0;
            }

            if ((ee77[i][0] == 0.0) && (ee77[i][1] == 0.0) && (ee77[i][2] == 0.0) &&
                (ee77[i][3] == 0.0) && (ee77[i][4] == 0.0) && (ee77[i][5] == 0.0) &&
                (ee77[i][6] == 0.0)) {
                break;
            }
        }
        if (fabs(ee77[5][0])>1e-5) {
        //determine f6
            temp1 = ee77[4][0] / ee77[5][0];
            temp2 = (ee77[4][1] - temp1*ee77[5][1]) / ee77[5][0];
            ee77[6][0] = temp2*ee77[5][1] - ee77[4][2];
            if (fabs(ee77[6][0])<1e-5) {
                ee77[6][0] = 0.0;
            }
        }
        else {
            ee77[6][0] = 0.0;
        }

        for (int j=0; j<7; j++) {

            fun_a[0] = fun_a[0] + ee77[0][j]*pow( aaa, (6-j));
            fun_b[0] = fun_b[0] + ee77[0][j]*pow( bbb, (6-j));

        }
        for (int i=1; i<7; i++) {
            for (int j=0; j<(7-i); j++) {

                fun_a[i] = fun_a[i] + ee77[i][j]*pow( aaa, ((6-i)-j));
                fun_b[i] = fun_b[i] + ee77[i][j]*pow( bbb, ((6-i)-j));

            }
            if ((fun_a[i]*fun_a[i-1]) < 0) {
                change_num_a ++;
            }
            if ((fun_b[i]*fun_b[i-1]) < 0) {
                change_num_b ++;
            }
        }
        num = abs(change_num_a - change_num_b);
        return true;
}

如上修改之后，程序运行结果如下：