ファイルが読み込めてない…？

[kohaA]

ex7.c

(819 バイト) ダウンロード数: 145 回

ラグランジュ多項式を解けるはず…なのですが出力されるファイルに書き込まれていません
ex7.txtには下のように数字が入っています
4 1
2 5
6 0
3 2
1 5

コード:

#include<stdio.h>
#include<math.h>
#define max 100

double lag(int , int , double , double*);

int main(void)
{

int i,n;
double x,y;
double datax[max];
double datay[max];

FILE *input;
FILE *output;

input=fopen("ex7.txt","r");
output=fopen("ex7ans.txt","w");

 for(i=0;i<max;i++){
   datax[i]=0.0;
   datay[i]=0.0;
 }

 i=0;

 while(fscanf(input,"%lf%lf",&datax[i],&datay[i])!=EOF){
   i++;
   n=i;
 }

 for(x=datax[0];x<datax[n-1];x+=0.1){
   y=0.0;
   for(i=0;i<n;i++){
     y=0.0;
     y+=datay[i]*lag(i,n,x,datax);
   }
   fprintf(output,"%f %f\n",x,y);
 }

 fclose(input);
 fclose(output);

 return 0;

}




double lag(int i,int n,double x,double datax[]){

  int j;
  double ans=1.0;

  for(j=0;j<n;j++){

    if(i!=j){
      ans*=(x-datax[j])/(datax[i]-datax[j]);
    }
  }
  return 0;
}

[box]

実験も何もしておらず、コードを見ただけですが、
36行目の
y=0.0;
は本当に必要なのか？という疑問があります。

[みけCAT]

実行はしていませんが、

ex7.txtには下のように数字が入っています
4 1
2 5
6 0
3 2
1 5

このデータを読み込むとdatax[0]は4、datax[n-1]は1になるので、最初から33行目のx<datax[n-1]が偽となり、
39行目のfprintfは1回も実行されず、その結果「出力されるファイルに書き込まれていません」という事象が観測されたのだと思います。

[みけCAT]

オフトピック

コード:

 while(fscanf(input,"%lf%lf",&datax[i],&datay[i])!=EOF){

不正な入力を与えると無限ループになるリスクがあるので、

コード:

 while(fscanf(input,"%lf%lf",&datax[i],&datay[i])==2){

の方がいいかもしれないと思いました。

[box]

別の疑問が…。
63行目の
return 0;
は本当に正しいでしょうか？
return ans;
のような気がするけど…。

[kohaA]

はいreturn 0;は間違いでした…

しかし未だ出力されず格闘中・・・

[みけCAT]

はいreturn 0;は間違いでした…

しかし未だ出力されず格闘中・・・

試しに入力を

コード:

4 1
2 5
6 0

にしたところ、(正しさは知りませんが)出力が出ました。
やはりこれが原因でしょう。

このデータを読み込むとdatax[0]は4、datax[n-1]は1になるので、最初から33行目のx<datax[n-1]が偽となり、
39行目のfprintfは1回も実行されず、その結果「出力されるファイルに書き込まれていません」という事象が観測されたのだと思います。

[kohaA]

お二人とも有難うございました
精度があっているのかは明日の解説で理解しようと思います

返戻値を間違えるくらいはいい加減気をつけます。。。

みけCATさんには毎度お世話になってます・・・なさけなひ