Нахождение обратной матрицы

Sub Test_invmatr()
Dim X(1 To 3, 1 To 3) As Double
Dim R() As Double
 
    X(1, 1) = 1
    X(1, 2) = 2
    X(2, 1) = -2
    X(2, 2) = 1
    X(3, 3) = 1
 
    R = Inv_Matr(X)
 
    If Ubound(R,1)=1 then
       MsgBox "Обратная матрица не существует. Опр=0"
    else
       For i% = 1 To 3
            For j% = 1 To 3
                 Debug.Print R(i%, j%); " ";
            Next j%
            Debug.Print
       Next i%
    End if
End Sub

Function Inv_Matr(matr() As Double) As Double()
Dim A() As Double
Dim R() As Double
 
    n% = UBound(matr, 1)
    
    ReDim A(1 To n%, 1 To 2 * n%) As Double
 
    M% = 2 * n%
    s% = 0
    
    For i% = 1 To n%
        For j% = 1 To M%
            If j% <= n% Then
               A(i%, j%) = matr(i%, j%)
            Else
               If j% = n% + i% Then
                  A(i%, j%) = 1
               Else
                  A(i%, j%) = 0
               End If
            End If
        Next j%
    Next i%
 
    For i% = 1 To n%
    
        k% = i%
        
        s% = 0
        
        Do
        
           If Abs(A(k%, i%)) <= 0.00000000000001 Then
              s% = 1
              If k% < n% Then
                 k% = k% + 1
              Else
                 ReDim R(1 To 1) As Double
                 Inv_Matr = R
                 Exit Function
              End If
           Else
              Exit Do
           End If
                
        Loop
    
        If s% = 1 Then
        
           For j% = 1 To M%
               T# = A(k%, j%)
               A(k%, j%) = A(i%, j%)
               A(i%, j%) = T#
           Next j%
        
        End If
        
        For j% = M% To i% Step -1
            A(i%, j%) = A(i%, j%) / A(i%, i%)
        Next j%
        
        For k% = 1 To n%
            If k% <> i% Then
               For j% = M% To i% Step -1
                   A(k%, j%) = A(k%, j%) - A(i%, j%) * A(k%, i%)
               Next j%
            End If
        Next k%
    
    Next i%
    
    ReDim R(1 To n%, 1 To n%) As Double
    
    For i% = 1 To n%
        For j% = 1 To n%
            R(i%, j%) = A(i%, j% + n%)
        Next j%
    Next i%
    
    Inv_Matr = R
 
End Function
 
'##### Тест - матрица Гильберта ######
 
Function Hilbert(n As Integer) As Double()
Dim X() As Double
 
    ReDim X(1 To n, 1 To n) As Double
    
    p% = 0
    For i% = 1 To n
        p% = p% + 1
        q% = p%
        For j% = 1 To n
            X(i%, j%) = 1 / q%
            q% = q% + 1
        Next j%
    Next i%
 
    Hilbert = X
 
End Function
 
 
Sub Test_invmatr()
Dim R() As Double
Dim X() As Double
 
    n% = 10
 
    X = Hilbert(n%)
 
    R = Inv_Matr(X)
    
    For i% = 1 To n%
        For j% = 1 To n%
            Debug.Print R(i%, j%); " ";
        Next j%
        Debug.Print
    Next i%
 
End Sub