lecpxgaussmainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblAR(4, 4) As Double, dblAI(4, 4) As Double
    Dim dblBR(4) As Double, dblBI(4) As Double
    Dim sCpx1 As String, sCpx2 As String, sCpx3 As String, sCpx4 As String
    
    '系数矩阵
    dblAR(1, 1) = 1
    dblAR(1, 2) = 3
    dblAR(1, 3) = 2
    dblAR(1, 4) = 13
    dblAR(2, 1) = 7
    dblAR(2, 2) = 2
    dblAR(2, 3) = 1
    dblAR(2, 4) = -2
    dblAR(3, 1) = 9
    dblAR(3, 2) = 15
    dblAR(3, 3) = 3
    dblAR(3, 4) = -2
    dblAR(4, 1) = -2
    dblAR(4, 2) = -2
    dblAR(4, 3) = 11
    dblAR(4, 4) = 5
    
    dblAI(1, 1) = 3
    dblAI(1, 2) = -2
    dblAI(1, 3) = 1
    dblAI(1, 4) = 6
    dblAI(2, 1) = -2
    dblAI(2, 2) = 7
    dblAI(2, 3) = 5
    dblAI(2, 4) = 8
    dblAI(3, 1) = 9
    dblAI(3, 2) = -3
    dblAI(3, 3) = 15
    dblAI(3, 4) = 1
    dblAI(4, 1) = -2
    dblAI(4, 2) = -2
    dblAI(4, 3) = 7
    dblAI(4, 4) = 6
    '常数向量
    dblBR(1) = 2
    dblBR(2) = 7
    dblBR(3) = 3
    dblBR(4) = 9
    
    dblBI(1) = 1
    dblBI(2) = 2
    dblBI(3) = -2
    dblBI(4) = 3
    
    '求解
    If LECpxGauss(4, dblAR, dblAI, dblBR, dblBI) = True Then
        sCpx1 = dblBR(1) & " + " & dblBI(1) & "j"
        sCpx2 = dblBR(2) & " + " & dblBI(2) & "j"
        sCpx3 = dblBR(3) & " + " & dblBI(3) & "j"
        sCpx4 = dblBR(4) & " + " & dblBI(4) & "j"
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & "x1 = " & sCpx1 & Chr$(13) & "x2 = " & sCpx2 & Chr$(13) & "x3 = " & sCpx3 & Chr$(13) & "x4 = " & sCpx4
    Else
        MsgBox "求解失败！"
    End If
    
End Sub