mlumainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim mtxA(4, 4) As Double
    Dim mtxL(4, 4) As Double
    Dim mtxU(4, 4) As Double
    
    ' 原矩阵
    mtxA(1, 1) = 2:  mtxA(1, 2) = 4:  mtxA(1, 3) = 4:  mtxA(1, 4) = 2
    mtxA(2, 1) = 3:  mtxA(2, 2) = 3:  mtxA(2, 3) = 12:  mtxA(2, 4) = 6
    mtxA(3, 1) = 2:  mtxA(3, 2) = 4:  mtxA(3, 3) = -1:  mtxA(3, 4) = 2
    mtxA(4, 1) = 4:  mtxA(4, 2) = 2:  mtxA(4, 3) = 1:  mtxA(4, 4) = 1
    
    '求解
    If MLU(4, mtxA, mtxL, mtxU) Then
        MsgBox "分解后的下三角矩阵L为" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 4, mtxL, "######0.0000000") & Chr$(13) & Chr$(13) & _
                "分解后的上三角矩阵U为" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 4, mtxU, "######0.0000000")
    Else
        MsgBox "求解失败"
    End If
    
End Sub