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