trapzint.cpp

清华大学郑莉《C++语言程序设计》课件、C++教材源代码、实验参考程序全集

//Trapzint.cpp	//文件二，类实现
#include "Trapzint.h"	//包含类的声明头文件
double Trapz::operator ()(double a,double b,double eps) const
{			//积分运算过程，重载为运算符()
	int done(0);	//是Trapz类的虚函数成员
	int n;
	double h,Tn,T2n;
	n=1;
	h=b-a;
	Tn=h*(f(a)+f(b))/2.0;	//计算n=1时的积分值
	while(!done)
	{
		double temp(0.0);
		for(int k=0;k<n;k++)
		{
		      double x=a+(k+0.5)*h;
		      temp+=f(x);
		}
		T2n=(Tn+h*temp)/2.0;	//变步长梯形法计算
		if(fabs(T2n-Tn)<eps) done=1;	//判断积分误差
		else	//进行下一步计算
		{
		      Tn=T2n;
		       n*=2;
		       h/=2;
		}
	}
	return T2n;
}
//End of file Trapzint.cpp