قاعدة سمبسون (بالإنجليزية : Simpson's rule ) في التحليل العددي هي طريقة من طرق التكامل العددي و هي في الحقيقة حالة خاصة من صيغ نيوتن-كوتس المغلقة لتقريب تكامل الدالة f باستخدام كثيرة الحدود التربيعية وهي طريقة محسنة لطريقة شبه المنحرف كما أنها أسرع تقارباً وأدق ويفسر ذلك من خلال أن قاعدة سمبسون تحتوي على نقطة المنتصف التي توفر توازن أفضل للتقريب. لأنه كلما زادت عدد التقسيمات في الفترة الجزئية كانت الطريقة أدق.
و تحسب صيغة شبه المنحرف القيمة الفعلية للتكامل عندما تكون f دالة كثيرة حدود من الدرجة الأولى على الأكثر . بينما صيغة سمبسون فإنها تحسب القيمة الفعلية للتكامل إذا كانت f دالة كثيرة الحدود من الدرجة الثالثة أو أقل.[ 1]
في الواقع تعتبر صيغة سمبسون من أكثر الصيغ استخداما حيث انها تستخدم على نطاق واسع لحل المسائل التطبيقية التي تتضمن تكاملات محدودة لدقتها الحسابية وسهولة استخدامها.
يمكن أن تستمد قاعدة سمبسون بتقريب المساحة المطلوبة تحت المنحنى الدالة f باللون الأزرق بالمساحة تحت منحنى حدودية من الدرجة الثانية لP باللون الأحمر
يمكن استنتاج قاعدة سمبسون بطرق مختلفة.
ولكي نحصل على صيغة سمبسون لابد بالتعويض في الصيغة الإستكمالية
P
2
(
x
)
=
c
∘
+
c
1
x
+
c
2
x
2
{\displaystyle P_{2}(x)=c_{\circ }+c_{1}x+c_{2}x^{2}}
∫
a
b
f
(
x
)
d
x
=
∫
x
∘
x
2
[
(
x
−
x
1
)
(
x
−
x
2
)
(
x
∘
−
x
1
)
(
x
∘
−
x
2
)
f
(
x
∘
)
+
(
x
−
x
∘
)
(
x
−
x
2
)
(
x
1
−
x
∘
)
(
x
1
−
x
2
)
f
(
x
1
)
+
(
x
−
x
∘
)
(
x
−
x
1
)
(
x
2
−
x
∘
)
(
x
2
−
x
1
)
f
(
x
2
)
]
d
x
+
∫
x
∘
x
2
(
x
−
x
∘
)
(
x
−
x
1
)
(
x
−
x
2
)
6
f
(
3
)
(
ξ
(
x
)
)
d
x
{\displaystyle \int _{a}^{b}f(x)dx={\int _{x_{\circ }}^{x_{2}}}[{\frac {(x-x_{1})(x-x_{2})}{(x_{\circ }-x_{1})(x_{\circ }-x_{2})}}f(x_{\circ })+{\frac {(x-x_{\circ })(x-x_{2})}{(x_{1}-x_{\circ })(x_{1}-x_{2})}}f(x_{1})+{\frac {(x-x_{\circ })(x-x_{1})}{(x_{2}-x_{\circ })(x_{2}-x_{1})}}f(x_{2})]dx+{\int _{x_{\circ }}^{x_{2}}}{\frac {(x-x_{\circ })(x-x_{1})(x-x_{2})}{6}}f^{(3)}(\xi (x))dx}
f
(
x
)
=
f
(
x
1
)
+
f
′
(
x
1
)
(
x
−
x
1
)
+
f
″
(
x
1
)
2
(
x
−
x
1
)
2
+
f
‴
(
x
1
)
6
(
x
−
x
1
)
3
+
f
(
4
)
(
ξ
(
x
)
24
(
x
−
x
1
)
4
{\displaystyle f(x)=f(x_{1})+f'(x_{1})(x-x_{1})+{\frac {f''(x_{1})}{2}}(x-x_{1})^{2}+{\frac {f'''(x_{1})}{6}}(x-x_{1})^{3}+{\frac {f^{(4)}(\xi (x)}{24}}(x-x_{1})^{4}}
and
∫
x
∘
x
2
f
(
x
)
d
x
−
[
f
(
x
1
)
(
x
−
x
1
)
+
f
′
(
x
1
)
2
(
x
−
x
1
)
2
+
f
″
(
x
1
)
6
(
x
−
x
1
)
3
+
f
‴
(
ξ
(
x
)
24
(
x
−
x
1
)
4
]
x
∘
x
2
+
1
24
∫
x
∘
x
2
f
(
4
)
(
ξ
(
x
)
)
(
x
−
x
1
)
4
d
x
{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx-\left[f(x_{1})(x-x_{1})+{\frac {f'(x_{1})}{2}}(x-x_{1})^{2}+{\frac {f''(x_{1})}{6}}(x-x_{1})^{3}+{\frac {f'''(\xi (x)}{24}}(x-x_{1})^{4}\right]_{x_{\circ }}^{x_{2}}+{\frac {1}{24}}\int _{x_{\circ }}^{x_{2}}f^{(4)}(\xi (x))(x-x_{1})^{4}dx}
1
24
∫
x
∘
x
2
f
(
4
)
(
ξ
(
x
)
)
(
x
−
x
1
)
4
d
x
=
f
(
4
)
(
ξ
1
)
24
∫
x
∘
x
2
(
x
−
x
1
)
4
d
x
=
[
f
(
4
)
(
ξ
1
)
120
(
x
−
x
1
)
5
]
x
∘
x
2
{\displaystyle {\frac {1}{24}}\int _{x_{\circ }}^{x_{2}}f^{(4)}(\xi (x))(x-x_{1})^{4}dx={\frac {f^{(4)}(\xi _{1})}{24}}\int _{x_{\circ }}^{x_{2}}(x-x_{1})^{4}dx=\left[{\frac {f^{(4)}(\xi _{1})}{120}}(x-x_{1})^{5}\right]_{x_{\circ }}^{x_{2}}}
لبعض قيم
ξ
1
{\displaystyle \xi _{1}}
(
x
∘
,
x
1
)
{\displaystyle (x_{\circ },x_{1})}
أيضا
h
=
x
2
−
x
1
=
x
1
−
x
∘
{\displaystyle h=x_{2}-x_{1}=x_{1}-x_{\circ }}
بالتالي :
(
x
2
−
x
1
)
2
−
(
x
∘
−
x
1
)
2
=
(
x
2
−
x
1
)
4
−
(
x
∘
−
x
1
)
4
=
0
{\displaystyle (x_{2}-x_{1})^{2}-(x_{\circ }-x_{1})^{2}=(x_{2}-x_{1})^{4}-(x_{\circ }-x_{1})^{4}=0}
(
x
2
−
x
1
)
3
−
(
x
∘
−
x
1
)
3
=
2
h
a
n
d
(
x
2
−
x
1
)
5
−
(
x
∘
−
x
1
)
5
=
2
h
5
{\displaystyle (x_{2}-x_{1})^{3}-(x_{\circ }-x_{1})^{3}=2hand(x_{2}-x_{1})^{5}-(x_{\circ }-x_{1})^{5}=2h^{5}}
∫
x
∘
x
2
f
(
x
)
d
x
=
2
h
f
(
x
1
)
+
h
3
3
f
″
(
x
1
)
+
f
(
4
)
(
ξ
1
)
60
h
5
{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx=2hf(x_{1})+{\frac {h^{3}}{3}}f''(x_{1})+{\frac {f^{(4)}(\xi _{1})}{60}}h^{5}}
وبتعويض قيمة
f
″
(
x
1
)
{\displaystyle f''(x_{1})}
∫
x
∘
x
2
f
(
x
)
d
x
=
2
h
f
(
x
1
)
+
h
3
3
{
1
h
2
[
f
(
x
∘
)
−
2
f
(
x
1
)
+
f
(
x
2
)
]
−
h
2
12
f
(
4
)
(
ξ
2
)
}
+
f
(
4
)
(
ξ
1
)
60
h
5
{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx=2hf(x_{1})+{\frac {h^{3}}{3}}\left\{{\frac {1}{h^{2}}}[f(x_{\circ })-2f(x_{1})+f(x_{2})]-{\frac {h^{2}}{12}}f^{(4)}(\xi _{2})\right\}+{\frac {f^{(4)}(\xi _{1})}{60}}h^{5}}
إذاً
∫
x
∘
x
2
f
(
x
)
d
x
=
h
3
[
f
(
x
∘
)
+
4
f
(
x
1
)
+
f
(
x
2
)
]
−
h
5
12
[
1
3
f
(
4
)
(
ξ
2
)
−
1
5
f
(
4
)
(
ξ
1
)
]
{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx={\frac {h}{3}}[f(x_{\circ })+4f(x_{1})+f(x_{2})]-{\frac {h^{5}}{12}}\left[{\frac {1}{3}}f^{(4)}(\xi _{2})-{\frac {1}{5}}f^{(4)}(\xi _{1})\right]}
[ 1]
صيغة سمبسون البسيطة :
الفكرة التي قامت عليها صيغة سمبسون البسيطة:
تقوم على أساس تقريب المساحة المطلوبة تحت منحنى الدالة
f
(
x
)
{\displaystyle f(x)}
P
2
(
x
)
{\displaystyle P_{2}(x)}
(
x
∘
,
f
(
x
∘
)
,
(
x
1
,
f
(
x
1
)
,
(
x
2
,
f
(
x
2
)
{\displaystyle (x_{\circ },f(x_{\circ }),(x_{1},f(x_{1}),(x_{2},f(x_{2})}
=
h
3
[
f
(
x
0
)
+
4
f
(
x
1
)
+
f
(
x
2
)
]
{\displaystyle ={\frac {h}{3}}[f(x_{0})+4f(x_{1})+f(x_{2})]}
∫
x
0
x
2
f
(
x
)
d
x
{\displaystyle \,\int _{x_{0}}^{x_{2}}f(x)\,dx}
بحيث تكون المساحة محصورة بين الخطين [ 1]
a
=
x
∘
,
b
=
x
2
{\displaystyle a=x_{\circ },b=x_{2}}
استخدام قاعدة سمبسون البسيطة لإيجاد قيمة التكامل
∫
0
2
x
2
d
x
{\displaystyle \,\int _{0}^{2}x^{2}\,dx}
الحل:
∫
x
0
x
2
f
(
x
)
d
x
{\displaystyle \,\int _{x_{0}}^{x_{2}}f(x)\,dx}
=
h
3
[
f
(
x
0
)
+
4
f
(
x
1
)
+
f
(
x
2
)
]
{\displaystyle ={\frac {h}{3}}[f(x_{0})+4f(x_{1})+f(x_{2})]}
x
0
=
0
{\displaystyle x_{0}=0}
x
1
=
1
{\displaystyle x_{1}=1}
x
2
=
2
{\displaystyle x_{2}=2}
h
=
x
1
−
x
0
=
1
−
0
=
1
{\displaystyle h=x_{1}-x_{0}=1-0=1}
→
∫
0
2
x
2
d
x
=
1
3
[
f
(
0
)
+
4
f
(
1
)
+
f
(
2
)
]
{\displaystyle \to \,\int _{0}^{2}x^{2}\,dx={\frac {1}{3}}[f(0)+4f(1)+f(2)]}
=
1
3
[
(
0
)
2
+
4
(
1
)
2
+
(
2
)
2
]
{\displaystyle ={\frac {1}{3}}[(0)^{2}+4(1)^{2}+(2)^{2}]}
=
1
3
[
8
]
=
8
3
{\displaystyle ={\frac {1}{3}}[8]={\frac {8}{3}}}
[ 1]
قد لا يكون استخدام صيغ نيوتن-كوتس عمليا في الكثير من الحالات خصوصا إذا كانت فترة التكامل كبيرة نسبيا ؛ حيث أننا نضطر ففي مثل هذه الحالات إلى استخدام كثيرات حدود ذات درجات عالية من الذبذبة وهذا بدوره يترك أثر سيئا على دقة الحلول العددية، وللتغلب على هذه المشكلة وتقليل الخطأ الناتج عن تطبيق صيغة تكامل ذات رتبة منخفضة أي قيمة n صغيرة فإننا نقسم فترة التكامل [a,b] إلى فترات أصغر ونطبق صيغة التكامل على كل فترة جزئية على حدة . والصيغة الناتجة من التطبيق المتكرر لصيغة ذات رتبة منخفضة تسمى صيغة تكامل مركبة . وفيما يلي نحصل على الصيغة المركبة الناتجة من التطبيق المتكرر لصيغة سمبسون البسيطة عندما (n=2) . ولكن من شروط قاعدة سمبسون ان تكون n عدد زوجي
وكي نكرر تطبيق قاعدة سمبسون عدد r من المرات نحتاج معرفة قيم الدالة عند عدد 2r+1 من نقاط الأساس
x
0
,
x
1
,
x
2
.
.
.
.
.
x
2
r
{\displaystyle x_{0},x_{1},x_{2}.....x_{2}r}
[ 2]
وتكون قاعدة سمبسون المركبة كالتالي :
∫
a
b
f
(
x
)
d
x
≈
h
3
[
f
(
x
0
)
+
2
∑
j
=
1
n
/
2
−
1
f
(
x
2
j
)
+
4
∑
j
=
1
n
/
2
f
(
x
2
j
−
1
)
+
f
(
x
n
)
]
{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\tfrac {h}{3}}{\bigg [}f(x_{0})+2\sum _{j=1}^{n/2-1}f(x_{2j})+4\sum _{j=1}^{n/2}f(x_{2j-1})+f(x_{n}){\bigg ]}}
∫
−
0.5
0.5
x
l
n
(
x
+
1
)
d
x
{\displaystyle \int _{-0.5}^{0.5}xln(x+1)dx\,}
n
=
6
{\displaystyle n=6}
الحل بستخدام قاعدة سمبسون المركبة :
x
0
=
−
0.5
{\displaystyle x_{0}=-0.5}
x
1
=
−
0.333
{\displaystyle x_{1}=-0.333}
x
2
=
−
0.166
{\displaystyle x_{2}=-0.166}
x
3
=
0
{\displaystyle x_{3}=0}
x
4
=
0.166
{\displaystyle x_{4}=0.166}
x
5
=
0.333
{\displaystyle x_{5}=0.333}
x
6
=
0.5
{\displaystyle x_{6}=0.5}
h
=
0.5
−
(
−
0.5
)
6
=
1
6
{\displaystyle h={\frac {0.5-(-0.5)}{6}}={\frac {1}{6}}}
∫
−
0.5
0.5
x
l
n
(
x
+
1
)
d
x
=
1
6
3
[
f
(
−
0.5
)
+
4
f
(
−
0.333
)
+
2
f
(
−
0.166
)
+
4
f
(
0
)
+
2
f
(
0.166
)
+
4
f
(
0.333
)
+
f
(
0.5
)
]
{\displaystyle \int _{-0.5}^{0.5}xln(x+1)dx\,={\frac {\frac {1}{6}}{3}}[f(-0.5)+4f(-0.333)+2f(-0.166)+4f(0)+2f(0.166)+4f(0.333)+f(0.5)]}
f
(
−
0.5
)
=
(
−
0.5
)
l
n
(
−
0.5
+
1
)
=
0.346
{\displaystyle f(-0.5)=(-0.5)ln(-0.5+1)=0.346}
f
(
−
0.333
)
=
(
−
0.333
)
l
n
(
−
0.333
+
1
)
=
0.134
{\displaystyle f(-0.333)=(-0.333)ln(-0.333+1)=0.134}
f
(
−
0.166
)
=
(
−
0.166
)
l
n
(
−
0.166
+
1
)
=
0.030
{\displaystyle f(-0.166)=(-0.166)ln(-0.166+1)=0.030}
f
(
0
)
=
(
0
)
l
n
(
0
+
1
)
=
0
{\displaystyle f(0)=(0)ln(0+1)=0}
f
(
0.166
)
=
(
0.166
)
l
n
(
0.166
+
1
)
=
0.025
{\displaystyle f(0.166)=(0.166)ln(0.166+1)=0.025}
f
(
0.333
)
=
(
0.333
)
l
n
(
0.333
+
1
)
=
0.095
{\displaystyle f(0.333)=(0.333)ln(0.333+1)=0.095}
f
(
0.5
)
=
(
0.5
)
l
n
(
0.5
+
1
)
=
0.202
{\displaystyle f(0.5)=(0.5)ln(0.5+1)=0.202}
∫
−
0.5
0.5
x
l
n
(
x
+
1
)
d
x
{\displaystyle \int _{-0.5}^{0.5}xln(x+1)dx\,}
=
1
18
[
0.346
+
0.536
+
0.06
+
0
+
0.05
+
0.38
+
0.202
]
{\displaystyle ={\frac {1}{18}}[0.346+0.536+0.06+0+0.05+0.38+0.202]}
=
1.574
18
=
0.0874
{\displaystyle ={\frac {1.574}{18}}=0.0874}
[ 1]
وفي بعض التكاملات على الفترات الكبيرة كان الحل بصيغة سمبسون يحتوي نسبة خطأ أكبر من لو قام بتقسيم الفترة
نلاحظ أن الخطأ يقل بعد تقسيم الفترة ومن هذا المنطلق نتجت صيغة سمبسون المركبة ولها العديد من الصيغ
∫
a
b
f
(
x
)
d
x
=
h
3
[
f
(
x
∘
)
+
2
∑
j
=
1
(
n
2
)
−
1
f
(
x
2
j
)
+
4
∑
j
=
1
n
2
f
(
x
2
j
−
1
)
+
f
(
x
n
)
]
−
h
5
90
∑
j
=
1
n
2
f
(
4
)
(
ξ
j
)
{\displaystyle \int _{a}^{b}f(x)dx={\frac {h}{3}}\left[f(x_{\circ })+2\sum _{j=1}^{({\frac {n}{2}})-1}f(x_{2_{j}})+4\sum _{j=1}^{\frac {n}{2}}f(x_{2_{j-1}})+f(x_{n})\right]-{\frac {h^{5}}{90}}\sum _{j=1}^{\frac {n}{2}}f^{(4)}(\xi _{j})}
h
=
b
−
a
n
{\displaystyle h={\frac {b-a}{n}}}
والخطأ مرتبط بالتقريب :
E
(
f
)
=
−
h
5
90
∑
j
=
1
n
2
f
(
4
)
(
ξ
j
)
{\displaystyle E(f)=-{\frac {h^{5}}{90}}\sum _{j=1}^{\frac {n}{2}}f^{(4)}(\xi _{j})}
[ 1]
حسب نظرية القيمة المتوسطة يوجد
μ
∈
(
a
,
b
)
{\displaystyle \mu \in (a,b)}
f
(
4
)
(
μ
)
=
n
2
∑
j
=
1
n
2
f
(
4
)
(
ξ
j
)
{\displaystyle f^{(4)}(\mu )={\frac {n}{2}}\sum _{j=1}^{\frac {n}{2}}f^{(4)}(\xi _{j})}
بالتالي
E
(
f
)
=
−
h
5
90
∑
j
=
1
n
2
f
(
4
)
(
ξ
j
)
=
−
h
5
180
n
f
(
4
)
(
μ
)
{\displaystyle E(f)=-{\frac {h^{5}}{90}}\sum _{j=1}^{\frac {n}{2}}f^{(4)}(\xi _{j})=-{\frac {h^{5}}{180}}nf^{(4)}(\mu )}
وبما أن
h
=
b
−
a
n
{\displaystyle h={\frac {b-a}{n}}}
و أيضاً
E
(
f
)
=
−
b
−
a
180
h
4
n
f
(
4
)
(
μ
)
{\displaystyle E(f)=-{\frac {b-a}{180}}h^{4}nf^{(4)}(\mu )}
^ ا ب ج د ه و “Numerical Analysis”, Richard L. Burden & J. Douglas Faires, Brooks/Cole, Cengage Learning,(2011) ISBN 0-534-38216-9
^ كتاب الطرق العددية والتحليل العددي - أ.د. أبو بكر أحمد السيد، جامعة الكويت ISBN9789957171353
كتاب الطرق العددية والتحليل العددي - أ.د. أبو بكر أحمد السيد، جامعة الكويت ISBN9789957171353
“Numerical Analysis”, Richard L. Burden & J. Douglas Faires, Brooks/Cole, Cengage Learning,(2011) ISBN 0-534-38216-9
![{\displaystyle \int _{a}^{b}f(x)dx={\int _{x_{\circ }}^{x_{2}}}[{\frac {(x-x_{1})(x-x_{2})}{(x_{\circ }-x_{1})(x_{\circ }-x_{2})}}f(x_{\circ })+{\frac {(x-x_{\circ })(x-x_{2})}{(x_{1}-x_{\circ })(x_{1}-x_{2})}}f(x_{1})+{\frac {(x-x_{\circ })(x-x_{1})}{(x_{2}-x_{\circ })(x_{2}-x_{1})}}f(x_{2})]dx+{\int _{x_{\circ }}^{x_{2}}}{\frac {(x-x_{\circ })(x-x_{1})(x-x_{2})}{6}}f^{(3)}(\xi (x))dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/97c22dfe6e763295afa2c6c567b47dec9f5da391)
![{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx-\left[f(x_{1})(x-x_{1})+{\frac {f'(x_{1})}{2}}(x-x_{1})^{2}+{\frac {f''(x_{1})}{6}}(x-x_{1})^{3}+{\frac {f'''(\xi (x)}{24}}(x-x_{1})^{4}\right]_{x_{\circ }}^{x_{2}}+{\frac {1}{24}}\int _{x_{\circ }}^{x_{2}}f^{(4)}(\xi (x))(x-x_{1})^{4}dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/44530110e4878469b99cbcbaa094bb3bb0af982d)
![{\displaystyle {\frac {1}{24}}\int _{x_{\circ }}^{x_{2}}f^{(4)}(\xi (x))(x-x_{1})^{4}dx={\frac {f^{(4)}(\xi _{1})}{24}}\int _{x_{\circ }}^{x_{2}}(x-x_{1})^{4}dx=\left[{\frac {f^{(4)}(\xi _{1})}{120}}(x-x_{1})^{5}\right]_{x_{\circ }}^{x_{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc483a4d273af909dc953e09bd50564fc79f1de6)
![{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx=2hf(x_{1})+{\frac {h^{3}}{3}}\left\{{\frac {1}{h^{2}}}[f(x_{\circ })-2f(x_{1})+f(x_{2})]-{\frac {h^{2}}{12}}f^{(4)}(\xi _{2})\right\}+{\frac {f^{(4)}(\xi _{1})}{60}}h^{5}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf7f16e09b06b562c77106f6bfac74ab024a0e37)
![{\displaystyle \int _{x_{\circ }}^{x_{2}}f(x)dx={\frac {h}{3}}[f(x_{\circ })+4f(x_{1})+f(x_{2})]-{\frac {h^{5}}{12}}\left[{\frac {1}{3}}f^{(4)}(\xi _{2})-{\frac {1}{5}}f^{(4)}(\xi _{1})\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0076db5325ed8f2bbeacb627b77df969db76947)
![{\displaystyle ={\frac {h}{3}}[f(x_{0})+4f(x_{1})+f(x_{2})]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a8619d0fc05f082c63d44a86dd656bf73bb8aad)
![{\displaystyle \to \,\int _{0}^{2}x^{2}\,dx={\frac {1}{3}}[f(0)+4f(1)+f(2)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/85d06629237c7f205769741bd0f2186b7f2c9ac9)
![{\displaystyle ={\frac {1}{3}}[(0)^{2}+4(1)^{2}+(2)^{2}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e03c3ecbc32e15345bc5f9bd7866ba58f0077890)
![{\displaystyle ={\frac {1}{3}}[8]={\frac {8}{3}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc8747db3347773eceef9bb54550e476575ab557)
![{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\tfrac {h}{3}}{\bigg [}f(x_{0})+2\sum _{j=1}^{n/2-1}f(x_{2j})+4\sum _{j=1}^{n/2}f(x_{2j-1})+f(x_{n}){\bigg ]}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53796a1aa92d5950c5378f5e30672fa0929a2b0d)
![{\displaystyle \int _{-0.5}^{0.5}xln(x+1)dx\,={\frac {\frac {1}{6}}{3}}[f(-0.5)+4f(-0.333)+2f(-0.166)+4f(0)+2f(0.166)+4f(0.333)+f(0.5)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e7da15984c9652c93fbe6e640f004a60c937ee94)
![{\displaystyle ={\frac {1}{18}}[0.346+0.536+0.06+0+0.05+0.38+0.202]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c5ac2b82a1fee3d4acff83551cccfd2482321785)
![{\displaystyle \int _{a}^{b}f(x)dx={\frac {h}{3}}\left[f(x_{\circ })+2\sum _{j=1}^{({\frac {n}{2}})-1}f(x_{2_{j}})+4\sum _{j=1}^{\frac {n}{2}}f(x_{2_{j-1}})+f(x_{n})\right]-{\frac {h^{5}}{90}}\sum _{j=1}^{\frac {n}{2}}f^{(4)}(\xi _{j})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0ca6cc6319a1eea67833f544d3087fd66e9a2b9f)
