یک برنامهٔ C++ بنویسید که تعداد قسمت‌های تقسیم N را از آرگومان خط فرمان با std::atoi خوانده و با روش مستطیل ساده مقدار تقریبی π را از انتگرال ∫_0^1 4/(1+x^2) dx محاسبه و چاپ کند.

فایل را با #include و استفاده از std::atoi برای تبدیل argv[1] به int بنویسید؛ سپس اگر N>0 نبود پیام خطا نمایش دهید. برای محاسبه از double استفاده کنید: h = 1.0 / N، در یک حلقه از i=0 تا N-1 میانهٔ هر زیربازه x = (i+0.5)*h را محاسبه کرده و جمعِ 4.0/(1+x*x) را بر حسب h ضرب کنید تا π حاصل شود. نکات: ورودی را اعتبارسنجی کنید، از double برای دقت بیشتر استفاده کنید و برای بهبود کارایی در برنامه‌های بزرگ می‌توانید حلقه را بعداً با OpenMP موازی‌سازی کنید.