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

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

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

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

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

рд╕рдорд╛рдХрд▓ рдХрд╛ рдорд╛рди
1.9999999545
рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (Tanh-Sinh) рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рджреНрд╡рд╛рд░рд╛ рдЕрдиреБрдорд╛рдирд┐рдд
рд╕рдорд╛рдХрд▓реНрдп f(x) 1/sqrt(x)
рдЕрдВрддрд░рд╛рд▓ [ 0 , 1 ]
рд▓рдХреНрд╖рд┐рдд рдЕрдВрдХ 15
рд╡рд┐рдзрд┐ DE / Tanh-Sinh рд╕рдорд▓рдВрдм рдирд┐рдпрдо

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

рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рдПрдХ-рдЪрд░ рдлрд▓рди \(f(x)\) рдХрд╛ рдХрд┐рд╕реА рдкрд░рд┐рдорд┐рдд рдЕрдВрддрд░рд╛рд▓ \([a, b]\) рдкрд░ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓рди рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (DE) рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рд╕реЗ рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИ, рдЬрд┐рд╕реЗ Tanh-Sinh рд╡рд┐рдзрд┐ рднреА рдХрд╣рд╛ рдЬрд╛рддрд╛ рд╣реИред рдкрд░рд┐рдорд┐рдд рдЕрдВрддрд░рд╛рд▓реЛрдВ рдХреЗ рд▓рд┐рдП DE рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рд╕рдмрд╕реЗ рднрд░реЛрд╕реЗрдордВрдж рд╕рд╛рдорд╛рдиреНрдп-рдЙрджреНрджреЗрд╢реНрдп рд╡рд┐рдзрд┐рдпреЛрдВ рдореЗрдВ рд╕реЗ рдПрдХ рд╣реИ рдФрд░ рдЙрди рд╕рдорд╛рдХрд▓реНрдпреЛрдВ рдХреЛ рд╕рдВрднрд╛рд▓рдиреЗ рдореЗрдВ рдЦрд╛рд╕ рддреМрд░ рдкрд░ рдорд╛рд╣рд┐рд░ рд╣реИ рдЬреЛ рдХрд┐рд╕реА рд╕рд┐рд░реЗ рдкрд░ рдЕрдирдВрдд рд╣реЛ рдЬрд╛рддреЗ рд╣реИрдВ тАФ рдЬреИрд╕реЗ \(1/\sqrt{x}\) рдпрд╛ \(\log(x)\) тАФ рдЬрд╣рд╛рдБ рд╕рд╛рдорд╛рдиреНрдп рдЧреЙрд╕рд┐рдпрди рдпрд╛ рд╕рд┐рдореНрдкрд╕рди рдирд┐рдпрдо рд▓рдбрд╝рдЦрдбрд╝рд╛ рдЬрд╛рддреЗ рд╣реИрдВред

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

рд╕рдорд╛рдХрд▓реНрдп f(x) рдмреЙрдХреНрд╕ рдореЗрдВ рдлрд▓рди рдХреЛ рд╕рд╛рдорд╛рдиреНрдп рдЧрдгрд┐рддреАрдп рд╕рдВрдХреЗрддрди рдореЗрдВ рд▓рд┐рдЦреЗрдВ: + - * / ^, рдХреЛрд╖реНрдардХ, рдФрд░ рдлрд▓рди sqrt, exp, log, ln, sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, abs рдХреЗ рд╕рд╛рде-рд╕рд╛рде рдЕрдЪрд░ pi рдФрд░ eред рдирд┐рдЪрд▓реА рд╕реАрдорд╛ \(a\) рдФрд░ рдКрдкрд░реА рд╕реАрдорд╛ \(b\) рджрд░реНрдЬ рдХрд░реЗрдВ, рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ рддрдХ рд╕рдЯреАрдХрддрд╛ рдЪрд╛рд╣рд┐рдП рд╡рд╣ рдЪреБрдиреЗрдВ, рдФрд░ рд╕рдмрдорд┐рдЯ рдХрд░реЗрдВред рд╡рд┐рдЪрд┐рддреНрд░рддрд╛рдПрдБ (singularities) рдХреЗрд╡рд▓ рд╕рд┐рд░реЛрдВ \(a\) рдФрд░ \(b\) рдкрд░ рд╣реА рд╕реНрд╡реАрдХрд╛рд░реНрдп рд╣реИрдВ; рдЗрд╕рдХреЗ рдЕрд▓рд╛рд╡рд╛ рдлрд▓рди рдЦреБрд▓реЗ рдЕрдВрддрд░рд╛рд▓ \((a, b)\) рдкрд░ рд╡рд┐рд╢реНрд▓реЗрд╖рдХ (analytic) рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП рдФрд░ рдЖрд╡рд░реНрддреА (periodic) рдирд╣реАрдВ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдПред

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

