C言語なら、俺に聞け! <9>
■ このスレッドは過去ログ倉庫に格納されています
0331名無しさん@お腹いっぱい。
NGNG#include <math.h>
int N;
double M=10.0+h,h,pi=M_PI; ←ここのhと
void integ();
double f(double);
main()
{
h=1.5;
integ();
h=0.5;
integ();
void integ() ←ここが
{
double daikei=0.,simpson=0.;
int i;
N=M/h;
for(i=-N;i<=N;i++){daikei +=f(i*h);}
daikei *=h;
for(i=-N;i<=N;i+=2){simpson +=4.*f(i*h);}
for(i=-N;i<=N-1;i+=2){simpson +=2.*f(i*h);}
simpson *=h/3.0;
printf("\n h=%lf\n",h);
printf("sqrt(pi)=%15.10e\n",sqrt(pi));
printf("daikei =%15.10e\n",daikei);
printf("simpson =%15.10e\n",simpson);
}
}
double f(double x)
{
return exp(-x*x);
}
コンパイルしてもうまくいかないんですけど、どうしてでしょうか?
よろしかったら教えてください。
■ このスレッドは過去ログ倉庫に格納されています