在C ++中执行朴素高斯消元时,我可以得到一个上三角矩阵,但当方程组的系数大于5 * 5时,它不起作用。知道为什么会这样吗?我使用double作为系数的数据类型。
double factor =0;
for(int k=0;k<(count-1);k++){
for(int i=(k+1);i<count;i++){
factor = (doubleCoefficient[i][k])/(doubleCoefficient[k][k]);
for(int j=0;j<count;j++){
doubleCoefficient[i][j] -= (doubleCoefficient[k][j]*factor);
}
}
}