MCP рдХреЗ рдорд╛рдзреНрдпрдо рд╕реЗ рдХрдиреЗрдХреНрдЯ рдХрд░реЗрдВ тЖТ

рдЧрдгрдирд╛ рджрд░реНрдЬ рдХрд░реЗрдВ

рд╕реВрддреНрд░ (рдлреЙрд░реНрдореВрд▓рд╛)

рд╡рд┐рдЬреНрдЮрд╛рдкрди

рдкрд░рд┐рдгрд╛рдо

Simpson's rule result (n = 64)
3.1415926535892149
рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓рди рдХрд╛ рдЕрдиреБрдорд╛рди
рдЯреНрд░реИрдкреЗрдЬрд╝реЙрдЗрдбрд▓ рдирд┐рдпрдо
3.1415519634856541
рдорд┐рдбрдкреЙрдЗрдВрдЯ рдирд┐рдпрдо
3.1416129986418473
рд╕рд┐рдореНрдкрд╕рди рдирд┐рдпрдо
3.1415926535892149
n рдЯреНрд░реИрдкреЗрдЬрд╝реЙрдЗрдбрд▓ рдорд┐рдбрдкреЙрдЗрдВрдЯ рд╕рд┐рдореНрдкрд╕рди
2 3.1000000000000001 3.1623529411764704 3.1333333333333333
4 3.1311764705882354 3.1468005183939427 3.1415686274509804
8 3.1389884944910889 3.1428947295916889 3.1415925024587068
16 3.1409416120413889 3.1419181743085600 3.1415926512248218
32 3.1414298931749749 3.1416740337963374 3.1415926535528360
64 3.1415519634856541 3.1416129986418473 3.1415926535892149

рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ

рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рдлрд╝рдВрдХреНрд╢рди \(f(x)\) рдХреЗ рдПрдХ рд╕реАрдорд┐рдд рдЕрдВрддрд░рд╛рд▓ \([a, b]\) рдкрд░ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓рди (definite integral) рдХрд╛ рдЕрдиреБрдорд╛рди рддреАрди рдХреНрд▓рд╛рд╕рд┐рдХ рд╕рдВрдпреБрдХреНрдд рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдирд┐рдпрдореЛрдВ рд╕реЗ рд▓рдЧрд╛рддрд╛ рд╣реИ: рдЯреНрд░реИрдкреЗрдЬрд╝реЙрдЗрдбрд▓ рдирд┐рдпрдо, рдорд┐рдбрдкреЙрдЗрдВрдЯ рдирд┐рдпрдо рдФрд░ рд╕рд┐рдореНрдкрд╕рди рдирд┐рдпрдоред рдпрд╣ рдХреЗрд╡рд▓ рдПрдХ рд╣реА рд╕рдВрдЦреНрдпрд╛ рдирд╣реАрдВ рджреЗрддрд╛, рдмрд▓реНрдХрд┐ рд╣рд░ рдирд┐рдпрдо рдХреЛ 2, 4, 8, 16, ... рдЙрдкрдЦрдВрдбреЛрдВ рдкрд░ тАФ рдпрд╛рдиреА рд╣рд░ рдмрд╛рд░ рджреЛрдЧреБрдирд╛ рдХрд░рддреЗ рд╣реБрдП рдЖрдкрдХреА рдЪреБрдиреА рд╣реБрдИ рдЕрдзрд┐рдХрддрдо рд╕рдВрдЦреНрдпрд╛ \(N\) рддрдХ тАФ рд╣рд▓ рдХрд░рддрд╛ рд╣реИ рдФрд░ рдПрдХ рдЕрднрд┐рд╕рд░рдг рддрд╛рд▓рд┐рдХрд╛ рджрд┐рдЦрд╛рддрд╛ рд╣реИред рдЗрд╕рд╕реЗ рдЖрдк рджреЗрдЦ рд╕рдХрддреЗ рд╣реИрдВ рдХрд┐ рдЕрдиреБрдорд╛рди рдХреИрд╕реЗ рд╕реНрдерд┐рд░ рд╣реЛрддреЗ рдЬрд╛рддреЗ рд╣реИрдВ рдФрд░ рд╕рдЯреАрдХрддрд╛ рдХрд╛ рдЕрдВрджрд╛рдЬрд╝рд╛ рдЦреБрдж рд▓рдЧрд╛ рд╕рдХрддреЗ рд╣реИрдВред

