使用全选主元高斯消去法求解线性方程组程序.doc

使用全选主元高斯消去法求解线性方程组程序 参数 n - Integer型变量，线性方程组的阶数 dblA - Double型 n x n 二维数组，线性方程组的系数矩阵 dblB - Double型长度为 n 的一维数组，线性方程组的常数向量，返回方程组的解向量 返回值：Boolean型，求解成功为True，无解或求解失败为False====函数代码如下====Public Function LEGauss(n As Integer, dblA() As Double, dblB() As Double) As Boolean 局部变量 Dim i As Integer, j As Integer, k As Integer Dim nIs As Integer ReDim nJs(n) As Integer Dim d As Double, t As Double  开始求解 For k = 1 To n - 1 d = 0# 归一 For i = k To n For j = k To n t = Abs(dblA(i, j)) If t d Then d = t nJs(k) = j nIs = i End If Next j Next i 无解，返回 If d + 1# = 1# Then LEGauss = False Exit Function End If 消元 If nJs(k) k Then For i = 1 To n t = dblA(i, k) dblA(i, k) = dblA(i, nJs(k)) dblA(i, nJs(k)) = t Next i End If  If nIs k Then For j = k To n t = dblA(k, j) dblA(k, j) = dblA(nIs, j) dblA(nIs, j) = t Next j t = dblB(k) dblB(k) = dblB(nIs) dblB(nIs) = t End If  d = dblA(k, k) For j = k + 1 To n dblA(k, j) = dblA(k, j) / d Next j  dblB(k) = dblB(k) / d For i = k + 1 To n For j = k + 1 To n dblA(i, j) = dblA(i, j) - dblA(i, k) * dblA(k, j) Next j dblB(i) = dblB(i) - dblA(i, k) * dblB(k) Next i Next k  d = dblA(n, n)  无解，返回 If Abs(d) + 1# = 1# Then LEGauss = False Exit Function End If  回代 dblB(n) = dblB(n) / d For i = n - 1 To 1 Step -1 t = 0# For j = i + 1 To n t =