با زبان C++ تابعی بنویسید که با استفاده از std::accumulate و قاعده ذوزنقه انتگرال عددی یک تابع را روی بازه [a, b] با مقادیر نمونه‌شده در یک std::vector محاسبه کند.

روش حل: با فرض داشتن برداری از مقادیر f(x) که هم‌ردیف نقاط روی بازه [a,b] نمونه‌برداری شده‌اند، گام h را برابر (b-a)/(n-1) محاسبه کنید و مجموع مقادیر داخلی را با std::accumulate روی بازه iterators از begin()+1 تا end()-1 بگیرید؛ سپس مقدار انتگرال برابر h * ( (f0 + fn)/2 + sum_interior ) خواهد بود. در C++ هدر مورد نیاز و است و از double برای دقت استفاده کنید؛ برای کارهای بزرگِ محاسباتی می‌توانید بعدها به جای std::accumulate از std::reduce با execution::par یا تکنیک‌های OpenMP برای اجرای موازی فکر کنید تا عملکرد بهتری در محیط‌های HPC داشته باشید.