lecholeskymainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblA(4, 4) As Double
    Dim dblB(4, 2) As Double
    
    '系数矩阵
    dblA(1, 1) = 5:   dblA(1, 2) = 7:    dblA(1, 3) = 6:    dblA(1, 4) = 5
    dblA(2, 1) = 7:   dblA(2, 2) = 10:   dblA(2, 3) = 8:    dblA(2, 4) = 7
    dblA(3, 1) = 6:   dblA(3, 2) = 8:    dblA(3, 3) = 10:   dblA(3, 4) = 9
    dblA(4, 1) = 5:   dblA(4, 2) = 7:    dblA(4, 3) = 9:    dblA(4, 4) = 10
    '常数矩阵
    dblB(1, 1) = 23:  dblB(2, 1) = 32:   dblB(3, 1) = 33:   dblB(4, 1) = 31
    dblB(1, 2) = 92:  dblB(2, 2) = 128:  dblB(3, 2) = 132:  dblB(4, 2) = 124
    
    '求解
    If LECholesky(4, 2, dblA, dblB) = True Then
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
                "x11 = " & Round(dblB(1, 1), 5) & "       x12 = " & Round(dblB(1, 2), 5) & Chr$(13) & _
                "x21 = " & Round(dblB(2, 1), 5) & "       x22 = " & Round(dblB(2, 2), 5) & Chr$(13) & _
                "x31 = " & Round(dblB(3, 1), 5) & "       x32 = " & Round(dblB(3, 2), 5) & Chr$(13) & _
                "x41 = " & Round(dblB(4, 1), 5) & "       x42 = " & Round(dblB(4, 2), 5)
    Else
        MsgBox "求解失败！"
    End If
    
End Sub