mprove.txt

来自「本书共有数值计算中常用的VisualBasic子过程近200个」· 文本 代码 · 共 15 行

Sub MPROVE(A(), ALUD(), N, INDX(), B(), X())
    NMAX = 100
    Dim R(100)
    For I = 1 To N
        SDP = -B(I)
        For J = 1 To N
            SDP = SDP + A(I, J) * X(J)
        Next J
        R(I) = SDP
    Next I
    Call LUBKSB(ALUD(), N, INDX(), R())
    For I = 1 To N
        X(I) = X(I) - R(I)
    Next I
End Sub