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;
        }
    }

楼主可否发一个稍微完整点的程序啊，否则实在是看不懂！！

哇哦，这个表达式看起来就很爽，大兵沃恩你怎么看？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);

怎么没有看到关于s1的语句呢?程序应该不完整吧

呵呵，等他把代码都发出来可能更爽！！

还有这样问问题的？什么都不说清楚还想知道答案！

哎，楼主敢发全部代码么、

试想让你看这么一段，而且没有注释。

楼主什么赶脚？

反正我一看那个最长的一行就不想看了。。

特别是在没有注释的情况下。。

#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环境下运行的，应该没错吧