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

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

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

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

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

[a, рдЕрдирдВрдд) рдкрд░ рд╕рдорд╛рдХрд▓рди
рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (DE) рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рджреНрд╡рд╛рд░рд╛ рдкрд░рд┐рдХрд▓рд┐рдд
рд╡рд┐рдзрд┐ рдЖрдзреА-рд░реЗрдЦрд╛ DE рд░реВрдкрд╛рдВрддрд░рдг x = a + exp((pi/2) sinh t)
рдЕрдВрддрд┐рдо рдЪрд░рдг h
рдЕрднрд┐рд╕рд░рд┐рдд No (check convergence / decay)
рд╕рдорд▓рдВрдм рдпреЛрдЧ рдЕрдиреБрд░реЛрдзрд┐рдд рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рддрдХ рд╕реНрдерд┐рд░ рдирд╣реАрдВ рд╣реБрдЖред рд╣реЛ рд╕рдХрддрд╛ рд╣реИ рдХрд┐ рд╕рдорд╛рдХрд▓реНрдп рдХреНрд╖рдп рди рд╣реЛ рд░рд╣рд╛ рд╣реЛ (рдЕрдкрд╕рд╛рд░реА рд╕рдорд╛рдХрд▓рди), рджреЛрд▓рдирд╢реАрд▓/рдЖрд╡рд░реНрддреА рд╣реЛ, рдпрд╛ рдЗрд╕реЗ рдХрд┐рд╕реА рдЕрд▓рдЧ DE рдореИрдк рдХреА рдЖрд╡рд╢реНрдпрдХрддрд╛ рд╣реЛред рдРрд╕реЗ рдЗрдирдкреБрдЯ рдХреЗ рдкрд░рд┐рдгрд╛рдо рдЕрд╡рд┐рд╢реНрд╡рд╕рдиреАрдп рд╣реЛрддреЗ рд╣реИрдВред

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

рдпрд╣ рдЯреВрд▓ рдЕрд░реНрдз-рдЕрдирдВрдд рдЕрдВрддрд░рд╛рд▓ рдкрд░ рдПрдХ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓рди рдХрд╛ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рдорд╛рди рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИ, \(I = \int_{a}^{\infty} f(x)\,dx\), рдЬрд┐рд╕рдХреЗ рд▓рд┐рдП рдпрд╣ рдЖрдзреА-рд░реЗрдЦрд╛ (half-line) рдХреЗ рд▓рд┐рдП рд╡рд┐рд╢реЗрд╖ рд░реВрдк рд╕реЗ рдмрдиреА рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (DE) рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред DE рд╡рд┐рдзрд┐ рдПрдХ рд╕рд╛рд░реНрд╡рднреМрдорд┐рдХ рдФрд░ рдЙрдЪреНрдЪ-рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рд╡рд╛рд▓реА рдпреЛрдЬрдирд╛ рд╣реИ, рдЬреЛ рд╡рд┐рд╢реЗрд╖ рд░реВрдк рд╕реЗ рдЙрди рд╕рдорд╛рдХрд▓реНрдпреЛрдВ рдХреЗ рд▓рд┐рдП рдЙрдкрдпреБрдХреНрдд рд╣реИ рдЬреЛ рдЕрдирдВрдд рдХреА рдУрд░ рдмреАрдЬрдЧрдгрд┐рддреАрдп рд░реВрдк рд╕реЗ (рдкрд╛рд╡рд░-рд▓реЙ, рдЬреИрд╕реЗ \(1/x^{p}\)) рдХреНрд╖рдп рд╣реЛрддреЗ рд╣реИрдВ рдФрд░ рдЬрд┐рдирдХреЗ рдирд┐рдЪрд▓реЗ рд╕рд┐рд░реЗ \(x = a\) рдкрд░ рд╣рд▓реНрдХреА рд╡рд┐рдЪрд┐рддреНрд░рддрд╛ (singularity) рд╣реЛ рд╕рдХрддреА рд╣реИред рдпрд╣ рдЖрд╡рд░реНрддреА (periodic) рдпрд╛ рджреЛрд▓рдирд╢реАрд▓ (oscillatory) рд╕рдорд╛рдХрд▓реНрдпреЛрдВ рдХреЗ рд▓рд┐рдП рдирд╣реАрдВ рдмрдиреА рд╣реИред

