leseidelmainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblA(4, 4) As Double
    Dim dblB(4) As Double
    Dim x(4) As Double
    
    '系数矩阵
    '系数矩阵
    dblA(1, 1) = 7:   dblA(1, 2) = 2:    dblA(1, 3) = 1:    dblA(1, 4) = -2
    dblA(2, 1) = 9:   dblA(2, 2) = 15:   dblA(2, 3) = 3:    dblA(2, 4) = -2
    dblA(3, 1) = -2:   dblA(3, 2) = -2:    dblA(3, 3) = 11:   dblA(3, 4) = 5
    dblA(4, 1) = 1:   dblA(4, 2) = 3:    dblA(4, 3) = 2:    dblA(4, 4) = 13

    '常数向量
    dblB(1) = 4
    dblB(2) = 7
    dblB(3) = -1
    dblB(4) = 0
    
    '求解
    If LESeidel(4, dblA, dblB, x, 0.000001) = True Then
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
                "x1 = " & Round(x(1), 7) & Chr$(13) & _
                "x2 = " & Round(x(2), 7) & Chr$(13) & _
                "x3 = " & Round(x(3), 7) & Chr$(13) & _
                "x4 = " & Round(x(4), 7)
    Else
        MsgBox "求解失败！"
    End If
    
End Sub