یک برنامه پایتون بنویس که با استفاده از تابع random.randint تعداد n نمونه صحیح تصادفی در بازه [a, b] تولید کند و با روش مونت‌کارلو مقدار تقریبی انتگرال تابع f(x) = x**2 را در همان بازه محاسبه و نمایش دهد.

توجه: random.randint به صورت گسسته از بازه [a,b] نمونه می‌گیرد و این برای محاسبه انتگرال پیوسته نتیجه‌ای دقیق را تضمین نمی‌کند؛ اگر می‌خواهید نمونه‌گیری پیوسته باشد از random.uniform استفاده کنید یا از روش‌های گسسته با Δx مناسب استفاده کنید. برای بازه [a,b] و تابع f(x)=x^2، تقریب انتگرال برابر است با (b-a) ضربدر میانگین مقادیر نمونه‌برداری‌شده و با افزایش n خطای نمونه‌ای کاهش می‌یابد. برای بهبود سرعت می‌توانید از numpy برای محاسبات برداری بهره ببرید یا از multiprocessing/numba استفاده کنید و با seedکردن مقدار بازگشت‌پذیری را تضمین کنید.