//IntegralSimpson2D.cpp		 二重變步長辛卜生法求積

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

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

	cout << "dValue = " << dValue << endl << endl;
}

// 計算上、下限y1(x)、y0(x)值
template <class _Ty>
void FunctionBoundaryIS2D(_Ty x, valarray<_Ty>& y)
{
	y[0] = -sqrt(1.0 -x *x);
    y[1] = -y[0];
}

// 計算被積f(x,y)函數值
template <class _Ty>
_Ty FunctionValueIS2D(_Ty x, _Ty y)
{
	_Ty zValue;

	zValue = exp(x * x + y * y);
    
	return(zValue);
}