a рдФрд░ b рдХреЗ рдмреАрдЪ рдКрд░реНрдзреНрд╡рд╛рдзрд░ рдкрдЯреНрдЯрд┐рдпреЛрдВ рд╕реЗ рдЕрдиреБрдорд╛рдирд┐рдд рд╡рдХреНрд░ рдХреЗ рдиреАрдЪреЗ рдХрд╛ рдХреНрд╖реЗрддреНрд░рдлрд▓
рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕рдорд╛рдХрд▓рди рд╕рд░рд▓ рдкрдЯреНрдЯрд┐рдпреЛрдВ рдХреЛ рдЬреЛрдбрд╝рдХрд░ a рд╕реЗ b рддрдХ f(x) рдХреЗ рдиреАрдЪреЗ рдХрд╛ рдХреНрд╖реЗрддреНрд░рдлрд▓ рдЕрдиреБрдорд╛рдирд┐рдд рдХрд░рддрд╛ рд╣реИред

рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ

рд╕рдорд╛рдХрд▓реНрдп (integrand) рдХреЛ рдЪрд░ x рдореЗрдВ рдПрдХ рдЧрдгрд┐рддреАрдп рдПрдХреНрд╕рдкреНрд░реЗрд╢рди рдХреЗ рд░реВрдк рдореЗрдВ рд▓рд┐рдЦреЗрдВ (рдЬреИрд╕реЗ 4/(1+x^2) рдпрд╛ sin(x)*exp(-x))ред рдпрд╣рд╛рдБ + - * / ^ рдСрдкрд░реЗрдЯрд░ рдФрд░ рдХреЛрд╖реНрдардХ рд╕рдорд░реНрдерд┐рдд рд╣реИрдВ, рд╕рд╛рде рд╣реА sin, cos, tan, exp, log/ln, sqrt рдФрд░ abs рдЬреИрд╕реЗ рд╕рд╛рдорд╛рдиреНрдп рдлрд╝рдВрдХреНрд╢рди рддрдерд╛ pi рдФрд░ e рдЬреИрд╕реЗ рд╕реНрдерд┐рд░рд╛рдВрдХ рднреАред рдирд┐рдЪрд▓реА рд╕реАрдорд╛ \(a\) рдФрд░ рдКрдкрд░реА рд╕реАрдорд╛ \(b\) рднрд░реЗрдВ, рдЙрдкрдЦрдВрдбреЛрдВ рдХреА рдЕрдзрд┐рдХрддрдо рд╕рдВрдЦреНрдпрд╛ \(N\) рдЪреБрдиреЗрдВ (2 рдХреА рдШрд╛рдд рдореЗрдВ), рдФрд░ рддрдп рдХрд░реЗрдВ рдХрд┐ рдХрд┐рддрдиреЗ рдЕрдВрдХ рддрдХ рдкрд░рд┐рдгрд╛рдо рджрд┐рдЦрд╛рдирд╛ рд╣реИред рдореБрдЦреНрдп рдкрд░рд┐рдгрд╛рдо рдХреЗ рд░реВрдк рдореЗрдВ \(n = N\) рдкрд░ рд╕рд┐рдореНрдкрд╕рди рдХрд╛ рдЕрдиреБрдорд╛рди рджрд┐рдЦрд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдХреНрдпреЛрдВрдХрд┐ рдпрд╣ рдЖрдорддреМрд░ рдкрд░ рд╕рдмрд╕реЗ рддреЗрдЬрд╝реА рд╕реЗ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛрддрд╛ рд╣реИред

рд╕реВрддреНрд░реЛрдВ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛

