//IntegralRomberg.cpp	 龍貝格法求積

#include <iostream>		//輸入輸出流頭文件
#include "Integral.h"	//積分頭文件
using namespace std;	//名字空間

void main(void)
{	
	double eps = FLOATERROR, a(0), b(1);

	cout << "IntegralRomberg()" << endl << endl;

	double dIntegralValue = IntegralRomberg(a,b,eps);

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

// 計算被積函數(shù)值
template <class _Ty>
_Ty FunctionValueR(_Ty x)
{
	_Ty y;

	y = x / (4.0 + x * x);

    return(y);
}