یک برنامه به زبان C بنویسید که با استفاده از srand از stdlib.h و تابع rand یک محاسبهٔ مونت‌کارلو برای تخمین عدد π انجام دهد: N جفت نقطه تصادفی در مربع [0,1]×[0,1] تولید کند، تعداد نقاط داخل چارک واحد را بشمارد و مقدار تقریبی π را چاپ کند.

برای تولید اعداد تصادفی در بازه [0,1)، از (double)rand() / RAND_MAX استفاده کنید؛ اما بعضی پیاده‌سازی‌ها توزیع کامل [0,1) را تضمین نمی‌کنند. برای پوشش دقیق‌تر بازه، می‌توانید از (double)rand() / (RAND_MAX + 1.0) استفاده کنید و همواره srand(time(NULL)) را قبل از استفاده فراخوانی کنید. در حلقه N بار دو عدد x و y تولید کنید و اگر x*x + y*y <= 1، شمارنده داخل چارک واحد را افزایش دهید و در پایان pi ≈ 4.0 * inside / N را چاپ کنید. هر چه N بزرگ‌تر باشد، خطای آماری تقریب pi کمتر و به ∼1/√N می‌رسد.