integraltrapeziaselfadapt.cpp

来自「计算机常用数值计算算法与程序 C++版 各章代码」· C++ 代码 · 共 31 行

//IntegralTrapeziaSelfAdapt.cpp		
//自适应梯形法求积

#include <iostream>			//输入输出流头文件
#include "Integral.h"		//积分头文件
using namespace std;		//名字空间

void main(void)
{	
	double eps = FLOATERROR, a(-1), b(1), d(0.0001);
            
	cout << "IntegralTrapeziaSelfAdapt()" << endl << endl;

	double dIntegralValue = IntegralTrapeziaSelfAdapt(a,b,eps,d);

	cout << "IntegralValue = " << dIntegralValue << endl << endl;
}

// 计算被积函数值
template <class _Ty>
_Ty FunctionValueITSA(_Ty x)
{
	_Ty y;

	y = 1.0 / (1.0 + 25.0 * x * x);

    return(y);
}