tridag.txt

来自「用VB实现线性方程组的求解。包括：高斯消元法」· 文本 代码 · 共 16 行

Sub TRIDAG(A(), B(), C(), R(), U(), N)
    NMAX = 100
    Dim GAM(100)
    If B(1) = 0# Then Exit Sub
    BET = B(1)
    U(1) = R(1) / BET
    For J = 2 To N
        GAM(J) = C(J - 1) / BET
        BET = B(J) - A(J) * GAM(J)
        If BET = 0# Then Exit Sub
        U(J) = (R(J) - A(J) * U(J - 1)) / BET
    Next J
    For J = N - 1 To 1 Step -1
        U(J) = U(J) - GAM(J + 1) * U(J + 1)
    Next J
End Sub