//IntegralSimpson1D.cpp		 變步長辛卜生法求積

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

void main(void)
{	
	double eps = FLOATERROR, a(0), b(1);
                
	cout << "IntegralSimpson1D()" << endl << endl;

	double t = IntegralSimpson1D(a, b, eps);

	cout << "t = " << t << endl << endl;
}

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

	y = log(1.0 + x) / (1.0 + x * x);
    
	return(y);
}