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

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

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

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

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

рдЯреЗрдмрд▓ рдХрд╛ рдкрд╣рд▓рд╛ рдорд╛рди I_v(x), x = рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди рдкрд░
1
51 rows generated
x I_v(x)
0 1
0.1 1.0025015629
0.2 1.0100250278
0.3 1.0226268794
0.4 1.0404017822
0.5 1.0634833707
0.6 1.0920453643
0.7 1.1263030183
0.8 1.1665149229
0.9 1.2129851657
1 1.2660658778
1.1 1.3261601837
1.2 1.3937255841
1.3 1.4692777979
1.4 1.5533950997
1.5 1.6467231898
1.6 1.7499806397
1.7 1.8639649621
1.8 1.9895593566
1.9 2.1277401941
2 2.2795853023
2.1 2.4462831294
2.2 2.6291428636
2.3 2.8296056006
2.4 3.049256658
2.5 3.2898391441
2.6 3.5532689042
2.7 3.8416509766
2.8 4.1572977035
2.9 4.5027486613
3 4.8807925859
3.1 5.2944914897
3.2 5.7472071872
3.3 6.2426304652
3.4 6.7848131604
3.5 7.3782034322
3.6 8.0276845471
3.7 8.7386175242
3.8 9.5168880261
3.9 10.3689579167
4 11.3019219521
4.1 12.323570116
4.2 13.4424561633
4.3 14.6679729918
4.4 16.0104355249
4.5 17.4811718556
4.6 19.0926234795
4.7 20.8584555266
4.8 22.7936779931
4.9 24.9147790758
5 27.2398718236

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

рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рдирд┐рд╢реНрдЪрд┐рдд рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдХреЛрдЯрд┐ v рдХреЗ рд▓рд┐рдП рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░ рдХреЗ рд╕реБрдзрд╛рд░рд┐рдд рдмреЗрд╕реЗрд▓ рдлрд▓рди \(I_{v}(x)\) рдХреА рдЯреЗрдмрд▓ рдмрдирд╛рддрд╛ рд╣реИ, рдЬреЛ x рдорд╛рдиреЛрдВ рдХреЗ рдПрдХ рдХреНрд░рдо рдкрд░ рдкрд░рд┐рдХрд▓рд┐рдд рд╣реЛрддрд╛ рд╣реИред рдЖрдк рдХреЛрдЯрд┐, x рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди, рд╡реГрджреНрдзрд┐ (step) рдФрд░ рдХрд┐рддрдиреА рдкрдВрдХреНрддрд┐рдпрд╛рдБ рдЪрд╛рд╣рд┐рдП тАФ рдпреЗ рджрд░реНрдЬ рдХрд░рддреЗ рд╣реИрдВ; рдХреИрд▓рдХреБрд▓реЗрдЯрд░ \(x_{i} = \text{start} + i \cdot \text{step}\) рдХреА рд╕реВрдЪреА рдмрдирд╛рддрд╛ рд╣реИ рдФрд░ рд╣рд░ рдмрд┐рдВрджреБ рдкрд░ \(I_{v}(x_{i})\) рдХреА рдЧрдгрдирд╛ рдХрд░рдХреЗ рдЯреЗрдмрд▓ рдФрд░ рдЧреНрд░рд╛рдл рджреЛрдиреЛрдВ рджреЗрддрд╛ рд╣реИред рдпрд╣ рд╢реБрджреНрдз рдЧрдгрд┐рдд рдХрд╛ рд╡рд┐рд╢реЗрд╖-рдлрд▓рди рдЯреВрд▓ рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рд╕рдорд╛рди рд░реВрдк рд╕реЗ рд▓рд╛рдЧреВ рд╣реЛрддрд╛ рд╣реИ (рдЗрд╕рдореЗрдВ рдХреЛрдИ рдХреНрд╖реЗрддреНрд░реАрдп рдирд┐рдпрдо рдпрд╛ рдЗрдХрд╛рдЗрдпрд╛рдБ рдирд╣реАрдВ рд╣реИрдВ)ред

рдХрдИ рдХреЛрдЯрд┐рдпреЛрдВ рдХреЗ рд▓рд┐рдП рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░ рдХреЗ рд╕рдВрд╢реЛрдзрд┐рдд рдмреЗрд╕реЗрд▓ рдлрд▓рди рдХреА рдПрдХрджрд┐рд╖реНрдЯ рд░реВрдк рд╕реЗ рдмрдврд╝рддреА рдЪрд░рдШрд╛рддрд╛рдВрдХреА-рдЬреИрд╕реА рд╡рдХреНрд░реЛрдВ рдХрд╛ рд╕рдореВрд╣
рд╕рдВрд╢реЛрдзрд┐рдд рдмреЗрд╕реЗрд▓ рдлрд▓рди \(I_v(x)\) рдХреЛрдЯрд┐ v = 0, 1, 2, 3 рдХреЗ рд▓рд┐рдП x рдХреЗ рд╕рд╛рде рддреЗрдЬрд╝реА рд╕реЗ рдмрдврд╝рддреЗ рд╣реИрдВред

рд╕реВрддреНрд░

рд╕реБрдзрд╛рд░рд┐рдд рдмреЗрд╕реЗрд▓ рдлрд▓рди \(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\)ред

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