Не более чем иллюстрации для - применение метода Монте-Карло к вычислению объема шара единичного радиуса (двойной интеграл по области x^2+y^2 <= 1 функции 2*sqrt(1-x^2-y^2)).
Для наглядности она из полученного значения вычисляет число "пи".
#include <random>
#include <iostream>
using namespace std;
int main(int argc, const char * argv[])
{
default_random_engine e;
uniform_real_distribution r(-1.0,1.0);
unsigned long long total = 0, vol = 0;
for(unsigned long long i = 1; ; ++i)
{
double x = r(e), y = r(e), z = r(e);
++total;
if (x*x+y*y > 1.0) continue; // Выход за область интегрирования
if (z*z < 1 - x*x - y*y) ++vol; // Точка входит в область
if (i%1000000 == 0)
cout << 6.0*vol/total << endl;
}
}
Запускайте и смотрите сами :)