//IntegralMonteCarlo2D.cpp	 多重蒙特卡洛法求積

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

void main(void)
{	
	double a[3] = {1, 1, 1};
    double b[3] = {2, 2, 2};

	valarray<double> da(a,3);
	valarray<double> db(b,3);
    
	double dValue = IntegralMonteCarlo2D(da, db);	//調用求積函數
	    
	cout << "IntegralMonteCarlo2D()" << endl << endl;

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

//計算被積函數值
template <class _Ty>
_Ty FunctionValueMC2D(int n, valarray<_Ty>& x)		
{
	_Ty dValue(0);

	for(int i=0; i<n; i++) dValue = dValue + x[i] * x[i];

    return(dValue);
}