Тема: Есть вопрос
Хочу уточнить код алгоритма "Вычисление определителя матрицы методом Гаусса", а именно:
for (int j=i+1; j<=n; ++j)
a[i][j] /= a[i][i];
Тогда как в самом начале алгоритма указано:
vector < vector<double> > a (n, vector<double> (n));
.
IMHO вышеприведенный цикл следует заменить на следующий:
for (int j=i+1; j<n; ++j)
a[i][j] /= a[i][i];
т. к. происходит выход за границы массива при j==n.