自适应梯形公式（变步长）.cpp

来自「提供常用数值计算的c++代码源程序。方便初学者进行练习！」· C++ 代码 · 共 46 行

#include<iostream.h>
#include<math.h>
int n;
float f(float);
void main()
{
	float s;
	float AutoTrap(float (*)(float),float,float);
	s=AutoTrap(f,0.0,1.0);
	cout<<"T("<<n<<")="<<s<<endl;
}
float AutoTrap(float (*)(float),float a,float b)
{
	int i;
	float x,s,h=b-a;
	float t1,t2=h/2.0*(f(a)+f(b));
	n=1;
	do{
		s=0.0;
		t1=t2;		
		for(i=0;i<=n-1;i++)
		{
			x=a+i*h+h/2;
			s+=f(x);
		}
		t2=(t1+s*h)/2.0;
		n*=2;
		h/=2.0;
	}while( fabs(t2-t1)>1e-6);
	return t2;
}

float f(float x)
{
	return 1/(1+x*x);
}
//用变步长梯形法计算积分积分只要将函数定义为

/*float f(float x)
{
	if(x==0)
		return 1;
	else 
		return sin(x)/x;
}*/