C++博客-CoolWalfhttp://www.cppblog.com/CoolWalf/zh-cnTue, 09 Jun 2026 23:57:52 GMTTue, 09 Jun 2026 23:57:52 GMT60c++实现变步长积分http://www.cppblog.com/CoolWalf/archive/2008/10/19/64398.htmlCoolWalfCoolWalfSun, 19 Oct 2008 03:42:00 GMThttp://www.cppblog.com/CoolWalf/archive/2008/10/19/64398.htmlhttp://www.cppblog.com/CoolWalf/comments/64398.htmlhttp://www.cppblog.com/CoolWalf/archive/2008/10/19/64398.html#Feedback0http://www.cppblog.com/CoolWalf/comments/commentRss/64398.htmlhttp://www.cppblog.com/CoolWalf/services/trackbacks/64398.html#include<math.h>
double f(double x)
{
 if(x==0)
  return 1;
 else
 return sin(x)/x;
}
double ff(double a,double b,double c)
{
 double h=b-a;
 double t1,t2=0,x,s=0;
 t1=(h/2)*(f(a)+f(b));
 while(fabs(t2-t1)>c)
 {
  s=0;x=a+h/2;
  t1=t2;
  while(x<b)
  {
   s+=f(x);x+=h;
  }
  t2=t1/2+(h/2)*s;
  h=h/2;
 }
    return t2; 
}
void main()
{
 double a,b,c;
 cout<<"请分别输入上下限"<<endl;
 cin>>a>>b;
 cout<<"请输入误差限度："<<endl;
  cin>>c;
  cout<<ff(a,b,c)<<endl;
}

CoolWalf 2008-10-19 11:42 发表评论
]]>