###前提・実現したいこと
C言語プログラミングの課題で逆行列を求めるプログラムを組んでいます。
コンパイル自体は通ったのですが、実行した際行列の次数を入力すると
永遠と入力できてしまい、C-Zで中断したところセグメントエラー(coreを出力しました)と表示されました。
###発生している問題・エラーメッセージ
^Zセグメントエラー (coreを出力しました)
###該当のソースコード
#include<stdio.h> #include<stdlib.h> int main(void){ int i,j,k,n,N; printf("行列の次数を入力 n: "); scanf("%d",N); double **M; //N*2N行列 M = (double **)malloc(N*sizeof(double *)); M[0] = (double *)malloc(N*(2*N)* sizeof(double)); for(i=1;i<N;i++) M[i]=M[0]+i*(2*N); for(i=0;i<N;i++){ for(j=0;j<N;j++){ //行列の左半分を入力 printf("行列の成分を入力(%d,%d):",i,j); scanf("%lf",M[i][j]); //行列の右半分を入力 if(i==j) M[i][j+N]=1; else M[i][j+N]=0; } } double pivot, mul; // 対角成分が1で正規化された階段行列を作る for (i = 0; i < N; ++i) { // 対角成分の選択、この値で行成分を正規化 pivot = M[i][i]; for (j = 0; j < N * 2; ++j) { M[i][j] = (1 / pivot) * M[i][j]; } // 階段行列を作る為に、現在の行より下の行につ // i列目の成分が0になるような基本変形をする for (k = i + 1; k < N; ++k) { mul = M[k][i]; for (n = i; n < N * 2; ++n) { M[k][n] = M[k][n] - mul * M[i][n]; } } } // 下から上に向かって簡約化 for (i = N - 1; i > 0; --i) { for (k = i - 1; k >= 0; --k) { mul = M[k][i]; for (n = i; n < N * 2; ++n) { M[k][n] = M[k][n] - mul * M[i][n]; } } } printf("求める逆行列は\n"); for(i=0;i<N;i++){ for(j=N;j<2*N;j++){ printf("%f",M[i][j]); } } return 0; }
###補足情報(言語/FW/ツール等のバージョンなど)
gcc -oでコンパイルしています。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2016/12/02 09:24