Как вычислить определитель матрицы для NxN? ⇐ C#

[Anonymous]

Мне нужно найти определитель матрицы NxN, который должен работать как минимум для 5x5. У меня есть функция, переведенная из Паскаля, но есть ИСКЛЮЧЕНИЕ INDEX OUT OF RANGE. Может ли кто-нибудь мне помочь?
Вот мой код:
public static double DET(double[,] a, int n)
{
int i, j, k;
double det = 0;
for (i = 0; i < n - 1; i++)
{
for (j = i + 1; j < n + 1; j++)
{
det = a[j, i] / a[i, i];
for (k = i; k < n; k++)
a[j, k] = a[j, k] - det * a[i, k]; // Here's exception
}
}
det = 1;
for (i = 0; i < n; i++)
det = det * a[i, i];
return det;
}


Подробнее здесь: https://stackoverflow.com/questions/505 ... nt-for-nxn