lecpxgaussjordenmainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblAR(4, 4) As Double, dblAI(4, 4) As Double
    Dim dblBR(4, 2) As Double, dblBI(4, 2) As Double
    Dim sCpx11 As String, sCpx12 As String, sCpx21 As String, sCpx22 As String
    Dim sCpx31 As String, sCpx32 As String, sCpx41 As String, sCpx42 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, 1) = 2
    dblBR(1, 2) = -2
    dblBR(2, 1) = 7
    dblBR(2, 2) = 3
    dblBR(3, 1) = 3
    dblBR(3, 2) = 2
    dblBR(4, 1) = 9
    dblBR(4, 2) = 1
    
    dblBI(1, 1) = 1
    dblBI(1, 2) = 3
    dblBI(2, 1) = 2
    dblBI(2, 2) = 7
    dblBI(3, 1) = -2
    dblBI(3, 2) = 9
    dblBI(4, 1) = 3
    dblBI(4, 2) = 2
    
    '求解
    If LECpxGaussJordan(4, 2, dblAR, dblAI, dblBR, dblBI) = True Then
        sCpx11 = Round(dblBR(1, 1), 5) & " + " & Round(dblBI(1, 1), 5) & "j"
        sCpx12 = Round(dblBR(1, 2), 5) & " + " & Round(dblBI(1, 2), 5) & "j"
        sCpx21 = Round(dblBR(2, 1), 5) & " + " & Round(dblBI(2, 1), 5) & "j"
        sCpx22 = Round(dblBR(2, 2), 5) & " + " & Round(dblBI(2, 2), 5) & "j"
        sCpx31 = Round(dblBR(3, 1), 5) & " + " & Round(dblBI(3, 1), 5) & "j"
        sCpx32 = Round(dblBR(3, 2), 5) & " + " & Round(dblBI(3, 2), 5) & "j"
        sCpx41 = Round(dblBR(4, 1), 5) & " + " & Round(dblBI(4, 1), 5) & "j"
        sCpx42 = Round(dblBR(4, 2), 5) & " + " & Round(dblBI(4, 2), 5) & "j"
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
                "x11 = " & sCpx11 & Chr$(13) & _
                "x21 = " & sCpx21 & Chr$(13) & _
                "x31 = " & sCpx31 & Chr$(13) & _
                "x41 = " & sCpx41 & Chr$(13) & Chr$(13) & _
                " x12 = " & sCpx12 & Chr$(13) & _
                " x22 = " & sCpx22 & Chr$(13) & _
                " x32 = " & sCpx32 & Chr$(13) & _
                " x42 = " & sCpx42
    Else
        MsgBox "求解失败！"
    End If
    
End Sub