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

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

This rebuilt version computes integer order v тЙе 0 for real x > 0. (The original tool also supports complex x and arbitrary real order.)

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

рд╕реВрддреНрд░ (рдлреЙрд░реНрдореВрд▓рд╛): рдЧреЛрд▓реАрдп рдмреЗрд╕реЗрд▓ рдлрд▓рди рдХреИрд▓рдХреБрд▓реЗрдЯрд░

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

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

рдЧреЛрд▓реАрдп рдмреЗрд╕реЗрд▓ jv(x)
0.4546487134
рдкрд╣рд▓рд╛ рдкреНрд░рдХрд╛рд░
рдлрд▓рди рдорд╛рди
jv(x) тАФ first kind 0.4546487134
yv(x) тАФ second kind 0.2080734183
j'v(x) тАФ derivative -0.435397775
y'v(x) тАФ derivative 0.3506120043

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

рдпрд╣ рдЯреВрд▓ рдкрд╣рд▓реЗ рдкреНрд░рдХрд╛рд░ рдХреЗ рдЧреЛрд▓реАрдп рдмреЗрд╕реЗрд▓ рдлрд▓рди jv(x), рджреВрд╕рд░реЗ рдкреНрд░рдХрд╛рд░ рдХреЗ рдЧреЛрд▓реАрдп рдмреЗрд╕реЗрд▓ рдлрд▓рди yv(x), рддрдерд╛ рдЗрдирдХреЗ рдкреНрд░рдердо рдЕрд╡рдХрд▓рдЬреЛрдВ j'v(x) рдФрд░ y'v(x) рдХреА рдЧрдгрдирд╛ рдХрд░рддрд╛ рд╣реИред рдпреЗ рдлрд▓рди рдЧреЛрд▓реАрдп рдирд┐рд░реНрджреЗрд╢рд╛рдВрдХреЛрдВ рдореЗрдВ рддрд░рдВрдЧ рд╕рдореАрдХрд░рдг рдФрд░ рд╣реЗрд▓реНрдорд╣реЛрд▓реНрдЯреНрдЬрд╝ рд╕рдореАрдХрд░рдг рдХреЗ рддреНрд░рд┐рдЬреНрдп (рд░реЗрдбрд┐рдпрд▓) рд╣рд▓ рд╣реИрдВ, рдФрд░ рднреМрддрд┐рдХреА рдореЗрдВ рд╣рд░ рдЬрдЧрд╣ рджрд┐рдЦрд╛рдИ рджреЗрддреЗ рд╣реИрдВ тАФ рдкреНрд░рдХреАрд░реНрдгрди рд╕рд┐рджреНрдзрд╛рдВрдд, рд╡рд┐рджреНрдпреБрддрдЪреБрдВрдмрдХреАрдп рдПрд╡рдВ рдзреНрд╡рдирд┐рдХ рд╡рд┐рдХрд┐рд░рдг, рддрдерд╛ рдХреНрд╡рд╛рдВрдЯрдо рдпрд╛рдВрддреНрд░рд┐рдХреА (рдореБрдХреНрдд-рдХрдг рдЖрдВрд╢рд┐рдХ рддрд░рдВрдЧреЗрдВ)ред рдЗрд╕ рдирдП рд╕рдВрд╕реНрдХрд░рдг рдореЗрдВ рдкреВрд░реНрдгрд╛рдВрдХ рдХреЛрдЯрд┐ \(v \ge 0\) рдФрд░ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рддрд░реНрдХ \(x > 0\) рдХреЛ рд╕рдВрднрд╛рд▓рд╛ рдЬрд╛рддрд╛ рд╣реИред

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

