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

برای حل، از هدر و std::accumulate برای جمع‌زدن مقادیر تابع در نقاط داخلی استفاده کنید: گام h = (b-a)/n را محاسبه کرده، مقادیر f(x_i) را در یک بردار ذخیره کنید و سپس جمع داخلی‌ها را با std::accumulate بگیرید و فرمول ذوزنقه را اعمال کنید: integral ≈ h*( (f(a)+f(b))/2 + sum_internal ). از double برای دقت بهتر استفاده کرده و می‌توانید f را به‌صورت یک lambda تعریف کنید؛ در زمینه محاسبات علمی و HPC توجه به تخصیص حافظه (reserve برای بردار) و دقت عددی مهم است.