Алгоритм нахождения обратной мтарицы

[AlinAA]

проверьте,пожалуйста,правильный ли... самой нет пока возможности проверить..задействован метод гаусса

#include <stdio.h>
#include <math.h>
#include <stdlib.h>

/* получение обратной матрицы */
extern int invert(double* matrix, int* nstr, int* nstb)
{ int N = *nstr, i, j, k;
double **e;
double *tmp;

if (*nstr != *nstb)
{
puts("matrica ne kvadratnaia");
return 0;
}



/* Выделяем память под единичную матрицу */
e= (double**)malloc(N*sizeof(double));
for (k = 0; k < N; k++)
e[k] = (double*)malloc(N*sizeof(double));
/* ------------------------ */

/* Создадим единичную матрицу */
for (i=0;i<N;i++)
for (j=0;j<N;j++)
e[i][j]=(i==j?1:0);

/* Действия с обратной матрицей */
for(i=0;i<N;i++)
{

tmp = (double*)malloc(sizeof(double));
*tmp=*(matrix+i*N+i);
for(j=N-1;j>=0;j--)
{
e[i][j]/=*tmp;
*(matrix+i*N+j)/=*tmp;
}

/* excluding i-th element from each row except i-th one */
for(j=0;j<N;j++)
if (j!=i)
{
*tmp =*(matrix + j*N + i);

for(k=N-1;k>=0;k--)
{
e[j][k]-=e[i][k] * (*tmp);
*(matrix+j*N+k) -= *(matrix+i*N+k) * (*tmp);
}
}
free((void*)tmp);
}

/* Копируем единичную матрицу к исходной */
for(i=0;i<N;i++)
for(j=0;j<N;j++)
*(matrix+i*N+j)=e[i][j];

/* освобождаем память */
for (i = 0; i != N; ++i)
free((void*)e[i]);
free((void**)e);

return 1;
}

/*получение нормы матрицы */
extern double norma(double* matrix, int* nstr, int* nstb)
{
int N = *nstr, i, j; double max_stb = *(matrix); double sum_ = 0.0;

for(i=0;i<N;i++)
{
for(j=0;j<N;j++) /* считаем сумму в столбце в цикле */
{
sum_ += *(matrix+j*N + i);
}
if (sum_ > max_stb) /* определяем максимальный элемент в столбце */
max_stb = sum_;
}
return fabs(max_stb); /* возвращяем модуль максимальной суммы в столбцах */
}

[como]

Для начала, вместо:
e= (double**)malloc(N*sizeof(double));
должно быть:
e = (double**)malloc(N*sizeof(double *));