рд╕рдмрд╕реЗ рдкрд╣рд▓реЗ рдЕрдВрддрд░рд╛рд▓ рдХреЛ \(x(u) = \frac{(b+a)+(b-a)u}{2}\) рдХреЗ рджреНрд╡рд╛рд░рд╛ \([-1, 1]\) рдкрд░ рдореИрдк рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред рдлрд┐рд░ DE рд░реВрдкрд╛рдВрддрд░рдг \(u = \tanh\!\left(\frac{\pi}{2}\sinh t\right)\) рд░реЗрдЦрд╛ рдХреЛ рдЗрд╕ рддрд░рд╣ рдЦреАрдВрдЪрддрд╛ рд╣реИ рдХрд┐ рдЬреИрд╕реЗ-рдЬреИрд╕реЗ \(t\) рдмрдврд╝рддрд╛ рд╣реИ, \(u\) рд╕рд┐рд░реЛрдВ рдХреА рдУрд░ рдЕрддрд┐-рдШрд╛рддрд╛рдВрдХреАрдп (super-exponentially) рдЧрддрд┐ рд╕реЗ рдкрд╣реБрдБрдЪрддрд╛ рд╣реИ, рдЬрдмрдХрд┐ рднрд╛рд░ \(g'(t)\) рдЙрддрдиреА рд╣реА рддреЗрдЬрд╝реА рд╕реЗ рд╢реВрдиреНрдп рдХреА рдУрд░ рдЧрд┐рд░ рдЬрд╛рддрд╛ рд╣реИред рдЪреВрдБрдХрд┐ рдиреЛрдбреНрд╕ рдХрднреА рдареАрдХ \(a\) рдпрд╛ \(b\) рдкрд░ рдирд╣реАрдВ рдкрдбрд╝рддреЗ, рдЗрд╕рд▓рд┐рдП рд╕рд┐рд░реЗ рдХреА рд╡рд┐рдЪрд┐рддреНрд░рддрд╛ рдХрд╛ рдХрднреА рд╕реАрдзреЗ рдореВрд▓реНрдпрд╛рдВрдХрди рд╣реА рдирд╣реАрдВ рд╣реЛрддрд╛ тАФ рд╡рд╣ "рд╕рд╛рдз" рд▓реА рдЬрд╛рддреА рд╣реИред рдлрд┐рд░ рд░реВрдкрд╛рдВрддрд░рд┐рдд рд╕рдорд╛рдХрд▓ рдХреЛ рдЪрд░рдг \(h\) рдХреЗ рд╕рд╛рдзрд╛рд░рдг рд╕рдорд▓рдВрдм (trapezoidal) рдирд┐рдпрдо рд╕реЗ рдЬреЛрдбрд╝рд╛ рдЬрд╛рддрд╛ рд╣реИ, рдФрд░ \(h\) рдХреЛ рддрдм рддрдХ рдЖрдзрд╛ рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ рдЬрдм рддрдХ рдЙрддреНрддрд░ рдмрджрд▓рдирд╛ рдмрдВрдж рди рдХрд░ рджреЗред

$$\int_{a}^{b} f(x)\,dx \;\approx\; \frac{b-a}{2}\,h\sum_{k} w_k\, f\!\left(x_k\right)$$

$$\text{where}\quad \left\{ \begin{aligned} a &= \text{Lower limit} \\ b &= \text{Upper limit} \\ x_k &= \tfrac{b+a}{2} + \tfrac{b-a}{2}\tanh\!\left(\tfrac{\pi}{2}\sinh(k h)\right) \\ w_k &= \dfrac{\tfrac{\pi}{2}\cosh(k h)}{\cosh^{2}\!\left(\tfrac{\pi}{2}\sinh(k h)\right)} \end{aligned} \right.$$

t рдореЗрдВ рд╕рдорд╛рди рджреВрд░реА рд╡рд╛рд▓реЗ рдиреЛрдбреНрд╕, рдЬреЛ рд╕рд┐рд░реЛрдВ a рдФрд░ b рдХреЗ рдкрд╛рд╕ рдЧреБрдЪреНрдЫрд┐рдд рдиреЛрдбреНрд╕ рдореЗрдВ рдмрджрд▓ рдЬрд╛рддреЗ рд╣реИрдВ
рджреНрд╡рд┐-рдШрд╛рддрд╛рдВрдХреАрдп рдЪрд░ рдкрд░рд┐рд╡рд░реНрддрди рд╕рдорд╛рди рдиреЛрдбреНрд╕ рдХреЛ рдЙрди рдмрд┐рдВрджреБрдУрдВ рдореЗрдВ рдмрджрд▓ рджреЗрддрд╛ рд╣реИ рдЬреЛ рд╕рд┐рд░реЛрдВ a рдФрд░ b рдХреЗ рдкрд╛рд╕ рдЗрдХрдЯреНрдард╛ рд╣реЛрддреЗ рд╣реИрдВред

рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг

\(f(x) = 1/\sqrt{x}\) рдХреЗ рд▓рд┐рдП \([0, 1]\) рдкрд░ рд╕рдЯреАрдХ рдорд╛рди \([2\sqrt{x}]\) рдХрд╛ \(0\) рд╕реЗ \(1\) рддрдХ \(= 2\) рд╣реЛрддрд╛ рд╣реИред рдПрдХ рдореЛрдЯрд╛ 7-рдмрд┐рдВрджреБ DE рдЧреНрд░рд┐рдб (\(h = 0.5\)) рд╣реА рд▓рдЧрднрдЧ \(1.94\) рджреЗ рджреЗрддрд╛ рд╣реИ; \(h\) рдХреЛ рдФрд░ рдорд╣реАрди рдХрд░рдиреЗ рдкрд░ рдЕрдиреБрдорд╛рди \(2.000000000000000\) рддрдХ рдкрд╣реБрдБрдЪ рдЬрд╛рддрд╛ рд╣реИред рдПрдХ рдЧреИрд░-рд╡рд┐рдЪрд┐рддреНрд░ рдЬрд╛рдБрдЪ, \(f(x) = x^2\) рдкрд░ \([0, 1]\), рдкрд░рд┐рдгрд╛рдо \(1/3 = 0.3333333333333\) рджреЗрддреА рд╣реИред

a рдФрд░ b рдХреЗ рдмреАрдЪ рд╡рдХреНрд░ f(x) рдХреЗ рдиреАрдЪреЗ рдЫрд╛рдпрд╛рдВрдХрд┐рдд рдХреНрд╖реЗрддреНрд░рдлрд▓, рдЬрд┐рд╕рдореЗрдВ рд╕рд┐рд░реЛрдВ рдХреЗ рдкрд╛рд╕ рд╕рдШрди рдкреНрд░рддрд┐рдЪрдпрди рдмрд┐рдВрджреБ рд╣реИрдВ
рд╕рдорд╛рдХрд▓ рдЫрд╛рдпрд╛рдВрдХрд┐рдд рдХреНрд╖реЗрддреНрд░рдлрд▓ рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯрди рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рд╡рд┐рдЪрд┐рддреНрд░рддрд╛рдУрдВ рдХреЛ рд╕рдВрднрд╛рд▓рдиреЗ рдХреЗ рд▓рд┐рдП рдиреЛрдбреНрд╕ рд╕рд┐рд░реЛрдВ рдХреЗ рдкрд╛рд╕ рд╕рдШрди рд╣реЛрддреЗ рд╣реИрдВред

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

рдХреНрдпрд╛ рдпрд╣ рдЕрдВрддрд░рд╛рд▓ рдХреЗ рднреАрддрд░ рдХреА рд╡рд┐рдЪрд┐рддреНрд░рддрд╛ рд╕рдВрднрд╛рд▓ рд╕рдХрддрд╛ рд╣реИ? рдирд╣реАрдВ тАФ DE рдХреЗрд╡рд▓ рд╕рд┐рд░реЛрдВ рдкрд░ рд╣реА рд╡рд┐рдЪрд┐рддреНрд░рддрд╛рдПрдБ рд╕рд╣рди рдХрд░рддрд╛ рд╣реИред рдпрджрд┐ рднреАрддрд░ рдХрд┐рд╕реА рдмрд┐рдВрджреБ \(c\) рдкрд░ рд╡рд┐рдЪрд┐рддреНрд░рддрд╛ рд╣реЛ, рддреЛ рд╕рдорд╛рдХрд▓ рдХреЛ \([a, c]\) рдФрд░ \([c, b]\) рдореЗрдВ рдмрд╛рдБрдЯ рд▓реЗрдВ рдФрд░ рджреЛрдиреЛрдВ рдкрд░рд┐рдгрд╛рдореЛрдВ рдХреЛ рдЬреЛрдбрд╝ рджреЗрдВред

рдЖрд╡рд░реНрддреА (periodic) рдлрд▓рдиреЛрдВ рдХреЗ рд▓рд┐рдП рдпрд╣ рдХреНрдпреЛрдВ рдЦрд░рд╛рдм рд╣реИ? рдЪрд┐рдХрдиреЗ рдЖрд╡рд░реНрддреА рд╕рдорд╛рдХрд▓реНрдпреЛрдВ рдХреЗ рд▓рд┐рдП рд╕рд╛рдзрд╛рд░рдг рд╕рдорд▓рдВрдм рдирд┐рдпрдо рдкрд╣рд▓реЗ рд╕реЗ рд╣реА рдШрд╛рддрд╛рдВрдХреАрдп рд░реВрдк рд╕реЗ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП DE рдЪрд░-рдкрд░рд┐рд╡рд░реНрддрди рдЪреАрдЬрд╝реЛрдВ рдХреЛ рдХреЗрд╡рд▓ рдзреАрдорд╛ рд╣реА рдХрд░рддрд╛ рд╣реИред

"рдЕрдВрдХ" рд╡рд╛рд▓реА рд╕реЗрдЯрд┐рдВрдЧ рдХреНрдпрд╛ рдХрд░рддреА рд╣реИ? рдпрд╣ рд╡рд╣ рд╕рд╛рдкреЗрдХреНрд╖ рд╕рд╣рдирд╢реАрд▓рддрд╛ (relative tolerance) рддрдп рдХрд░рддреА рд╣реИ рдЬреЛ рдирд┐рд░реНрдгрдп рдХрд░рддреА рд╣реИ рдХрд┐ рдорд╣реАрдирдХрд░рдг рдХрдм рд░реБрдХреЗ, рдФрд░ рдкреНрд░рджрд░реНрд╢рд┐рдд рдорд╛рди рдХреЛ рдЙрд╕реА рдХреЗ рдЕрдиреБрд╕рд╛рд░ рдкреВрд░реНрдгрд╛рдВрдХрд┐рдд рдХрд░рддреА рд╣реИред

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