lemivmainmodule.bas

来自「清华大学2002年出版的《科学与工程数值计算算法VB》配套源码」· 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