ファイルから読み取った5×5行列の行列式を計算したい

前提・実現したいこと

C言語でmatrix,datというテキストファイルで与えられた5×5行列の行列式を計算しようとしています。
ファイルから行列を読みとり配列に格納しましたが、
その後三角行列を作成し対角部分の積を計算したさいに計算結果が合いません。
どの部分がおかしいのでしょうか？
正しい答えは9262です。

該当のソースコード

C言語

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

int main(void){
int matrix[5][5] ;
int det=1,x; 
int n=5;   
int a,i,j,k; 
FILE *fp;


fp=fopen("matrix.dat","r");

if(fp==NULL){
    puts("matrix.datは開けません");
}
    
for(a=0;a<5;a++){
    fscanf(fp,"%d %d %d %d %d",&matrix[a][0],&matrix[a][1],&matrix[a][2],&matrix[a][3],&matrix[a][4]);
    printf("%d,%d,%d,%d,%d\n",matrix[a][0],matrix[a][1],matrix[a][2],matrix[a][3],matrix[a][4]);
}    

    
for(i=0;i<n;i++){ 
for(j=0;j<n;j++){ 
if(i<j){ 
x=matrix[j][i]/matrix[i][i]; 
for(k=0;k<n;k++){ 
matrix[j][k]=matrix[j][k]-matrix[i][k]*x;
} 
} 
} 
} 


for(i=0;i<n;i++){ 
det=det*matrix[i][i]; 
} 
     printf("%d ",det);
            fclose(fp);
            return 0;
}
-------------------------------

行動規範の内容に同意します

回答2件

ベストアンサー

こちらに似たサンプルがありました。確認したところ int でなく double が使われているようですね。

■ 行列式の値の求め方、逆行列の作り方の C 言語プログラム
https://thira.plavox.info/blog/2008/06/_c.html

c
1#include<stdio.h>
2#include<math.h>
3#include <stdlib.h>
4
5int main(void) {
6    double matrix[5][5];
7    double det = 1, x;
8    int n = 5;
9    int a, i, j, k;
10    FILE * fp;
11
12
13    fp = fopen("matrix.dat", "r");
14
15    if (fp == NULL) {
16        puts("matrix.datは開けません");
17    }
18
19    for (a = 0; a < 5; a++) {
20        fscanf(fp, "%lf %lf %lf %lf %lf", &matrix[a][0], &matrix[a][1], &matrix[a][2], &matrix[a][3], &matrix[a][4]);
21        printf("%f,%f,%f,%f,%f\n", matrix[a][0], matrix[a][1], matrix[a][2], matrix[a][3], matrix[a][4]);
22    }
23
24
25    for (i = 0; i < n; i++) {
26        for (j = 0; j < n; j++) {
27            if (i < j) {
28                x = matrix[j][i] / matrix[i][i];
29                for (k = 0; k < n; k++) {
30                    matrix[j][k] = matrix[j][k] - matrix[i][k] * x;
31                }
32            }
33        }
34    }
35
36
37    for (i = 0; i < n; i++) {
38        det = det * matrix[i][i];
39    }
40    printf("%lf ", det);
41    fclose(fp);
42    return 0;
43}

以下、実行結果です。

5.000000,0.000000,2.000000,8.000000,-1.000000
-3.000000,9.000000,10.000000,0.000000,0.000000
-7.000000,1.000000,1.000000,0.000000,2.000000
0.000000,9.000000,7.000000,-1.000000,5.000000
2.000000,-2.000000,3.000000,8.000000,0.000000
9262.000000

投稿2021/08/21 04:40