یک برنامه‌ی ساده به زبان C++ بنویسید که مقدار انتگرال تقریبی تابع f(x)=sin(x^2) را روی بازه [0,1] با استفاده از روش مونت‌کارلو و تولید اعداد شبه‌تصادفی از هدر <random> محاسبه کند.

برای این کار از std::mt19937 با بذر گرفته‌شده از std::random_device و یک توزیع یکنواخت real در [0,1] استفاده کنید تا N نمونه تصادفی بسازید. مقدار f(x) = sin(x^2) را برای هر نمونه محاسبه کنید و میانگین این مقادیر را بگیرید؛ از آنجا مقدار انتگرال روی [0,1] تقریباً برابر با میانگین ضرب در طول بازه است (در این بازه طول بازه برابر با 1 است). هر چه N بزرگ‌تر باشد، خطای تقریبی ~1/√N کاهش می‌یابد و برای دقت بهتر از double استفاده کنید. برای بهبود کارایی می‌توانید از موازی‌سازی ساده مانند OpenMP یا threading بهره ببرید تا محاسبات نمونه‌ها به‌صورت هم‌زمان انجام شود.