рдХреЛрдЯрд┐ v (рдХреЛрдИ рдЛрдгреЗрддрд░ рдкреВрд░реНрдгрд╛рдВрдХ рдЬреИрд╕реЗ 0, 1, 2) рдФрд░ рддрд░реНрдХ x (рдХреЛрдИ рдзрдирд╛рддреНрдордХ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛) рджрд░реНрдЬ рдХрд░реЗрдВред рдЪрд╛рд░реЛрдВ рдорд╛рди рдкрд╛рдиреЗ рдХреЗ рд▓рд┐рдП "рдЧрдгрдирд╛ рдХрд░реЗрдВ" рджрдмрд╛рдПрдБред рдзреНрдпрд╛рди рджреЗрдВ рдХрд┐ рдЬрдм x рд╢реВрдиреНрдп рдХреА рдУрд░ рдмрдврд╝рддрд╛ рд╣реИ рддреЛ yv(x) рдФрд░ y'v(x) рдЕрдкрд╕рд░рд┐рдд (diverge) рд╣реЛрддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП рдЯреВрд▓ \(x = 0\) рдкрд░ рдЗрдиреНрд╣реЗрдВ рдЕрдирдВрдд (infinite) рджрд┐рдЦрд╛рддрд╛ рд╣реИ; рдЬрдмрдХрд┐ рд╕реАрдорд╛ рдХреЗ рд░реВрдк рдореЗрдВ j0(0) рдХрд╛ рдорд╛рди 1 рд╣реЛрддрд╛ рд╣реИред

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

рдпреЗ рдлрд▓рди рд╕рдореАрдХрд░рдг $$x^2 w'' + 2x\,w' + (x^2 - v(v+1))w = 0$$ рдХреЛ рд╕рдВрддреБрд╖реНрдЯ рдХрд░рддреЗ рд╣реИрдВред рдмрдВрдж рд░реВрдкреЛрдВ \(j_0 = \sin(x)/x\) рдФрд░ \(y_0 = -\cos(x)/x\) рд╕реЗ рд╢реБрд░реВ рдХрд░рдХреЗ, рдЙрдЪреНрдЪрддрд░ рдХреЛрдЯрд┐рдпрд╛рдБ рддреАрди-рдкрдж рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ $$f_{n+1} = \frac{2n+1}{x}\,f_n - f_{n-1}$$ рд╕реЗ рдкреНрд░рд╛рдкреНрдд рд╣реЛрддреА рд╣реИрдВред yv рдХреЗ рд▓рд┐рдП рдКрдкрд░ рдХреА рдУрд░ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕реНрдерд┐рд░ рд░рд╣рддреА рд╣реИ, рдкрд░рдВрддреБ jv рдХреЗ рд▓рд┐рдП рдЬрдм \(n > x\) рд╣реЛ рддреЛ рдпрд╣ рдЕрд╕реНрдерд┐рд░ рд╣реЛ рдЬрд╛рддреА рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╣рдо рдорд┐рд▓рд░ рдХреА рдиреАрдЪреЗ рдХреА рдУрд░ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдХрд╛ рдкреНрд░рдпреЛрдЧ рдХрд░рддреЗ рд╣реИрдВ: рдХрд┐рд╕реА рдЙрдЪреНрдЪ рдХреЛрдЯрд┐ рд╕реЗ f рдХреЛ 0 рдФрд░ 1 рд░рдЦрдХрд░ рд╢реБрд░реВ рдХрд░реЗрдВ, рдиреАрдЪреЗ рдХреА рдУрд░ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдХрд░реЗрдВ, рдлрд┐рд░ рд╣рд░ рдорд╛рди рдХреЛ рдЗрд╕ рдкреНрд░рдХрд╛рд░ рдкреБрдирдГ рдорд╛рдкреЗрдВ рдХрд┐ рдХреЛрдЯрд┐-0 рдХрд╛ рдкрдж \(\sin(x)/x\) рд╕реЗ рдореЗрд▓ рдЦрд╛рдПред рдЕрд╡рдХрд▓рдЬ \(j'_v = j_{v-1} - \frac{v+1}{x} j_v\) рд╕реЗ рдирд┐рдХрд╛рд▓реЗ рдЬрд╛рддреЗ рд╣реИрдВред

рдХреНрд╖реАрдг рд╣реЛрддреЗ рджреЛрд▓рдиреЛрдВ рдХреЗ рд░реВрдк рдореЗрдВ рдкрд╣рд▓реА рддрд░рд╣ рдХреЗ рдЧреЛрд▓рд╛рдХрд╛рд░ рдмреЗрд╕рд▓ рдлрд▓рдиреЛрдВ рдХрд╛ рдЧреНрд░рд╛рдл
рдкрд╣рд▓реА рддрд░рд╣ рдХреЗ рдЧреЛрд▓рд╛рдХрд╛рд░ рдмреЗрд╕рд▓ рдлрд▓рди jv(x) x рдмрдврд╝рдиреЗ рдкрд░ рджреЛрд▓рди рдХрд░рддреЗ рдФрд░ рдХреНрд╖реАрдг рд╣реЛрддреЗ рд╣реИрдВред

рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг (v = 0, x = 2)

\(j_0(2) = \sin(2)/2 = 0.4546487134\)ред \(y_0(2) = -\cos(2)/2 = 0.2080734183\)ред рдЪреВрдБрдХрд┐ \(j'_0 = -j_1\), рдЗрд╕рд▓рд┐рдП \(j'_0(2) = -0.4353977750\), рдФрд░ \(y'_0 = -y_1\) рд╕реЗ \(0.3506120043\) рдкреНрд░рд╛рдкреНрдд рд╣реЛрддрд╛ рд╣реИред

j0(x) = sin x / x рдХрд╛ рдЧреНрд░рд╛рдл рдЬрд┐рд╕рдореЗрдВ x=2 рдХреЗ рдкрд╛рд╕ рдПрдХ рдмрд┐рдВрджреБ рдЪрд┐рд╣реНрдирд┐рдд рд╣реИ
рд╢реВрдиреНрдп рдХреЛрдЯрд┐ рдХрд╛ рд╡рдХреНрд░ j0(x) = sin(x)/x, рд╣рд▓ рдХрд┐рдП рдЧрдП рдЙрджрд╛рд╣рд░рдг рдореЗрдВ x = 2 рдкрд░ рдореВрд▓реНрдпрд╛рдВрдХрд┐рддред

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

рдХреНрдпрд╛ рдпрд╣ рд╕рдореНрдорд┐рд╢реНрд░ (complex) x рдХрд╛ рд╕рдорд░реНрдерди рдХрд░рддрд╛ рд╣реИ? рдирд╣реАрдВред рдореВрд▓ рдкреГрд╖реНрда рд╕рдореНрдорд┐рд╢реНрд░ рддрд░реНрдХ рд╕реНрд╡реАрдХрд╛рд░ рдХрд░рддрд╛ рдерд╛; рдпрд╣ рдирдпрд╛ рд╕рдВрд╕реНрдХрд░рдг рд╕реНрдкрд╖реНрдЯрддрд╛ рдФрд░ рдЧрддрд┐ рдХреЗ рд▓рд┐рдП рдХреЗрд╡рд▓ рд╡рд╛рд╕реНрддрд╡рд┐рдХ \(x > 0\) рддрдХ рд╕реАрдорд┐рдд рд╣реИред

x = 0 рдкрд░ yv рдЕрдирдВрдд рдХреНрдпреЛрдВ рд╣реЛрддрд╛ рд╣реИ? рджреВрд╕рд░реЗ рдкреНрд░рдХрд╛рд░ рдХреЗ рдлрд▓рдиреЛрдВ рдХрд╛ рдореВрд▓рдмрд┐рдВрджреБ рдкрд░ рдПрдХ рдзреНрд░реБрд╡ (pole) рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП x рдХреЗ рд╢реВрдиреНрдп рдХреА рдУрд░ рдмрдврд╝рдиреЗ рдкрд░ рдЗрдирдХреЗ рдорд╛рди рдмрд┐рдирд╛ рдХрд┐рд╕реА рд╕реАрдорд╛ рдХреЗ рдмрдврд╝рддреЗ рдЬрд╛рддреЗ рд╣реИрдВред

рдпрд╣ рдХрд┐рддрдирд╛ рд╕рдЯреАрдХ рд╣реИ? рдЧрдгрдирд╛рдПрдБ рдбрдмрд▓ рдкреНрд░рд┐рд╕реАрдЬрд╝рди рдореЗрдВ рд╣реЛрддреА рд╣реИрдВ, рдЬрд┐рд╕рд╕реЗ рд▓рдЧрднрдЧ 15 рд╕рд╛рд░реНрдердХ рдЕрдВрдХ рдорд┐рд▓рддреЗ рд╣реИрдВ тАФ рд╕рд╛рдорд╛рдиреНрдп рдЗрдВрдЬреАрдирд┐рдпрд░рд┐рдВрдЧ рдФрд░ рднреМрддрд┐рдХреА рдХреЗ рдХрд╛рд░реНрдпреЛрдВ рдХреЗ рд▓рд┐рдП рдпрд╣ рдкрд░реНрдпрд╛рдкреНрдд рд╕реЗ рдЕрдзрд┐рдХ рд╣реИред

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