Ｃ言語のプログラムで質問です。

質問

Ｃ言語のプログラムで質問です。
下のプログラム（最小二乗法の計算）を実行したところ
－１．#IND00 というエラーが出てしまいます。
どこを直せばいいのでしょうか、教えてください。

#include <stdio.h>
#include <math.h>
/* gauss33.c */
#define N 3
main(){
double A[N][N],Aa[N][N];
double b[N],x[N], bb[N], e[N];
int n=N;
int i, j, k;
double akk, aik, s;

double y[N];
double xx,yy;

for(i=0;i<n;i++){ /*変数の初期化*/
x[i]=y[i]=0;
for(j=0;j<n;j++)
A[i][j]=0;
}

for(i=0;i<5;i++){ /*データ点は５点*/
printf("\n(x,y)=");
scanf("%lf,%lf",&xx,&yy);
A[0][0]+=xx*xx*xx*xx; /*Σx^4*/
A[0][1]+=xx*xx*xx; /*Σx^3*/
A[0][2]+=xx*xx; /*Σx^2*/
A[0][1]=A[1][0];
A[0][2]=A[1][1]=A[2][0];
A[1][2]+=xx; /*Σx*/
A[1][2]=A[2][1];
A[2][2]=n;
y[0]+=xx*xx*yy; /*Σx^2y*/
y[1]+=xx*yy; /*Σxy*/
y[2]+=yy; /*Σy*/
}

/* save original coefficients */

for(i=0; i<n; i++){
for(j=0; j<n; j++){
Aa[i][j]=A[i][j]; }
bb[i]=b[i];
}

/* forward operation */

for(k=0; k<n-1; k++){
akk=1/A[k][k];
for (i=k+1; i<n; i++){
aik=-A[i][k]*akk;
for (j=k+1; j<n; j++){
A[i][j]+=aik*A[k][j]; }
b[i]+=aik*b[k];
}
for(j=k+1; j<n; j++){
A[k][j]*=akk; }
b[k]*=akk;
}

/* backward operation */

x[n-1]=b[n-1]/A[n-1][n-1];
for(k=n-2; k>=0; k--){
s=0.0;
for (j=k+1; j<n; j++){ s+=A[k][j]*x[j]; }
x[k]=b[k]-s;
}

/* chek */

for(i=0; i<n; i++){
s=0.0;
for(j=0; j<n; j++){
s+=Aa[i][j]*x[j];}
e[i]=s-bb[i];

printf("\nx(%d)=%f error=%f\n",i, x[i], e[i]); }
}

回答No.2

asuncion · ベストアンサー

代入文は、右辺を左辺に代入するために使うことを
理解できていますか？

今のコードには、せっかく計算した結果を
直後の代入文で不用意に上書きしている箇所が
散見されます。

お礼

代入文があまり理解できてなかったようです。
直したら上手くいきました。ありがとうございます。

回答No.1

magicalpass

ロジックを追い切れていないので、全体的にどこがどうとは言えませんが、
パッと見たところでも代入の右辺と左辺の配置がおかしそうなところが見受けられます。

> A[0][1]=A[1][0];
> A[0][2]=A[1][1]=A[2][0];

> A[1][2]=A[2][1];

このあたり、一番右辺の項にはループ内で１回も代入が無いし、他の項もこの式しかないから、
ループをいくら回しても値が初期化した０のまま変わりません。
これが最終的に割り算の分母に入ってきて除算エラーになっているのでしょう。

基本に戻ってソースを見直してみてはいかがですか？

数式等のコメントでもあれば、もう少しアドバイスできるかもしれませんが。

お礼

上手くいきました。代入の方法がいまいち分かっていなかったようです。
有難うございます。