lemivmainmodule.bas

来自「求解线性方程组」· 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