//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);
}