mrinvmainmodule.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) = 0.2368:    mtxA(1, 2) = 0.2471:    mtxA(1, 3) = 0.2568:    mtxA(1, 4) = 1.2671
    mtxA(2, 1) = 1.1161:    mtxA(2, 2) = 0.1254:    mtxA(2, 3) = 0.1397:    mtxA(2, 4) = 0.149
    mtxA(3, 1) = 0.1582:    mtxA(3, 2) = 1.1675:    mtxA(3, 3) = 0.1768:    mtxA(3, 4) = 0.1871
    mtxA(4, 1) = 0.1968:    mtxA(4, 2) = 0.2071:    mtxA(4, 3) = 1.2168:    mtxA(4, 4) = 0.2271
    
    ' 备份原矩阵
    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