This is a test
double SimpleMonteCarlo() { double rangeMin = 0; double rangeMax = 3.14159265359; size_t numSamples = 10000; std::random_device rd; std::mt19937 mt(rd()); std::uniform_real_distribution dist(rangeMin, rangeMax); double ySum = 0.0; for (size_t i = 1; i <= numSamples; ++i) { double x = dist(mt); double y = sin(x)*sin(x); ySum += y; } double yAverage = ySum / double(numSamples); double width = rangeMax - rangeMin; double height = yAverage; return width * height; }
That was a test