روش سیمپسون

انتگرال‌گیری

فرض کنید هدف محاسبه‌ی تقریبی انتگرال ∫xixi+2f(x)dx\int_{x_i}^{x_{i+2}} f(x) dx است. f(x)f(x) را با چندجمله‌ای درونیاب آن در نقاط xix_i و xi+1x_{i+1} و xi+2x_{i+2} تقریب می‌زنیم، لذا داریم:

∫xixi+2f(x)dx≈h∫01[fi+sΔfi+s(s−1)2Δ2fi]ds=h[fi+23fi+1+16fi+2]\int_{x_i}^{x_{i+2}} f(x) dx \approx h \int_{0}^{1} \left[ f_i + s \Delta f_i + \frac{s(s-1)}{2} \Delta^2 f_i \right] ds = h \left[ f_i + \frac{2}{3} f_{i+1} + \frac{1}{6} f_{i+2} \right]

بنابراین:

∫xixi+2f(x)dx≈h3(fi+4fi+1+fi+2)\int_{x_i}^{x_{i+2}} f(x) dx \approx \frac{h}{3} \left( f_i + 4 f_{i+1} + f_{i+2} \right)

روش سیمپسون مرکب

برای محاسبه‌ی تقریبی انتگرال معین ∫abf(x)dx\int_{a}^{b} f(x) dx قرار می‌دهیم h=b−a2nh = \frac{b - a}{2n} و نقاط a=x0,x1,…,x2n=ba = x_0, x_1, \ldots, x_{2n} = b را در نظر می‌گیریم و در هر بازه‌ی [xi,xi+2][x_i, x_{i+2}] از رابطه‌ی قبل استفاده می‌کنیم. بنابراین نتیجه می‌شود:

∫x0x2nf(x)dx≈S(h)=h3(f0+4∑i=0n−1f2i+1+2∑i=1n−1f2i+f2n)\int_{x_0}^{x_{2n}} f(x) dx \approx S(h) = \frac{h}{3} \left( f_0 + 4 \sum_{i=0}^{n-1} f_{2i+1} + 2 \sum_{i=1}^{n-1} f_{2i} + f_{2n} \right)

فرمول به دست آمده را فرمول روش سیمپسون مرکب می‌نامند. نکته: روش سیمپسون را زمانی می‌توان به کار گرفت که تعداد نقاط به صورت 2n+12n + 1 باشد.

روش‌های مرتبط