با استفاده از تابع std::rand از هدر <cstdlib> برنامه‌ای بنویسید که با روش مونت‌کارلو مقدار عدد π را با تولید N نقطه تصادفی در مربع واحد محاسبه کند و مقدار تقریبی به‌همراه خطای نسبی را نمایش دهد.

از هدر <cstdlib> و تابع std::rand برای تولید اعداد تصادفی استفاده کنید؛ ابتدا با std::srand(time(nullptr)) بذر را مقداردهی کنید، سپس در حلقه‌ای N بار مختصات x و y را با تبدیل خروجی std::rand به double در بازه [0,1) بسازید (مثلاً double(x)/RAND_MAX). تعداد نقاطی که به داخل ربع‌دایره با شعاع 1 وارد می‌شوند را بشمارید و با نسبت شمار به N و ضرب در 4 مقدار π را تقریب بزنید. برای دقت بهتر N را افزایش دهید و خطای نسبی را با |approx - M_PI| / M_PI محاسبه کنید. در زمینه محاسبات علمی و HPC می‌توانید اشاره کنید که std::rand برای کاربردهای موازی مناسب نیست و برای مقیاس بالا از تولیدکننده‌های تصادفی با کیفیت بالاتر یا تکنیک‌هایی مانند OpenMP/رویه‌های تولید عدد تصادفی موازی استفاده کنید.