یک برنامه پایتون بنویسید که با استفاده از تابع random.choice از ماژول random و روش مونت‌کارلو انتگرال تقریب زده شدهٔ تابع f(x)=sin(x) را روی بازه [0, π] با n نمونه محاسبه کند و مقدار تقریبی انتگرال، میانگین نمونه‌ها و خطای برآورد شده را چاپ نماید.

برای کارایی خوب، می‌توانید ابتدا آرایه‌ای از نقاط یکنواخت در [0, π] بسازید و با random.choice به‌صورت با جایگزینی n نمونه از آن انتخاب کنید (یا از random.choices برای انجام همان کار به‌طور مستقیم استفاده کنید). مقدار انتگرال تقریبی برابر است با طول بازه π ضربدر میانگین مقادیر sin(x) روی این نمونه‌ها. خطای برآورد با تقریب استاندارد انحراف نمونه تقسیم بر جذر n محاسبه می‌شود؛ برای تکرارپذیری seeds هر دو کتابخانه random و numpy را مقداردهی کنید. استفاده از numpy برای محاسبه میانگین و واریانس سرعت محاسبات را افزایش می‌دهد، اما همین روش ساده نیز با n کافی نتیجه دقیق می‌دهد.