integralmontecarlo2d.cpp

何光渝所著《计算机常用数值算法与程序(C++版)》的随书源代码。 是想学习C++科学计算的人士的必看代码。

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