变步长复合梯形公式算法公式

#include "stdio.h" 
#include "math.h" 
double f(double x) 
{ 
double y; 
y=4/(1+x*x);
return(y); 
} 
main() 
{ 
double a,b,e,h,k,s,t[2000]; 
int n=1,m=0; 
printf("please input a,b,e\n"); 
scanf("%lf%lf%lf",&a,&b,&e); 

t[1]=(b-a)*(f(a)+f(b))/2; 
do 
{ 
h=(b-a)/n; 
s=0; 
for(k=0;k<n;k++) 
s+=f(a+(k+0.5)*h); 
t[2*n]=t[n]/2+h/2*s; 
n=n*2; 
m=m+1; 
} 
while(fabs(t[n]-t[n/2])>0.0000001); 
printf("%d\n",m); 
printf("f(x)*dx=%.7lf",t[n]);