mssgjmainmodule.bas

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

Attribute VB_Name = "MainModule"
Sub Main()
    Dim mtxA(4, 4) As Double
    Dim mtxB(4, 4) As Double
    Dim mtxC(4, 4) As Double
    
    ' 原矩阵
    mtxA(1, 1) = 5:    mtxA(1, 2) = 7:    mtxA(1, 3) = 6:    mtxA(1, 4) = 5
    mtxA(2, 1) = 7:    mtxA(2, 2) = 10:    mtxA(2, 3) = 8:    mtxA(2, 4) = 7
    mtxA(3, 1) = 6:    mtxA(3, 2) = 8:    mtxA(3, 3) = 10:    mtxA(3, 4) = 9
    mtxA(4, 1) = 5:    mtxA(4, 2) = 7:    mtxA(4, 3) = 9:    mtxA(4, 4) = 10
    
    ' 备份原矩阵
    For i = 1 To 4
        For j = 1 To 4
            mtxB(i, j) = mtxA(i, j)
        Next j
    Next i

    '求解
    If MRinv(4, mtxA) Then
        
        ' 矩阵乘法AA-
        Call MMul(4, 4, 4, mtxB, mtxA, mtxC)
        
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & _
                "矩阵A为" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 4, mtxB, "######0.0000000") & Chr$(13) & Chr$(13) & _
                "逆矩阵A-为" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 4, mtxA, "######0.0000000") & Chr$(13) & Chr$(13) & _
                "矩阵A与其逆矩阵A-之乘积矩阵为" & Chr$(13) & Chr$(13) & _
                MatrixToString(4, 4, mtxC, "######0.0000000")
    Else
        MsgBox "求解失败！"
    End If
    
End Sub