for(i=1;fabs(s1-s2)>e;i++)
    {
        double s=0;
        for(k=i;k<2*i;k++)
            s+=2*f(a+(k-i+1.0/4)*h)-f(a+(k-i+1.0/2)*h)+2*f(a+(k-i+3.0/4)*h);
        s2=0.5*s1+h*s/6.0;
        if(fabs(s1-s2)>e)
        {
            s1=s2;
            h=h/2.0;
        }
    }

#include<stdio.h>
#include<math.h>
double f(double x)
{
    return 1/(x*x*x);
}
void main()
{
    int i,k;
    double h,a,b,s1,s2=0,e=1.0e-6;
    a=1;b=2;h=b-a;
    s1=h*(f(a)+4*f((a+b)/2)+f(b))/6.0;
    for(i=1;fabs(s1-s2)>e;i++)
    {
        double s=0;
        for(k=i;k<2*i;k++)
            s=s+2*f(a+(k-i+1.0/4)*h)-f(a+(k-i+1.0/2)*h)+2*f(a+(k-i+3.0/4)*h);
        s2=0.5*s1+h*s/6.0;
        if(fabs(s1-s2)>e)
        {
            s1=s2;
            h=h/2.0;
        }
    }
            
    printf("%lf %lf\n",s1,s2);
}

这是在VC环境下运行的，应该没错吧