int main(void) {
double a[3][3];
a[0][0] = 1;
a[0][1] = 1;
a[0][2] = 1;
a[1][0] = 1;
a[1][1] = 2;
a[1][2] = 3;
a[2][0] = 2;
a[2][1] = 4;
a[2][2] = 5;

int i, j, k;
double p, q;
for (k = 0; k < 3; k++) {
p = a[k][k];
a[k][k] = 1.0;
for (j = 0; j < 3; j++) a[k][j] /= p;
for (i = 0; i < 3; i++) {
if (i != k) {
q = a[i][k];
a[i][k] = 0.0;
for (j = 0; j < 3; j++) a[i][j] -= q * a[k][j];
}
}
}

int c,d;
for(c=0;c<=2;c++) {
for(d=0;d<=2;d++) {
printf("inv a[%d][%d] = %d\n",c,d,a[c][d]);
}
}
return 0;
}

エラーは出ないが、正しい値が求まらない