lemivmainmodule.bas

科学与工程数值计算算法（Visual Basic版） 附赠的光盘包含了本书中全部的源代码

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