实现复合辛普森公式的简单C语言程序
#include
#include
#define f(x) (x)/(4+(x)*(x))
main()
{int k,n;
float a,b,h,s,x;
scanf("%f,%f,%d",&a,&b,&n);
h=(b-a)/(2*n);
x=a;
s=f(x)-f(b);
k=1;
do
{x=x+h;s=s+4*f(x);
x=x+h;s=s+2*f(x);
k++;
}while(k<=n);
s=(s*h)/3;
printf("s=%f\n",s);
}
1