Из-за деления на нуль программа выдает -nan в случае, если первый элемент матрицы нулевой. Как можно было бы исправить эту ошибку ?
#include <stdio.h>
int main () {
int n, i, j, k;
scanf("%d", &n);
double a[n][n], t, det;
for(i = 0; i < n; i++) {
for(j = 0; j < n; j++) {
scanf("%lf", &a[i][j]);
}
}
for(k = 0; k < n; k++) {
for(i = k+1; i < n; i++) {
t = a[i][k] / a[k][k];
a[i][k] = 0;
for(j = k+1; j < n; j++) {
a[i][j] = a[i][j] - t * a[k][j];
}
}
}
det = 1.0;
for(k = 0; k < n; k++) {
det = det * a[k][k];
}
printf ("det: %lf\n", det);
}