C言語を用いた45×４５の逆行列の表示について・・・・。。

質問

次のプログラムを動かすと逆行列を求めることができます。
１０×１０の計算までは、スムーズに計算する事ができたのですが、
私が知りたい４５×４５の行列の計算になるとなぜかうまくいきません。
そこでみなさんにお聞きしたいのは、計算がうまくいかない理由がプログラムにあるのか、単にＰＣの限界なのか、それ以外なのか。みなさんからアドバイスをいただけないでしょうか？
※下のプログラム、見づらくてすいません。

#include <stdio.h>
#define N 10 /* 元の数 */
#define N2 2*N

/* プロトタイプ宣言 */
void Print(int k);

/* 大域変数 */
double a[N][N2]=
{/* 行列 Aおよび単位行列 */
{
ここに45×45の正則行列＆45×45の単位行列を入れます・・・。。
};
int main(void)
{
double d;
int i,j,k;
for (k=0;k<=N-1;k++)
{
Print(k);
for(i=0;i<=N-1;i++)
{
if(i!=k)
{
d=a[i][k]/a[k][k];
for (j=k+1;j<=N2-1;j++)
a[i][j]=a[i][j]-a[k][j]*d;
a[i][k] = 0.0;
}
}
d=a[k][k];
for (j=k;j<=N2-1;j++)
a[k][j]=a[k][j]/d;
}
Print(N);
}

void Print(int k)
{
int i,j;
printf("\n ステップ %d \n",k); /* 小見出し */
for(i=0;i<=N-1;i++)
{ /* 第ｉ行を印刷 */
for(j=0;j<=N-1;j++)
printf("%6.4f,",a[i][j]);
printf(" : ");
for(j=N;j<=N2-1;j++)
printf("%6.4f,",a[i][j]);
printf("\n");
}

}

回答No.2

a-saitoh

ガウスの消去法でしょうか？だとしたらO(N^3)の計算量なので、N=10とN=45では所要時間が100倍くらい違うと予想できます。メモリ容量的にも余裕ですし、PCの限界ではありません。

回答No.1

foobar

少なくとも、2行目の
#define N 10
は、変える必要がありそうに思います。