Declare Function gelimd Lib "ged2.dll" (ByRef av As Double, ByRef b As Double, ByRef x As Double, ByVal N As Long) As Long
..
Function GE(a As Variant, b As Variant, N As Long) As Variant

' Set up 1D, base 0 arrays, aa and ba, from 2D variant arrays, a and b

a = a.Value
b = b.Value
For i = 1 To N
For j = 1 To N
aa((i - 1) * N + j - 1) = a(i, j)
Next j
ba(i - 1) = b(i, 1)
Next i

' Call gelimd function from ge2.dll

retn = gelimd(aa(0), ba(0), x(0), N)

/* ged.cpp -- Gaussian elimination linear equation solvers.
*
*  (C) 2001, 2004, C. Bond. All rights reserved.
*
*  Simple pivoting on zero diagonal element supported.
*   Eliminates unnecessary  zeroing of lower triangle.
*   Does not scale rows to unity pivot value.
*   Swaps b[] as well as a[][], so a pivot ID vector
*   is not required.
*/
#include "ged2.h"
#include <math.h>

#define EPS 1e-10

int gelimd(double* av, double* b, double* x, int n)
{
    double tmp,pvt,*t;
   
    int i,j,k,itmp;
   
   double ** a = new double*[n];
for( i = 0; i < n; i++)
   a[i] = new double[n];

   for (i=0;i<n;i++) {
      for (j=0;j<n;j++) {
         a[i][j] = av[i*n+j];
      }
   }
   

Declare Function vbrmatinv Lib "d:\Users\Doug\Documents\Visual Studio 2010\Projects\Alglib2\Release\Alglib2.dll" _
(ByRef A2 As Double, ByVal M As Long) As Long

Function RMatInvdll(a As Variant) As Variant
    Dim N As Long, M As Long, N2 As Long, Pivots() As Long
    Dim A2() As Double, i As Long, J As Long, K As Long, Rtn As Long

    If TypeName(a) = "Range" Then a = a.Value2
    M = UBound(a)
    N = UBound(a, 2)

    ' Copy 2D base 1 variant array to 1D base 0 double array
    ReDim A2(0 To M * N - 1)
    For i = 1 To M
        For J = 1 To N
            A2((i - 1) * M + J - 1) = a(i, J)
        Next J
    Next i

    Rtn = vbrmatinv(A2(0), M)

    ' Copy 1D base 0 array to 2D base 1 array
    For i = 1 To M
        For J = 1 To N
            a(i, J) = A2((i - 1) * M + J - 1)
        Next J
    Next i
    ' Assign array to function return value
    RMatInvdll = a
End Function

// vbmat.cpp -- Interface for VBA matrix functions.

#include "vbmat.h"

int vbrmatinv(double* av, int n)
{
    int i,j,info;
   matinvreport rep;
    ap::real_2d_array a2;
   
   a2.setcontent(0, n-1, 0, n-1, av);
   
  rmatrixinverse(a2, n, info, rep);

for( i = 0; i < n; i++){
     for (j=0;j<n;j++) {
   av[i*n+j] = a2(i,j);
      }
   }
return 0 ;
}