a рд╕реЗ рдЕрдирдВрдд рддрдХ рдЕрд░реНрдз-рдЕрдирдВрдд рдЕрдВрддрд░рд╛рд▓ рдкрд░ рд╡рдХреНрд░ рдХреЗ рдиреАрдЪреЗ рдХрд╛ рдХреНрд╖реЗрддреНрд░рдлрд▓
рдХреИрд▓рдХреБрд▓реЗрдЯрд░ f(x) рдХреЗ рдиреАрдЪреЗ рдХрд╛ рдХреНрд╖реЗрддреНрд░рдлрд▓ x = a рд╕реЗ рдЕрдирдВрдд рддрдХ рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИред

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

рд╕рдорд╛рдХрд▓реНрдп f(x) рдХреЛ рдЪрд░ x рдХреЗ рдкрдж рдореЗрдВ рдПрдХ рдЧрдгрд┐рддреАрдп рд╡реНрдпрдВрдЬрдХ рдХреЗ рд░реВрдк рдореЗрдВ рджрд░реНрдЬ рдХрд░реЗрдВ (рдпрд╣ + - * / ^, рдХреЛрд╖реНрдардХ, рдФрд░ sqrt, exp, log, log10, sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, abs рдХреЗ рд╕рд╛рде-рд╕рд╛рде рдЕрдЪрд░ pi рдФрд░ e рдХрд╛ рд╕рдорд░реНрдерди рдХрд░рддрд╛ рд╣реИ)ред рдкрд░рд┐рдорд┐рдд рдирд┐рдЪрд▓реА рд╕реАрдорд╛ a рджрд░реНрдЬ рдХрд░реЗрдВ, рдЪреБрдиреЗрдВ рдХрд┐ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ (significant digits) рджрд┐рдЦрд╛рдиреЗ рд╣реИрдВ, рдФрд░ рд╕рдмрдорд┐рдЯ рдХрд░реЗрдВред рдкрд░рд┐рдгрд╛рдо рд╕рдорд╛рдХрд▓рди рдХрд╛ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рдорд╛рди рд╣реЛрдЧрд╛ред

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

рдЪрд░-рдкрд░рд┐рд╡рд░реНрддрди \(\phi(t) = \exp\!\left(\tfrac{\pi}{2}\cdot\sinh t\right)\) рдХреЗ рджреНрд╡рд╛рд░рд╛ рдЖрдзреА-рд░реЗрдЦрд╛ рдХреЛ рдкреВрд░реА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдкрд░ рдореИрдк рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рддрд╛рдХрд┐ \(x = a + \phi(t)\) рдорд╛рди \(a\) (рдЬрдм \(t \to -\infty\)) рд╕реЗ \(+\infty\) рддрдХ рдлреИрд▓реЗред рдЗрд╕рдХрд╛ рдЕрд╡рдХрд▓рдЬ \(\phi^{\prime}(t) = \tfrac{\pi}{2}\cdot\cosh(t)\cdot\exp\!\left(\tfrac{\pi}{2}\cdot\sinh t\right)\) рд╣реЛрддрд╛ рд╣реИред рд░реВрдкрд╛рдВрддрд░рд┐рдд рд╕рдорд╛рдХрд▓реНрдп \(g(t) = f(a + \phi(t))\cdot\phi^{\prime}(t)\) рджреЛрдиреЛрдВ рд╕рд┐рд░реЛрдВ рдкрд░ рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ рд░реВрдк рд╕реЗ рдХреНрд╖рдп рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рдПрдХрд╕рдорд╛рди рдЧреНрд░рд┐рдб рдкрд░ рд╕рд░рд▓ рд╕рдорд▓рдВрдм рдирд┐рдпрдо (trapezoidal rule) рд▓рдЧрднрдЧ рд╕рд░реНрд╡реЛрддреНрддрдо рд╣реЛрддрд╛ рд╣реИ:

$$I \approx h\cdot\sum_{k=-N}^{N} g(kh)$$

рд╡реЗ рдиреЛрдб рдЬреЛ рдУрд╡рд░рдлреНрд▓реЛ рдХрд░рддреЗ рд╣реИрдВ (рдмрдбрд╝реЗ \(t\) рдкрд░) рдЫреЛрдбрд╝ рджрд┐рдП рдЬрд╛рддреЗ рд╣реИрдВ рдХреНрдпреЛрдВрдХрд┐ рд╡рд╣рд╛рдБ рд╕рдорд╛рдХрд▓реНрдп рд╡реНрдпрд╛рд╡рд╣рд╛рд░рд┐рдХ рд░реВрдк рд╕реЗ рд╢реВрдиреНрдп рд╣реЛрддрд╛ рд╣реИред

рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдХреЛ рдЕрд░реНрдз-рдЕрдирдВрдд рдЕрдВрддрд░рд╛рд▓ рдкрд░ рдореИрдк рдХрд░рдиреЗ рд╡рд╛рд▓рд╛ рджреНрд╡рд┐-рдШрд╛рддрд╛рдВрдХреАрдп рдЪрд░ рдкрд░рд┐рд╡рд░реНрддрди
DE рд░реВрдкрд╛рдВрддрд░рдг рдкреВрд░реА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдкрд░ t рдХреЛ [a, рдЕрдирдВрдд) рдореЗрдВ x рдкрд░ рдореИрдк рдХрд░рддрд╛ рд╣реИ, рдиреЛрдбреНрд╕ рдХреЛ рдХреБрд╢рд▓рддрд╛ рд╕реЗ рдЬрдорд╛рддреЗ рд╣реБрдПред

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

\(f(x) = 1/(1 + x^{3/2})\) рдФрд░ \(a = 0\) рдХреЗ рд▓рд┐рдП рд╕рдЯреАрдХ рдорд╛рди \(\tfrac{4\pi}{3\sqrt{3}} \approx 2.4183991523\) рд╣реИред \(h = 1/16\) рдХреЗ рд╕рд╛рде DE рдпреЛрдЧ рдЗрд╕реЗ рдХрдИ рдЕрдВрдХреЛрдВ рддрдХ рдкреБрдирдГ рдкреНрд░рд╕реНрддреБрдд рдХрд░ рджреЗрддрд╛ рд╣реИред \(t = 0\) рдкрд░: \(\phi = 1\), \(x = 1\), \(f = 0.5\), \(\phi^{\prime} = \pi/2 \approx 1.5708\), рдЗрд╕рд▓рд┐рдП \(g \approx 0.7854\); рд╕рднреА рднрд╛рд░рд┐рдд рдиреЛрдб рдХреЛ рдЬреЛрдбрд╝рдиреЗ рдкрд░ рдпрд╣ \(2.41840\) рдкрд░ рдЕрднрд┐рд╕рд░рд┐рдд (converge) рд╣реЛ рдЬрд╛рддрд╛ рд╣реИред

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

рдХреНрдпрд╛ рдореИрдВ рдЗрд╕реЗ рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ рд░реВрдк рд╕реЗ рдХреНрд╖рдпрд┐рдд рд╕рдорд╛рдХрд▓реНрдпреЛрдВ рдХреЗ рд▓рд┐рдП рдЙрдкрдпреЛрдЧ рдХрд░ рд╕рдХрддрд╛ рд╣реВрдБ? рдпрд╣ рддрдм рднреА рдХрд╛рдо рдХрд░рддрд╛ рд╣реИ, рд▓реЗрдХрд┐рди рд╢реБрджреНрдз рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ рдХреНрд╖рдп рдХреЗ рд▓рд┐рдП рдПрдХ рдЕрд▓рдЧ DE рдореИрдк рддреЗрдЬрд╝реА рд╕реЗ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛрддрд╛ рд╣реИ; рдпрд╣ рд╕рдВрд╕реНрдХрд░рдг рдмреАрдЬрдЧрдгрд┐рддреАрдп рдХреНрд╖рдп рдХреЗ рд▓рд┐рдП рд▓рдХреНрд╖рд┐рдд рд╣реИред

рджреЛрд▓рдирд╢реАрд▓ рд╕рдорд╛рдХрд▓реНрдп рдкрд░ рдпрд╣ рд╡рд┐рдлрд▓ рдХреНрдпреЛрдВ рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ? рдмрд┐рдирд╛ рдХреНрд╖рдп рд╡рд╛рд▓реЗ рджреЛрд▓рди рдХрд╛ рд╕рдорд▓рдВрдм рдпреЛрдЧ рдЕрднрд┐рд╕рд░рд┐рдд рдирд╣реАрдВ рд╣реЛрддрд╛, рдЗрд╕рд▓рд┐рдП рдЖрдзреА-рд░реЗрдЦрд╛ рдХреЗ рд▓рд┐рдП DE рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рд╡рд╣рд╛рдБ рдЙрдкрдпреБрдХреНрдд рдирд╣реАрдВ рд╣реИред

рдЕрдВрдХ (digits) рд╕реЗрдЯрд┐рдВрдЧ рд╕реЗ рдХреНрдпрд╛ рдмрджрд▓рддрд╛ рд╣реИ? рдХреЗрд╡рд▓ рджрд┐рдЦрд╛рдИ рдЬрд╛рдиреЗ рд╡рд╛рд▓реА рд░рд╛рдЙрдВрдбрд┐рдВрдЧ; рдЖрдВрддрд░рд┐рдХ рдЧрдгрдирд╛ рд╣рдореЗрд╢рд╛ рдкреВрд░реНрдг рдбрдмрд▓ рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ (double precision) рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреА рд╣реИред

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