یک برنامهٔ C++ بنویسید که انتگرال عددی تابع f(x)=sin(x) را روی بازهٔ [a,b] با روش ذوزنقه و با تقسیم‌بندی n قسمت محاسبه کند و مقدار تقریبی را چاپ کند.

برای بهبود دقت و کارایی، مقدار h و مقدار f(a) و f(b) را قبل از حلقه محاسبه کنید تا از محاسبه مکرر جلوگیری شود. همچنین اگر به سرعت اجرای بالاتر در محیط‌های HPC نیاز دارید، می‌توانید با درج #pragma omp parallel for reduction(+:sum) حلقه جمع را موازی کنید و بررسی صحت فرمان‌ها را با کاهش داده‌های به اشتراک گذاشته انجام دهید. مطمئن شوید ورودی‌ها a، b، n معتبر هستند (n>0) و از double برای همه مقادیر استفاده کنید تا سازگاری نوع داده حفظ شود. در نهایت خروجی مقدار تقریبی انتگرال است و می‌تواند همراه با ریز زمان اجرا چاپ شود.