- 169 名前：デフォルトの名無しさん [2007/11/09(金) 22:46:25 ]
- #include <stdio.h>
#include <math.h>
double f(double x){
return 1/sqrt(1+pow(x,4)); /*求める定積分*/
}
int main(void){
double x,sum=0,h;
printf("h:"); scanf("%lf",&h);
for(x=2*h;x<1;x+=2*h){ /*範囲指定*/
sum+=(f(x) + 4*f(x-h) + f(x-2*h)) * h/3; /*シンプソンの公式*/
printf("%f %f\n",sum,x);
}
}
