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

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

x, рд╕рдВрдХреНрд░рд┐рдпрдХ + - * / ^ рддрдерд╛ рдлрд▓рди exp, log, ln, sin, cos, tan, sqrt, abs, atan рдХрд╛ рдкреНрд░рдпреЛрдЧ рдХрд░реЗрдВред рдЕрдирдВрдд рдкрд░ рдХреНрд╖рдпрд┐рдд рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП (рдЧреИрд░-рдЖрд╡рд░реНрддреА)ред

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

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

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

Integral over (-тИЮ, тИЮ)
approximate value of тИл f(x) dx
рд╡рд┐рдзрд┐ рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (tanh-sinh рдкрд░рд┐рд╡рд╛рд░) рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░
рд░реВрдкрд╛рдВрддрд░рдг x = sinh((╧А/2) sinh t)
рд╡рд╛рдВрдЫрд┐рдд рдЕрдВрдХ
рдЕрднрд┐рд╕рд░рд┐рдд No (may diverge or need more levels)
DE рдпреЛрдЧ рд╡рд╛рдВрдЫрд┐рдд рд╕рд╣рдирд╢реАрд▓рддрд╛ рддрдХ рдирд╣реАрдВ рдкрд╣реБрдБрдЪрд╛ред рд╣реЛ рд╕рдХрддрд╛ рд╣реИ рд╕рдорд╛рдХрд▓реНрдп рдЕрдирдВрдд рдкрд░ рдХреНрд╖рдпрд┐рдд рди рд╣реЛрддрд╛ рд╣реЛ (рдЕрдкрд╕рд╛рд░реА рд╕рдорд╛рдХрд▓), рдЖрд╡рд░реНрддреА/рджреЛрд▓рд╛рдпрдорд╛рди рд╣реЛ, рдпрд╛ рдХрдо рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдХреА рдЖрд╡рд╢реНрдпрдХрддрд╛ рд╣реЛред DE рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдПрдХ рдЧреИрд░-рдЖрд╡рд░реНрддреА, рдЕрдВрддрдмрд┐рдВрджреБ-рд╡рд┐рд╢реНрд▓реЗрд╖рдгрд╛рддреНрдордХ рд╕рдорд╛рдХрд▓реНрдп рдорд╛рдирддрд╛ рд╣реИ рдЬрд┐рд╕рдХрд╛ рдХреНрд╖рдп рдмреАрдЬрдЧрдгрд┐рддреАрдп (algebraic) рд╣реЛред

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

рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рдлрд▓рди f(x) рдХрд╛ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓ рдкреВрд░реА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдкрд░, рдпрд╛рдиреА рдЛрдг рдЕрдирдВрдд рд╕реЗ рдзрди рдЕрдирдВрдд рддрдХ, рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (DE) рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдХреЗ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рдЬреНрдЮрд╛рдд рдХрд░рддрд╛ рд╣реИ тАФ рдЬрд┐рд╕реЗ tanh-sinh рдкрд░рд┐рд╡рд╛рд░ рдпрд╛ Takahasi-Mori рд╡рд┐рдзрд┐ рднреА рдХрд╣рддреЗ рд╣реИрдВред рдЪрд┐рдХрдиреЗ (smooth) рд╕рдорд╛рдХрд▓реНрдпреЛрдВ (integrands) рдХреЗ рд▓рд┐рдП DE рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рд╕рдмрд╕реЗ рдХреБрд╢рд▓ рд╕рд╛рдорд╛рдиреНрдп-рдЙрджреНрджреЗрд╢реНрдп рд╡рд┐рдзрд┐рдпреЛрдВ рдореЗрдВ рд╕реЗ рдПрдХ рд╣реИ рдФрд░ рдЗрд╕рдХрд╛ рдЕрднрд┐рд╕рд░рдг рдЖрд╢реНрдЪрд░реНрдпрдЬрдирдХ рд░реВрдк рд╕реЗ рддреЗрдЬрд╝ рд╣реЛрддрд╛ рд╣реИ: рд╕рд╣реА рдЕрдВрдХреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ рдирдореВрдирд╛ рдмрд┐рдВрджреБрдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ рдХреЗ рд▓рдЧрднрдЧ рд░реИрдЦрд┐рдХ рдЕрдиреБрдкрд╛рдд рдореЗрдВ рдмрдврд╝рддреА рд╣реИред

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

