msymtrimainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim mtxA(5, 5) As Double
    Dim mtxQ(5, 5) As Double
    Dim mtxT(5, 5) As Double
    Dim dblB(5) As Double, dblC(5) As Double
    
    ' 原矩阵
    mtxA(1, 1) = 10:  mtxA(1, 2) = 1: mtxA(1, 3) = 2:  mtxA(1, 4) = 3: mtxA(1, 5) = 4
    mtxA(2, 1) = 1:  mtxA(2, 2) = 9: mtxA(2, 3) = -1:  mtxA(2, 4) = 2: mtxA(2, 5) = -3
    mtxA(3, 1) = 2:  mtxA(3, 2) = -1: mtxA(3, 3) = 7:  mtxA(3, 4) = 3: mtxA(3, 5) = -5
    mtxA(4, 1) = 3:  mtxA(4, 2) = 2: mtxA(4, 3) = 3:  mtxA(4, 4) = 12: mtxA(4, 5) = -1
    mtxA(5, 1) = 4:  mtxA(5, 2) = -3: mtxA(5, 3) = -5:  mtxA(5, 4) = -1: mtxA(5, 5) = 15

    '求解
    Call MSymTri(5, mtxA, mtxQ, mtxT, dblB, dblC)
    
    MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
            "原矩阵A" & Chr$(13) & Chr$(13) & _
            MatrixToString(5, 5, mtxA, "######0.0000000") & Chr$(13) & Chr$(13) & _
            "矩阵Q" & Chr$(13) & Chr$(13) & _
            MatrixToString(5, 5, mtxQ, "######0.0000000") & Chr$(13) & Chr$(13) & _
            "对称三对角矩阵T" & Chr$(13) & Chr$(13) & _
            MatrixToString(5, 5, mtxT, "######0.0000000")

End Sub