Обратная матрица.Нужно создать функцию для обратной матрицы . Подскажите пожалуйста,как сделать

(P.s Необходимо работать только со статическими массивами).Я пытался сделать обратную матрицу через единичную,но не получается)

double dopolnit(double a[n][n],int nl1,int nc1)
{
cout<<'\n';
int i,j,k,N;
double vspom[n];
if (nl1!=nc1)
printf("Нарушено условие работы");
else
{
N=nl1=nc1;

for(k=0;k<N-1;k++)
{
for(i=1+k;i<N;i++)
{
for(j=0;j<N;j++)
{
vspom[j] = [/CPP]a[k][j]*a[i][k]/a[k][k];
//printf(" %d",vspom[j]);
}
for(j=0;j<N;j++)
{
a[i][j] = a[i][j]-vspom[j];

}
}
}
return a[n][n];
}
}
double obratn(double a[n][n],int nl1,int nc1)
{
int i,j,k,nc;
double e[n][n];
double pom[n];
double q[n][n];
double sum;
for (int i = 0; i < nl1; i++)
{
for (int j = 0; j < nl1; j++)
{
if (i == j) // Если диагональ
e[i][j] = 1;
else // а если не
e[i][j] = 0;
cout << e[i][j] << " ";
}
cout << endl;
}
cout<<'\n';
for(k=0;k<nl1-1;k++)
{
for(i=1+k;i<nl1;i++)
{
for(j=0;j<nl1;j++)
{
pom[j] = e[k][j]*a[i][k]/a[k][k];
e[i][j] = e[i][j]-pom[j];
}
for(i=0;i<nl1;i++)
{
for(j=0;j<nl1;j++)
{
cout << e[i][j] << " ";
}
cout << endl; //
}
}
}
dopolnit(a,nl1,nc1);
for(i=0;i<nl1;i++)
{
for(j=nl1;j>=0;j--)
{
sum=0;
for(k=nl1;k>=0;k--)
{
sum+=a[j][k]*e[k][i];
q[j][i]=(e[j][i]-sum)/a[j][i];
cout<<q[j][i];
}
putchar('\n');
}
}
return q[n][n];
}

пытался сделать,как здесь http://algorithmlib.org/gauss_inverse

0