lemivmainmodule.bas

来自「科学与工程数值计算算法（Visual Basic版）一书的源代码。对科学研究者自」· BAS 代码 · 共 35 行

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblA(4, 3) As Double
    Dim dblAP(3, 4) As Double
    Dim dblU(4, 4) As Double
    Dim dblV(3, 3) As Double
    Dim dblB(4) As Double
    Dim dblX(3) As Double
    
    '系数矩阵
    dblA(1, 1) = 1:   dblA(1, 2) = 1:    dblA(1, 3) = -1
    dblA(2, 1) = 2:   dblA(2, 2) = 1:   dblA(2, 3) = 0
    dblA(3, 1) = 1:    dblA(3, 2) = -1:   dblA(3, 3) = 0
    dblA(4, 1) = -1:    dblA(4, 2) = 2:    dblA(4, 3) = 1

    '常数向量
    dblB(1) = 2
    dblB(2) = -3
    dblB(3) = 1
    dblB(4) = 4
    
    '求解
    If LEMiv(4, 3, dblA, dblB, dblX, dblAP, dblU, dblV, 5, 0.000001) Then
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
                "x1 = " & Round(dblX(1), 6) & Chr$(13) & _
                "x2 = " & Round(dblX(2), 6) & Chr$(13) & _
                "x3 = " & Round(dblX(3), 6) & Chr$(13) & Chr$(13) & _
                "矩阵A+" & Chr$(13) & Chr$(13) & _
                MatrixToString(3, 4, dblAP, "######0.0000000")
    Else
        MsgBox "求解失败！"
    End If
    
End Sub