用梯形法求f(x)=x^2+1在[0,10]之间的定积分#include "math.h"double integ(double a,double b) { double s,x,h; int n=100,i; h = fabs(b-a)/n; s = (sin(a)+sin(b))/2.0; for (i=1; i<=n-1; i++) { x = a+i*h; s = s+sin(x); } s = s*h; return s; }main() { double s; s = integ(0.0,10.0); printf("s=%f\n",s); }
{ double s,x,h;
int n=100,i;
h = fabs(b-a)/n;
s = (sin(a)+sin(b))/2.0;
for (i=1; i<=n-1; i++)
{ x = a+i*h;
s = s+sin(x); }
s = s*h;
return s;
}main()
{ double s;
s = integ(0.0,10.0);
printf("s=%f\n",s);
}