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

برای پیاده‌سازی با decimal، مقدار دقت را با getcontext().prec تنظیم کن و ورودی‌های a، b و n را به Decimal تبدیل کن تا گام h و تمامی مقادیر میانی به‌طور دقیق محاسبه شوند. گام h به صورت Decimal محاسبه می‌شود و از طریق حلقهٔ ساده یا یک ژنراتور، جمع ف(a+i·h)ها را به‌همراه (f(a)+f(b))/2 انجام بده و نتیجه را با فرمول ذوزنقه‌ای ضرب کن تا مقدارIntegral به‌دست آید. برای تابع f، ورودی‌های Decimal را بپذیرد یا یک wrapper بساز تا خروجی Decimal بدهد تا دقت حفظ شود. همچنین به یاد داشته باش که هرچقدر n بزرگ باشد، زمان اجرا افزایش می‌یابد و در کاربردهای HPC معمولاً از برداری‌سازی یا روش‌های مناسب دیگر استفاده می‌شود.