با زبان C++ برنامه‌ای بنویسید که انتگرال عددی تابع f(x)=e^{-x^2} را روی بازه [a,b] با روش ذوزنقه‌ای (trapezoidal rule) و با n بخش محاسبه کرده و زمان اجرا را گزارش کند.

برای حل، مقدار a، b و n را از ورودی بگیرید، گام h = (b-a)/n را محاسبه کنید و با فرمول ذوزنقه مجموع را محاسبه نمایید: (h/2)*(f(a)+f(b)+2*sum_{i=1}^{n-1} f(a+i*h)). از نوع double برای دقت استفاده کنید و f را با std::exp(-x*x) پیاده‌سازی کنید. برای اندازه‌گیری زمان اجرا از std::chrono استفاده کنید. به‌عنوان نکات بهینه‌سازی می‌توانید حلقه جمع‌زنی را با OpenMP و گزینه reduction موازی کنید (با کامپایلر و فلگ مناسب مثل -fopenmp) و از دسترسی ترتیبی به حافظه برای کاهش کش میس‌ها استفاده نمایید.