рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рдирд┐рд╢реНрдЪрд┐рдд рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдХреЛрдЯрд┐ v рдХреЗ рд▓рд┐рдП рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░ рдХреЗ рд╕реБрдзрд╛рд░рд┐рдд рдмреЗрд╕реЗрд▓ рдлрд▓рди \(I_{v}(x)\) рдХреА рдЯреЗрдмрд▓ рдмрдирд╛рддрд╛ рд╣реИ, рдЬреЛ x рдорд╛рдиреЛрдВ рдХреЗ рдПрдХ рдХреНрд░рдо рдкрд░ рдкрд░рд┐рдХрд▓рд┐рдд рд╣реЛрддрд╛ рд╣реИред рдЖрдк рдХреЛрдЯрд┐, x рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди, рд╡реГрджреНрдзрд┐ (step) рдФрд░ рдХрд┐рддрдиреА рдкрдВрдХреНрддрд┐рдпрд╛рдБ рдЪрд╛рд╣рд┐рдП тАФ рдпреЗ рджрд░реНрдЬ рдХрд░рддреЗ рд╣реИрдВ; рдХреИрд▓рдХреБрд▓реЗрдЯрд░ \(x_{i} = \text{start} + i \cdot \text{step}\) рдХреА рд╕реВрдЪреА рдмрдирд╛рддрд╛ рд╣реИ рдФрд░ рд╣рд░ рдмрд┐рдВрджреБ рдкрд░ \(I_{v}(x_{i})\) рдХреА рдЧрдгрдирд╛ рдХрд░рдХреЗ рдЯреЗрдмрд▓ рдФрд░ рдЧреНрд░рд╛рдл рджреЛрдиреЛрдВ рджреЗрддрд╛ рд╣реИред рдпрд╣ рд╢реБрджреНрдз рдЧрдгрд┐рдд рдХрд╛ рд╡рд┐рд╢реЗрд╖-рдлрд▓рди рдЯреВрд▓ рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рд╕рдорд╛рди рд░реВрдк рд╕реЗ рд▓рд╛рдЧреВ рд╣реЛрддрд╛ рд╣реИ (рдЗрд╕рдореЗрдВ рдХреЛрдИ рдХреНрд╖реЗрддреНрд░реАрдп рдирд┐рдпрдо рдпрд╛ рдЗрдХрд╛рдЗрдпрд╛рдБ рдирд╣реАрдВ рд╣реИрдВ)ред
рд╕реВрддреНрд░
рд╕реБрдзрд╛рд░рд┐рдд рдмреЗрд╕реЗрд▓ рдлрд▓рди \(I_{v}(x)\) рд╕реБрдзрд╛рд░рд┐рдд рдмреЗрд╕реЗрд▓ рд╕рдореАрдХрд░рдг \(x^2 y'' + x y' - (x^2 + v^2) y = 0\) рдХрд╛ рд╣рд▓ рд╣реИред рдЗрд╕реЗ рдпрд╣рд╛рдБ рдЗрд╕рдХреА рдШрд╛рдд-рд╢реНрд░реЗрдгреА (power series) рд╕реЗ рдкрд░рд┐рдХрд▓рд┐рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ:
$$I_{v}(x) = \sum_{k=0}^{\infty} \frac{1}{k!\;\Gamma(v+k+1)}\left(\frac{x}{2}\right)^{v+2k}$$рдХреНрд░рдордЧреБрдгрд┐рдд (factorial) рдФрд░ рдЧрд╛рдорд╛ рдлрд▓рди рдХреА рд╡рдЬрд╣ рд╕реЗ v рдХреЛрдИ рднреА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реЛ рд╕рдХрддреА рд╣реИред рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕реНрдерд┐рд░рддрд╛ рдХреЗ рд▓рд┐рдП рд╣рд░ рдкрдж рдХреЛ log рд╕реНрдкреЗрд╕ рдореЗрдВ, \(\ln \Gamma\) рдХреЗ Lanczos рд╕рдиреНрдирд┐рдХрдЯрди рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдХреЗ рдкрд░рд┐рдХрд▓рд┐рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдФрд░ рддрдм рддрдХ рдЬреЛрдбрд╝рд╛ рдЬрд╛рддрд╛ рд╣реИ рдЬрдм рддрдХ рдкрдж рдирдЧрдгреНрдп рди рд╣реЛ рдЬрд╛рдПрдБред
рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рдХреЛрдЯрд┐ v (рдЬреИрд╕реЗ 0, 1, рдпрд╛ 2.5), x рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди, рд╣рд░ рдкрдВрдХреНрддрд┐ рдореЗрдВ x рдореЗрдВ рдЬреБрдбрд╝рдиреЗ рд╡рд╛рд▓реА рд╡реГрджреНрдзрд┐, рдФрд░ рдкрдВрдХреНрддрд┐рдпреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ рджрд░реНрдЬ рдХрд░реЗрдВред рдЧрдгрдирд╛ рдХрд░реЗрдВ рдкрд░ рдХреНрд▓рд┐рдХ рдХрд░рддреЗ рд╣реА рдЖрдкрдХреЛ x рдФрд░ \(I_{v}(x)\) рдХреА рджреЛ-рд╕реНрддрдВрдн рд╡рд╛рд▓реА рдЯреЗрдмрд▓ рдФрд░ рдЙрд╕реА рд░реЗрдВрдЬ рдкрд░ рдПрдХ рдЧреНрд░рд╛рдл рдорд┐рд▓ рдЬрд╛рдПрдЧрд╛ред
рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг
v = 0, start = 0, step = 0.5, count = 5 рд▓реЗрдиреЗ рдкрд░ рдЖрдкрдХреЛ x = 0, 0.5, 1, 1.5, 2 рдорд┐рд▓рддреЗ рд╣реИрдВ рдФрд░:
$$I_{0}(0) = 1,\quad I_{0}(0.5) \approx 1.0634834,\quad I_{0}(1) \approx 1.2660658,\quad I_{0}(1.5) \approx 1.6467232,\quad I_{0}(2) \approx 2.2795853$$рдпреЗ рдорд╛рди рдорд╛рдирдХ рд╕рдВрджрд░реНрдн рдЯреЗрдмрд▓реЛрдВ рд╕реЗ рдореЗрд▓ рдЦрд╛рддреЗ рд╣реИрдВред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди
рдХреНрдпрд╛ рдХреЛрдЯрд┐ рдЛрдгрд╛рддреНрдордХ рдпрд╛ рдЕрдкреВрд░реНрдгрд╛рдВрдХ рд╣реЛ рд╕рдХрддреА рд╣реИ? рд╣рд╛рдБред рдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХ рдХреЛрдЯрд┐ рдХреЗ рд▓рд┐рдП рд╕рд░реНрд╡рд╕рдорд┐рдХрд╛ \(I_{-n}(x) = I_{n}(x)\) рдХрд╛ рдЙрдкрдпреЛрдЧ рд╣реЛрддрд╛ рд╣реИред рдЕрдкреВрд░реНрдгрд╛рдВрдХ v рдХреЗ рд▓рд┐рдП \(x \geq 0\) рд╕рдорд░реНрдерд┐рдд рд╣реИ; рдЬрдмрдХрд┐ \(x < 0\) рдФрд░ рдЕрдкреВрд░реНрдгрд╛рдВрдХ v рд╣реЛрдиреЗ рдкрд░ рдорд╛рди рд╕рдореНрдорд┐рд╢реНрд░ (complex) рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП NaN рд▓реМрдЯрд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
\(I_{v}(x)\) рдЗрддрдиреА рддреЗрдЬрд╝реА рд╕реЗ рдХреНрдпреЛрдВ рдмрдврд╝рддрд╛ рд╣реИ? рд╕рд╛рдзрд╛рд░рдг рдмреЗрд╕реЗрд▓ рдлрд▓рди \(J_{v}\) рдХреА рддрд░рд╣ рджреЛрд▓рди рдХрд░рдиреЗ рдХреЗ рдмрдЬрд╛рдп, рд╕реБрдзрд╛рд░рд┐рдд рдлрд▓рди рдмрдбрд╝реЗ x рдХреЗ рд▓рд┐рдП рд▓рдЧрднрдЧ \(e^{x}/\sqrt{2\pi x}\) рдХреА рджрд░ рд╕реЗ рдмрдврд╝рддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рдмрд╣реБрдд рдмрдбрд╝реЗ x рдкрд░ рдорд╛рди рдЕрдирдВрдд (infinity) рддрдХ overflow рд╣реЛ рд╕рдХрддрд╛ рд╣реИред
\(I_{v}(0)\) рдХрд┐рддрдирд╛ рд╣реЛрддрд╛ рд╣реИ? \(I_{0}(0) = 1\), рдФрд░ v > 0 рдХреЗ рд▓рд┐рдП \(I_{v}(0) = 0\)ред