C语言 用矩阵计算抛物线方程

转载至C语言经典算法之高斯消元法_c语言高斯消元法化简矩阵-CSDN博客

这个文章的C语言代码我在使用的时候，出了点BUG，而且最后的值不对，所以略做更改

void gaussElimination(void) {
    for (int k = 0; k < N; ++k) {
        // 寻找主元，并交换行（如果需要）
        int maxIndex = k;
        float maxValue = fabs(A[k][k]);
        for (int i = k + 1; i < N; ++i) {
            if (fabs(A[i][k]) > maxValue) {
                maxIndex = i;
                maxValue = fabs(A[i][k]);
            }
        }
        if (maxIndex != k) {
            for (int j = k; j <= N; ++j) {
                float temp = A[k][j];
                A[k][j] = A[maxIndex][j];
                A[maxIndex][j] = temp;
            }
        }
 
				//归一化
				float ftemp =  A[k][k];
				for (int i = k; i <= N; ++i) {
						A[k][i] /= ftemp;
				}
        // 利用主元消元
        for (int i = k + 1; i < N; ++i) {
						if(fabs(A[i][k]) > 0.000001f){
							for (int j = k + 1; j <= N; ++j) {
									A[i][j] -= ( A[k][j] * A[i][k]);
							}
						}
            A[i][k] = 0; // 确保消元后主元下方元素为0
        }
    }
}

void backSubstitution(void) {
    for (int i = N - 1; i > 0; --i) {
        for (int j = i - 1; j >= 0; --j) {
            A[j][N] -= (A[i][N] * A[j][i]);
        }
    }
}