انتگرالگیری به روش سیمپسون

برای تقریب انتگرال یک تابع با روش سیمپسوت یک سوم در متلب به کد زیر توجه کنید:

در کد بالا داریم:

• بازه [a,b] به n زیربازه تقسیم شده است، بنابراین اندازه گام h = b-a/n است.
• مجموع با اولین و آخرین مقادیر تابع در نقاط پایانی a و b شروع می‌شود.
• حلقه روی نقاط داخلی بازه تکرار می‌شود. بسته به اینکه اندیس i زوج یا فرد باشد، وزن مربوطه (۴ یا ۲) به مقدار تابع در آن نقطه اعمال می‌شود.
• پس از حلقه، نتیجه در h/3 ضرب می‌شود تا تقریب انتگرال نهایی به دست آید.

با اجرا در پنجره فرمان متلب، برای یک مثال دلخواه خروجی به صورت زیر خواهد بود: