simp2.txt

来自「VB常用数值算法集 内含有解线性代数方程组」· 文本 代码 · 共 29 行

Sub SIMP2(A(), M, N, MP, NP, L2(), NL2, IP, KP, Q1)
    EPS = 0.000001
    IP = 0
    FLAG = 0
    For I = 1 To NL2
        If A(L2(I) + 1, KP + 1) < -EPS Then FLAG = 1
        If FLAG = 1 Then Exit For
    Next I
    If FLAG = 0 Then Exit Sub
    Q1 = -A(L2(I) + 1, 1) / A(L2(I) + 1, KP + 1)
    IP = L2(I)
    For I = I + 1 To NL2
        II = L2(I)
        If A(II + 1, KP + 1) < -EPS Then
            Q = -A(II + 1, 1) / A(II + 1, KP + 1)
            If Q < Q1 Then
                IP = II
                Q1 = Q
            ElseIf Q = Q1 Then
                For K = 1 To N
                    QP = -A(IP + 1, K + 1) / A(IP + 1, KP + 1)
                    Q0 = -A(II + 1, K + 1) / A(II + 1, KP + 1)
                    If Q0 <> QP Then Exit For
                Next K
                If Q0 < QP Then IP = II
            End If
        End If
    Next I
End Sub