C言語 N次正方行列の行列式を求めたいです

Question

前提・実現したいこと

C言語　N次正方行列の行列式を求めたいですがすべて０が返されてしまいます
コードの訂正をお願いしたいです

発生している問題・エラーメッセージ

すべて０が返されてしまう

該当のソースコード

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

int main(void) {
    int i,j,k,N;
    printf("Nを入力してください");
    scanf("%d",&N);
    double M[N][N],det=1.0,buf;
    for(i=1;i<=N;i++){
            for(j=1;j<=N;j++){
            M[i-1][j-1]=sqrt(i+j)/(i*j);
            
        }
    }
     
    //三角行列を作成
    for(i=0;i<N;i++){
     for(j=0;j<N;j++){
      if(i<j){
       buf=M[j][i]/M[i][i];
       for(k=0;k<N;k++){
       M[j][k]-=M[i][k]*buf;
       }
      }
     }
    }
     
    //対角部分の積
    for(i=0;i<N;i++){
     det*=M[i][i];
    }
     
    printf("%f\n",det);
    return 0;
 }

。

補足情報（FW/ツールのバージョンなど）

Answer 1

思考停止せずにデバッグして下さい。
途中経過を確認すると一発で分かりますよ。

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

int main(void) {
    int i,j,k,N;
    printf("Nを入力してください");
    scanf("%d",&N);
    printf("\n");
    double M[N][N],det=1.0,buf;
    for(i=1;i<=N;i++){
     for(j=1;j<=N;j++){
      M[i-1][j-1]=sqrt(i+j)/(i*j);
     }
    }

    //三角行列を作成
    for(i=0;i<N;i++){
     for(j=0;j<N;j++){
      if(i<j){
       buf=M[j][i]/M[i][i];
       for(k=0;k<N;k++){
        M[j][k]-=M[i][k]*buf;
       }
      }
     }
    }

    //対角部分の積
    for(i=0;i<N;i++){
     det*=M[i][i];
     printf("%d %lf\n", i, M[i][i]);
    }

    printf("%f\n",det);
    return 0;
 }

出力結果
Nを入力してください10
0 1.414214
1 -0.030330
2 -0.000388
3 -0.000011
4 -0.000000
5 -0.000000
6 -0.000000
7 -0.000000
8 -0.000000
9 -0.000000
-0.000000

Comments

[Hayato0511]

思考停止せずにデバックで自己解決できるよう努めていきたいと思います。
ありがとうございます。

[kaina]

頑張って下さい！
尚、細かい話ではありますが「デバック」では無く、「デバッグ」です。

[Hayato0511]

ご指摘ありがとうございます。

Answer 2

デバッガが使えるんですから、途中で実行を止めて変数の値を見ていきましょう。
そうすればどこでおかしくなってるのかわかるでしょう

Comments

[Hayato0511]

デバッガの存在を教えてくださりありがとうございます。
初心者脱却できそうです。

Answer 3

通常にコンパイルして実行するとゼロではない結果が返ってきました

$ gcc -lm -o det det.c
$ ./det
Nを入力してください1
1.414214
$ ./det
Nを入力してください2
-0.042893
$ ./det
Nを入力してください3
0.000017

Comments

[Hayato0511]

N=6,7を代入していたため、０が返されていたようです。
ご指摘ありがとうございます。