VB.net计算二维数组的行列式

‘根据行列式定义计算行列式’

Public Shared Function det(ByVal a(,) As Double)
        Dim b As Double : Dim p(,) As Double
        Dim i, j, k, l As Integer
        Dim c(1), m, n As Integer
        c = size(a) : m = c(0) : n = c(1) : l = m - 1 'l为子行列式的阶数
        ReDim p(l - 1, l - 1)
        If m <> n And m > 1 And n > 1 Then
            System.Windows.Forms.MessageBox.Show("Error of 'det':求解行列式的矩阵行数与列数不一致")
        ElseIf m = 1 Then
            b = a(0, 0)
        Else
            For k = 0 To m - 1
                For i = 0 To l - 1
                    For j = 0 To l - 1
                        If i < k Then
                            p(i, j) = a(i, j + 1) '初始化子行列式的值*/  
                        Else
                            p(i, j) = a(i + 1, j + 1) '初始化子行列式的值*/