در پایتون تابعی بنویسید که با استفاده از روش ذوزنقه‌ای (Trapezoidal rule) عددی انتگرال یک تابع دلخواه را روی بازه [a, b] با n بخش محاسبه کند و با استفاده از numpy برداریزه‌سازی را انجام دهد.

برای بهره‌گیری از برداری، x را با np.linspace(a, b, n+1) بسازید و f را روی آرایه x اعمال کنید تا مقادیر f(x) به دست آید. نتیجه نهایی با فرمول ذوزنقه‌ای: h * (0.5*f[0] + f[1:-1].sum() + 0.5*f[-1]) باشد که h=(b-a)/n است. برای کارایی، از dtype=float64 استفاده کنید و تابع خود را به صورت برداری پیاده‌سازی کنید یا در صورت لازم از np.vectorize سود ببرید؛ اگر f با آرایه‌ها سازگار نیست، تبدیل آن به آرایه را انجام دهید. برای اعتبارسنجی، نتیجه را با scipy.integrate.quad مقایسه کنید یا با افزایش n همگرایی را بررسی کنید و در صورت نیاز از numba یا موازی‌سازی برای سرعت بیشتر استفاده کنید.