f(x) рдХреЗ рд▓рд┐рдП рдПрдХ рдЧрдгрд┐рддреАрдп рд╡реНрдпрдВрдЬрдХ рд▓рд┐рдЦреЗрдВ рдЬрд┐рд╕рдореЗрдВ x, рд╕рдВрдХреНрд░рд┐рдпрдХ + - * / ^, рдХреЛрд╖реНрдардХ, рддрдерд╛ exp, log/ln, sin, cos, tan, sqrt, abs рдФрд░ atan рдЬреИрд╕реЗ рдорд╛рдирдХ рдлрд▓рдиреЛрдВ рдХрд╛ рдкреНрд░рдпреЛрдЧ рдХрд░ рд╕рдХрддреЗ рд╣реИрдВред рдлрд┐рд░ рддрдп рдХрд░реЗрдВ рдХрд┐ рдЖрдкрдХреЛ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ (6 рд╕реЗ 50) рдЪрд╛рд╣рд┐рдПред рдЕрдзрд┐рдХ рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдХреЗ рд▓рд┐рдП рдЕрдзрд┐рдХ рдмрд╛рд░реАрдХ рд╕реНрдЯреЗрдк рд╕рд╛рдЗрдЬрд╝ рдФрд░ рдЕрдзрд┐рдХ рдЪреМрдбрд╝рд╛ рдкрд░рд┐рдЪреНрдЫреЗрджрди (truncation) рдкрд░рд╛рд╕ рдкреНрд░рдпреБрдХреНрдд рд╣реЛрддрд╛ рд╣реИред рд╕рдорд╛рдХрд▓реНрдп рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдкрд░ рд╡рд┐рд╢реНрд▓реЗрд╖рдгрд╛рддреНрдордХ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП рдФрд░ x рдмрдврд╝рдиреЗ рдкрд░ рдХреНрд╖рдпрд┐рдд (decay) рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП; рдпрд╣ рдЖрд╡рд░реНрддреА (periodic) рдпрд╛ рдмрд┐рдирд╛ рдХреНрд╖рдп рдХреЗ рджреЛрд▓рд╛рдпрдорд╛рди рдирд╣реАрдВ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдПред

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

DE рд╡рд┐рдзрд┐ рдЪрд░-рдкрд░рд┐рд╡рд░реНрддрди \(x = \sinh\!\left(\tfrac{\pi}{2}\sinh t\right)\) рд▓рд╛рдЧреВ рдХрд░рддреА рд╣реИ, рдЬрд┐рд╕рдХрд╛ рдЕрд╡рдХрд▓рдЬ \(\phi'(t) = \tfrac{\pi}{2}\cosh t\,\cosh\!\left(\tfrac{\pi}{2}\sinh t\right)\) рд╣реИред рдкреНрд░рддрд┐рд╕реНрдерд╛рдкрди рдХреЗ рдмрд╛рдж рд╕рдорд╛рдХрд▓, t рдХреЗ рд╕рд╛рдкреЗрдХреНрд╖ \(f(\phi(t))\,\phi'(t)\) рдХрд╛ рд╕рдорд╛рдХрд▓ рдмрди рдЬрд╛рддрд╛ рд╣реИред рдЪреВрдБрдХрд┐ |t| рдмрдврд╝рдиреЗ рдкрд░ \(\phi'(t)\) рдбрдмрд▓ рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ рдХреА рддрд░рд╣ рдХреНрд╖рдпрд┐рдд рд╣реЛрддрд╛ рд╣реИ, рд╕рдорд╛рди рд╕реНрдЯреЗрдк h рд╡рд╛рд▓рд╛ рд╕рд░рд▓ рд╕рдорд▓рдВрдм (trapezoidal) рдирд┐рдпрдо,

$$I \approx h \cdot \sum f(\phi(kh))\,\phi'(kh)$$

рдЕрддреНрдпрдВрдд рд╕рдЯреАрдХ рд╣реЛ рдЬрд╛рддрд╛ рд╣реИред рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдПрдХ рдореЛрдЯреЗ рд╕реНрдЯреЗрдк рд╕реЗ рд╢реБрд░реВ рдХрд░рддрд╛ рд╣реИ, рдлрд┐рд░ h рдХреЛ рдмрд╛рд░-рдмрд╛рд░ рдЖрдзрд╛ рдХрд░рддрд╛ рд╣реИ (рдиреЛрдбреНрд╕ рдкреБрдирдГ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП) рдЬрдм рддрдХ рдХрд┐ рд▓рдЧрд╛рддрд╛рд░ рджреЛ рдЕрдиреБрдорд╛рди рд╡рд╛рдВрдЫрд┐рдд рд╕рд╣рдирд╢реАрд▓рддрд╛ рддрдХ рд╕рд╣рдордд рди рд╣реЛ рдЬрд╛рдПрдБред

рдПрдХ рд╕рдордорд┐рдд рдШрдВрдЯреА-рдЖрдХрд╛рд░ рд╡рдХреНрд░ рдЬреЛ рджреЛрдиреЛрдВ рдУрд░ рдмрд╣реБрдд рддреЗрдЬрд╝реА рд╕реЗ рдХреНрд╖рдп рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рднрд╛рд░ рдХреЛ рджрд░реНрд╢рд╛рддрд╛ рд╣реИ
рднрд╛рд░ ╧Ж'(t) рджреНрд╡рд┐-рдШрд╛рддрд╛рдВрдХреАрдп рд░реВрдк рд╕реЗ рдХреНрд╖рдп рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рджреВрд░ рдХреЗ рдкрдж рд▓рдЧрднрдЧ рдХреБрдЫ рднреА рдпреЛрдЧрджрд╛рди рдирд╣реАрдВ рджреЗрддреЗред
t рдореЗрдВ рд╕рдорд╛рди рджреВрд░реА рд╡рд╛рд▓реЗ рдмрд┐рдВрджреБ x рдореЗрдВ рд╢реВрдиреНрдп рдХреЗ рдкрд╛рд╕ рдЬрдорддреЗ рдФрд░ рджреЛрдиреЛрдВ рдЕрдирдВрдд рдХреА рдУрд░ рдлреИрд▓рддреЗ рдмрд┐рдВрджреБрдУрдВ рдкрд░ рдореИрдк рд╣реЛрддреЗ рд╣реБрдП
DE рд░реВрдкрд╛рдВрддрд░рдг t рдореЗрдВ рдПрдХрд╕рдорд╛рди рдЧреНрд░рд┐рдб рдХреЛ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдкрд░ рдореИрдк рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдмрд┐рдВрджреБ рдХреЗрдВрджреНрд░ рдХреЗ рдкрд╛рд╕ рд╕рдВрдХреЗрдВрджреНрд░рд┐рдд рд╣реЛрддреЗ рд╣реИрдВред

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

\(f(x) = \tfrac{1}{1+x^2}\) рдХреЗ рд▓рд┐рдП рд╕рдЯреАрдХ рд╕рдорд╛рдХрд▓

$$\int_{-\infty}^{\infty} \frac{1}{1+x^2}\,dx = \left[\arctan x\right]_{-\infty}^{\infty} = \pi \approx 3.14159265358979$$

рд╣реИред рд╕реНрдЯреЗрдк 0.5 рдФрд░ рдиреЛрдбреНрд╕ \(k = -8..8\) рдХреЗ рд╕рд╛рде рдПрдХ рдореЛрдЯрд╛ DE рдпреЛрдЧ рдкрд╣рд▓реЗ рд╣реА рд▓рдЧрднрдЧ \(3.15\) рджреЗрддрд╛ рд╣реИ; рд╕реНрдЯреЗрдк рдХреЛ рдмрд╛рд░реАрдХ рдХрд░рдиреЗ рдкрд░ рдпрд╣ pi рдХреЗ рдкреВрд░реНрдг-рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдорд╛рди рддрдХ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛ рдЬрд╛рддрд╛ рд╣реИред рдЗрд╕реА рддрд░рд╣ \(\exp(-x^2)\) рдХрд╛ рдкрд░рд┐рдгрд╛рдо \(\sqrt{\pi} \approx 1.77245385090552\) рдорд┐рд▓рддрд╛ рд╣реИред

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

рдпрд╣ "рдЕрднрд┐рд╕рд░рд┐рдд рдирд╣реАрдВ рд╣реБрдЖ" рдХреНрдпреЛрдВ рдХрд╣рддрд╛ рд╣реИ? рдЕрдзрд┐рдХрддрд░ рддрдм рдЬрдм рд╕рдорд╛рдХрд▓ рдЕрдкрд╕рд╛рд░реА (diverge) рд╣реЛрддрд╛ рд╣реИ (рд╕рдорд╛рдХрд▓реНрдп рдХреНрд╖рдпрд┐рдд рдирд╣реАрдВ рд╣реЛрддрд╛, рдЬреИрд╕реЗ \(f = 1\)) рдпрд╛ рдлрд▓рди рдЖрд╡рд░реНрддреА/рджреЛрд▓рд╛рдпрдорд╛рди рд╣реЛрддрд╛ рд╣реИред DE рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдПрдХ рдЧреИрд░-рдЖрд╡рд░реНрддреА, рдЕрдВрддрдмрд┐рдВрджреБ-рд╡рд┐рд╢реНрд▓реЗрд╖рдгрд╛рддреНрдордХ рд╕рдорд╛рдХрд▓реНрдп рдорд╛рдирдХрд░ рдЪрд▓рддрд╛ рд╣реИред

рдореБрдЭреЗ рдХрд┐рддрдиреА рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдЪреБрдирдиреА рдЪрд╛рд╣рд┐рдП? 15 рдЕрдВрдХ рдбрдмрд▓-рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдХреЗ рдмрд░рд╛рдмрд░ рд╣реИрдВ рдФрд░ рдПрдХ рдЕрдЪреНрдЫрд╛ рдбрд┐рдлрд╝реЙрд▓реНрдЯ рд╣реИрдВред рдбрдмрд▓-рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдЬрд┐рддрдирд╛ рд╕рдорд░реНрдерди рдХрд░ рд╕рдХреЗ рдЙрд╕рд╕реЗ рдХрд╣реАрдВ рдЕрдзрд┐рдХ рдЕрдВрдХ рдорд╛рдБрдЧрдиреЗ рд╕реЗ рдЙрддреНрддрд░ рдмреЗрд╣рддрд░ рдирд╣реАрдВ рд╣реЛрдЧрд╛ред

рдХреНрдпрд╛ рдореИрдВ рд╡рд┐рдЪрд┐рддреНрд░ (singular) рдлрд▓рдиреЛрдВ рдХрд╛ рд╕рдорд╛рдХрд▓рди рдХрд░ рд╕рдХрддрд╛ рд╣реВрдБ? рдпрд╣ рд╡рд┐рдзрд┐ рдЕрдВрддрдмрд┐рдВрджреБ рд╡реНрдпрд╡рд╣рд╛рд░ рдХреЛ рд╕рд╣ рд▓реЗрддреА рд╣реИ рдХреНрдпреЛрдВрдХрд┐ рдЪрд░рдо рдиреЛрдбреНрд╕ рдХрд╛ рднрд╛рд░ рд▓рдЧрднрдЧ рд╢реВрдиреНрдп рд╣реЛрддрд╛ рд╣реИ, рдкрд░рдиреНрддреБ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдЕрдХреНрд╖ рдкрд░ рдЖрдВрддрд░рд┐рдХ рдзреНрд░реБрд╡ (poles) рдЗрд╕реЗ рд╡рд┐рдлрд▓ рдХрд░ рджреЗрдВрдЧреЗред

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