یک برنامهٔ C++ بنویسید که با استفاده از std::rand (از هدر ) و روش مونت‌کارلو مقدار π را تخمین بزند و نتیجه را برای تعداد نمونهٔ دلخواه چاپ کند.

برای تولید عدد تصادفی در [0,1) بهتر است از تبدیل double(std::rand()) / (RAND_MAX + 1.0) استفاده کنید تا بازتولید 1.0 دست کم باشد. قبل از هر چیز با std::srand(std::time(nullptr)) مقداردهی اولیه انجام دهید و این کار را در ابتدای برنامه قرار دهید. هر نقطه (x,y) را از مربع [0,1]×[0,1] با x*x + y*y <= 1 بررسی کنید؛ نسبت نقاط داخل به کل، pi/4 است و در نتیجه pi ≈ 4*(داخل/کل). برای دقت و سرعت، از تعداد نمونه‌های بزرگ استفاده کنید و به‌سادگی می‌توانید با چند نخ موازی‌سازی کنید (مثلاً با OpenMP یا std::thread)، و اگر دقت بیشتری می‌خواهید، به جای std::rand از std::mt19937 به‌همراه توزیع یکنواخت استفاده کنید.