рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рдлрд▓рди 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 рдХреЛ рдмрд╛рд░-рдмрд╛рд░ рдЖрдзрд╛ рдХрд░рддрд╛ рд╣реИ (рдиреЛрдбреНрд╕ рдкреБрдирдГ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП) рдЬрдм рддрдХ рдХрд┐ рд▓рдЧрд╛рддрд╛рд░ рджреЛ рдЕрдиреБрдорд╛рди рд╡рд╛рдВрдЫрд┐рдд рд╕рд╣рдирд╢реАрд▓рддрд╛ рддрдХ рд╕рд╣рдордд рди рд╣реЛ рдЬрд╛рдПрдБред
рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг
\(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) рдЗрд╕реЗ рд╡рд┐рдлрд▓ рдХрд░ рджреЗрдВрдЧреЗред