\(n\) рдЙрдк-рдЕрдВрддрд░рд╛рд▓реЛрдВ рдХреЗ рд▓рд┐рдП рд╕реНрдЯреЗрдк \(h = (b - a)/n\) рд╣реЛрддрд╛ рд╣реИ рдФрд░ рдиреЛрдбреНрд╕ \(x_i = a + i\,h\) рд╣реЛрддреЗ рд╣реИрдВред рдЯреНрд░реИрдкреЗрдЬрд╝реЙрдЗрдбрд▓ рдирд┐рдпрдо рдлрд╝рдВрдХреНрд╢рди рдХреЗ рдорд╛рдиреЛрдВ рдХреЛ рдЬреЛрдбрд╝рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рджреЛрдиреЛрдВ рд╕рд┐рд░реЛрдВ (endpoints) рдХрд╛ рднрд╛рд░ рдЖрдзрд╛ рд░рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИред рдорд┐рдбрдкреЙрдЗрдВрдЯ рдирд┐рдпрдо рд╣рд░ рдЙрдк-рдЕрдВрддрд░рд╛рд▓ рдХреЗ рдмреАрдЪреЛрдВрдмреАрдЪ рдХрд╛ рдорд╛рди рд▓реЗрддрд╛ рд╣реИред рд╕рд┐рдореНрдкрд╕рди рдирд┐рдпрдо рдЗрди рджреЛрдиреЛрдВ рдХреЛ 1, 4, 2, 4, ..., 4, 1 рднрд╛рд░реЛрдВ рдХреЗ рд╕рд╛рде рдорд┐рд▓рд╛рддрд╛ рд╣реИ; рдпрд╣ рдШрди (cubic) рдлрд╝рдВрдХреНрд╢рдиреЛрдВ рдХреЗ рд▓рд┐рдП рдмрд┐рд▓рдХреБрд▓ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИ рдФрд░ рдЗрд╕рдХреА рддреНрд░реБрдЯрд┐ \(h^4\) рдХреНрд░рдо рдХреА рд╣реЛрддреА рд╣реИ, рдЬрдмрдХрд┐ рдмрд╛рдХреА рджреЛрдиреЛрдВ рдХреА рддреНрд░реБрдЯрд┐ \(h^2\) рдХреНрд░рдо рдХреАред

$$\int_{a}^{b} f(x)\,dx \approx \frac{h}{3}\left[ f(x_0) + 4\sum_{i\,\text{odd}} f(x_i) + 2\sum_{i\,\text{even}} f(x_i) + f(x_N) \right]$$

рд╡рдХреНрд░ рдХреЗ рдиреАрдЪреЗ рд╕рдорд▓рдВрдм, рдордзреНрдпрдмрд┐рдВрджреБ рдФрд░ рд╕рд┐рдореНрдкрд╕рди рдкрдЯреНрдЯреА рдЖрдХреГрддрд┐рдпреЛрдВ рдХреА рддреБрд▓рдирд╛ рдХрд░рддреЗ рддреАрди рдкреИрдирд▓
рддреАрдиреЛрдВ рдирд┐рдпрдо рдЗрд╕рдореЗрдВ рднрд┐рдиреНрди рд╣реИрдВ рдХрд┐ рд╣рд░ рдкрдЯреНрдЯреА рд╡рдХреНрд░ рдХрд╛ рдЕрдиреБрдорд╛рди рдХреИрд╕реЗ рд▓рдЧрд╛рддреА рд╣реИ: рд╕реАрдзреА рдЬреАрд╡рд╛, рдордзреНрдпрдмрд┐рдВрджреБ рдкрд░ рд╕рдкрд╛рдЯ рдЖрдпрдд, рдФрд░ рдлрд┐рдЯ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдкрд░рд╡рд▓рдпред

рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг

\(f(x) = 4/(1+x^2)\) рдХреЗ рд▓рд┐рдП \([0, 1]\) рдкрд░ рд╕рдЯреАрдХ рд╕рдорд╛рдХрд▓рди \(\pi = 3.14159265\ldots\) рд╣реИред \(n = 4\) рдФрд░ \(h = 0.25\) рдХреЗ рд╕рд╛рде рдЯреНрд░реИрдкреЗрдЬрд╝реЙрдЗрдбрд▓ рдЕрдиреБрдорд╛рди рд▓рдЧрднрдЧ \(3.131176\), рдорд┐рдбрдкреЙрдЗрдВрдЯ рдЕрдиреБрдорд╛рди рд▓рдЧрднрдЧ \(3.146801\) рдФрд░ рд╕рд┐рдореНрдкрд╕рди рдЕрдиреБрдорд╛рди рд▓рдЧрднрдЧ \(3.141569\) рдЖрддрд╛ рд╣реИ тАФ рдЬреЛ рдкрд╣рд▓реЗ рд╣реА рдкрд╛рдБрдЪ рджрд╢рдорд▓рд╡ рдЕрдВрдХреЛрдВ рддрдХ рд╕рдЯреАрдХ рд╣реИред рдЬреИрд╕реЗ-рдЬреИрд╕реЗ \(n\) рдмрдврд╝рддрд╛ рд╣реИ, рддреАрдиреЛрдВ рдЕрдиреБрдорд╛рди \(\pi\) рдХреЗ рдкрд╛рд╕ рдкрд╣реБрдБрдЪрддреЗ рдЬрд╛рддреЗ рд╣реИрдВред

рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рд╕рд╡рд╛рд▓

N рдХреЛ 2 рдХреА рдШрд╛рдд рд╣реА рдХреНрдпреЛрдВ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП? рдЙрдкрдЦрдВрдбреЛрдВ рдХреЛ рджреЛрдЧреБрдирд╛ рдХрд░рдиреЗ рд╕реЗ рд▓рдЧрд╛рддрд╛рд░ рдкрдВрдХреНрддрд┐рдпреЛрдВ рдХреА рддреБрд▓рдирд╛ рдЖрд╕рд╛рди рд╣реЛ рдЬрд╛рддреА рд╣реИ рдФрд░ рд╕рд┐рдореНрдкрд╕рди рдирд┐рдпрдо рдХреЗ рд▓рд┐рдП рдЬрд╝рд░реВрд░реА рд╕рдо (even) рд╕рдВрдЦреНрдпрд╛ рднреА рд╕реБрдирд┐рд╢реНрдЪрд┐рдд рд░рд╣рддреА рд╣реИред

рдХреМрди-рд╕реЗ рдлрд╝рдВрдХреНрд╢рди рдЗрд╕рдХреЗ рджрд╛рдпрд░реЗ рд╕реЗ рдмрд╛рд╣рд░ рд╣реИрдВ? рдЕрднрд┐рд╕рд░рдг рддрд╛рд▓рд┐рдХрд╛ рдорд╛рдирддреА рд╣реИ рдХрд┐ рд╕рдорд╛рдХрд▓реНрдп рдЪрд┐рдХрдирд╛ (analytic) рдФрд░ рдЕрдирд╛рд╡рд░реНрддреА (non-periodic) рд╣реИред рдЬрд┐рди рдлрд╝рдВрдХреНрд╢рдиреЛрдВ рдореЗрдВ \([a, b]\) рдХреЗ рднреАрддрд░ рд╡рд┐рдЪрд┐рддреНрд░рддрд╛ (singularity) рд╣реЛ тАФ рдЬреИрд╕реЗ рд╢реВрдиреНрдп рд╕реЗ рд╣реЛрдХрд░ рдЧреБрдЬрд╝рд░рддрд╛ \(1/x\) тАФ рд╡реЗ рдЕрдирдВрдд рдпрд╛ рдмреЗрдорддрд▓рдм рдкрд░рд┐рдгрд╛рдо рджреЗрдВрдЧреЗред

рдЕрдЧрд░ a рдФрд░ b рдмрд░рд╛рдмрд░ рд╣реЛрдВ, рдпрд╛ a, b рд╕реЗ рдмрдбрд╝рд╛ рд╣реЛ рддреЛ? рдЕрдЧрд░ \(a = b\) рд╣реИ рддреЛ рд╕рдорд╛рдХрд▓рди \(0\) рд╣реЛрдЧрд╛ред рдЕрдЧрд░ \(a > b\) рд╣реИ рддреЛ рдкрд░рд┐рдгрд╛рдо \([b, a]\) рдкрд░ рд▓рд┐рдП рдЧрдП рд╕рдорд╛рдХрд▓рди рдХрд╛ рдЪрд┐рд╣реНрдирд┐рдд (рдЛрдгрд╛рддреНрдордХ) рдорд╛рди рд╣реЛрдЧрд╛ред

рдЕрдВрддрд┐рдо рдЕрдкрдбреЗрдЯ: