lemqrmainmodule.bas

科学与工程数值计算算法（Visual Basic版）一书的源代码。对科学研究者自己开发领域内计算程序极具参考价值。

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblA(4, 3) As Double
    Dim dblQ(4, 4) As Double
    Dim dblB(4) 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 LEMqr(4, 3, dblA, dblB, dblQ) Then
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
                "x1 = " & Round(dblB(1), 6) & Chr$(13) & _
                "x2 = " & Round(dblB(2), 6) & Chr$(13) & _
                "x3 = " & Round(dblB(3), 6) & Chr$(13) & Chr$(13) & _
                "矩阵Q" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 4, dblQ, "######0.0000000") & Chr$(13) & Chr$(13) & _
                "矩阵R" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 3, dblA, "######0.0000000")
    Else
        MsgBox "求解失败！"
    End If
    
End Sub