s1与s2为什么始终相等啊求教？？？？
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;
}
}
回复 2楼 我叫沃恩
#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);
}