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

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

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

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

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

рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░ рдХрд╛ рдПрдпрд░реА рдлрд╝рдВрдХреНрд╢рди
0.1352924163
Ai(x)
Ai(x) (рдкреНрд░рдердо рдкреНрд░рдХрд╛рд░) 0.1352924163
Bi(x) (рджреНрд╡рд┐рддреАрдп рдкреНрд░рдХрд╛рд░) 1.207423595

рдПрдпрд░реА рдлрд╝рдВрдХреНрд╢рди рдХреНрдпрд╛ рд╣реИрдВ?

рдПрдпрд░реА рдлрд╝рдВрдХреНрд╢рди Ai(x) рдФрд░ Bi(x), рдПрдпрд░реА рдЕрд╡рдХрд▓ рд╕рдореАрдХрд░рдг \(y'' = x\cdot y\) (рдпрд╛рдиреА \(y'' - x\cdot y = 0\)) рдХреЗ рджреЛ рд░реИрдЦрд┐рдХ рд░реВрдк рд╕реЗ рд╕реНрд╡рддрдВрддреНрд░ рд╣рд▓ рд╣реИрдВред рдпреЗ рднреМрддрд┐рдХреА рдФрд░ рдЕрдиреБрдкреНрд░рдпреБрдХреНрдд рдЧрдгрд┐рдд рдореЗрдВ рдЬрдЧрд╣-рдЬрдЧрд╣ рджрд┐рдЦрд╛рдИ рджреЗрддреЗ рд╣реИрдВ: рдХреНрд╡рд╛рдВрдЯрдо рдпрд╛рдВрддреНрд░рд┐рдХреА рдореЗрдВ classical turning points рдХреЗ рдкрд╛рд╕ (WKB connection рд╕рдорд╕реНрдпрд╛ рдореЗрдВ), рдкреНрд░рдХрд╛рд╢рд┐рдХреА рдореЗрдВ caustics рдФрд░ рдЗрдВрджреНрд░рдзрдиреБрд╖ рдХреЗ рд╡рд░реНрдгрди рдореЗрдВ, рдФрд░ asymptotic рд╡рд┐рд╢реНрд▓реЗрд╖рдг рдореЗрдВред Ai(x) рд╡рд╣ рд╣рд▓ рд╣реИ рдЬреЛ x рдХреЗ рдмрдбрд╝реЗ рдФрд░ рдзрдирд╛рддреНрдордХ рд╣реЛрдиреЗ рдкрд░ рдХреНрд╖рдп рд╣реЛрддрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬрдмрдХрд┐ Bi(x) рдЙрд╕реА рд╕реАрдорд╛ рдореЗрдВ рдЪрд░рдШрд╛рддрд╛рдВрдХреА (exponential) рд░реВрдк рд╕реЗ рдмрдврд╝рддрд╛ рд╣реИред рдЛрдгрд╛рддреНрдордХ x рдХреЗ рд▓рд┐рдП рджреЛрдиреЛрдВ рджреЛрд▓рди (oscillate) рдХрд░рддреЗ рд╣реИрдВ рдФрд░ \(|x|^{-1/4}\) рдХреА рджрд░ рд╕реЗ рдзреАрд░реЗ-рдзреАрд░реЗ рдШрдЯрддреЗ рд╣реИрдВред

x рдХреЗ рд╕рд╛рдкреЗрдХреНрд╖ рдПрдпрд░реА рдлрд▓рди Ai(x) рдФрд░ Bi(x) рдХрд╛ рдЧреНрд░рд╛рдл
рдПрдпрд░реА рдлрд▓рди Ai(x) рдФрд░ Bi(x): рдЛрдгрд╛рддреНрдордХ x рдХреЗ рд▓рд┐рдП рджреЛрд▓рдирд╢реАрд▓, рдзрдирд╛рддреНрдордХ x рдХреЗ рд▓рд┐рдП Ai рдШрдЯрддрд╛ рд╣реИ рдФрд░ Bi рдмрдврд╝рддрд╛ рд╣реИред

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

x рдХрд╛ рдХреЛрдИ рднреА рдкрд░рд┐рдорд┐рдд рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдорд╛рди (рдзрдирд╛рддреНрдордХ, рдЛрдгрд╛рддреНрдордХ рдпрд╛ рд╢реВрдиреНрдп) рджрд░реНрдЬ рдХрд░реЗрдВ рдФрд░ Ai(x) рддрдерд╛ Bi(x) рдХреЗ рдорд╛рди рджреЗрдЦреЗрдВред рд╢реБрд░реБрдЖрдд рдХреЗ рд▓рд┐рдП рдбрд┐рдлрд╝реЙрд▓реНрдЯ рдорд╛рди x = 1.0 рджрд┐рдпрд╛ рдЧрдпрд╛ рд╣реИред рдпрд╣рд╛рдБ рдХреЛрдИ рдЗрдХрд╛рдИ рдирд╣реАрдВ рд╣реИ тАФ x рдПрдХ рд╢реБрджреНрдз, рд╡рд┐рдорд╛рд░рд╣рд┐рдд (dimensionless) рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реИред

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

рдЬрдм |x| рдордзреНрдпрдо рд╕реНрддрд░ рдХрд╛ рд╣реЛ, рддреЛ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рд╣рд░ рдЬрдЧрд╣ рдЕрднрд┐рд╕рд░рдг рдХрд░рдиреЗ рд╡рд╛рд▓реА рдШрд╛рдд рд╢реНрд░реЗрдгреА (power series) рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред рджреЛ рд╢реНрд░реЗрдгрд┐рдпрд╛рдБ f(x) рдФрд░ g(x) рдХреЛ рдПрдХ рд╕реНрдерд┐рд░ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ (recurrence) рдХреЗ рдЬрд╝рд░рд┐рдП рдЬреЛрдбрд╝рд╛ рдЬрд╛рддрд╛ рд╣реИ:

$$f(x)=\sum_{k\ge0}\frac{3^k(1/3)_k}{(3k)!}x^{3k},\quad g(x)=\sum_{k\ge0}\frac{3^k(2/3)_k}{(3k+1)!}x^{3k+1}$$

f рдХреЗ рд▓рд┐рдП, \(\text{term}_k = \text{term}_{k-1} \times x^3 / ((3k-1)(3k))\), рдЬрд┐рд╕рдХреА рд╢реБрд░реБрдЖрдд 1 рд╕реЗ рд╣реЛрддреА рд╣реИ; g рдХреЗ рд▓рд┐рдП, \(\text{term}_k = \text{term}_{k-1} \times x^3 / ((3k)(3k+1))\), рдЬрд┐рд╕рдХреА рд╢реБрд░реБрдЖрдд x рд╕реЗ рд╣реЛрддреА рд╣реИред рдЗрд╕рдХреЗ рдмрд╛рдж

$$\text{Ai}(x) = c_1 f(x) - c_2 g(x), \quad \text{Bi}(x) = \sqrt{3}\,\bigl(c_1 f(x) + c_2 g(x)\bigr)$$

рдЬрд╣рд╛рдБ \(c_1 = \text{Ai}(0) = 0.3550280539\) рдФрд░ \(c_2 = -\text{Ai}'(0) = 0.2588194038\) рд╣реИрдВред рдЬрдм |x| рд▓рдЧрднрдЧ 8 рд╕реЗ рдЬрд╝реНрдпрд╛рджрд╛ рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ, рддреЛ cancellation error рд╕реЗ рдмрдЪрдиреЗ рдХреЗ рд▓рд┐рдП рдХреЛрдб asymptotic рд╡рд┐рд╕реНрддрд╛рд░ рдкрд░ рд╕реНрд╡рд┐рдЪ рдХрд░ рджреЗрддрд╛ рд╣реИред

рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг (x = 1)

\(f(1) \approx 1.1722994\) рдФрд░ \(g(1) \approx 1.0853395\)ред рдЗрд╕рд▓рд┐рдП $$\text{Ai}(1) = 0.3550280539\times1.1722994 - 0.2588194038\times1.0853395 \approx 0.1352924,$$ рдФрд░ $$\text{Bi}(1) = \sqrt{3}\times(0.4161680 + 0.2808727) \approx 1.2074236.$$ рдпреЗ рдорд╛рди рдорд╛рдирдХ рд╕рдВрджрд░реНрдн рдорд╛рдиреЛрдВ рд╕реЗ рдореЗрд▓ рдЦрд╛рддреЗ рд╣реИрдВред

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

Ai(0) рдФрд░ Bi(0) рдХреНрдпрд╛ рд╣реИрдВ? \(\text{Ai}(0) = 0.3550280539\) рдФрд░ \(\text{Bi}(0) = \sqrt{3}\times\text{Ai}(0) = 0.6149266274\) тАФ рдпреЗ рд╕рдЯреАрдХ рдмрдВрдж-рд░реВрдк (closed form) рдорд╛рди рд╣реИрдВред

Bi(x) рдЗрддрдирд╛ рддреЗрдЬрд╝реА рд╕реЗ рдХреНрдпреЛрдВ рдмрдврд╝рддрд╛ рд╣реИ? рдмрдбрд╝реЗ рдзрдирд╛рддреНрдордХ x рдХреЗ рд▓рд┐рдП Bi(x), \(\exp((2/3)x^{3/2})\) рдХреА рджрд░ рд╕реЗ рдмрдврд╝рддрд╛ рд╣реИ рдФрд░ \(x \approx 100\) рдХреЗ рдЖрд╕рдкрд╛рд╕ double precision рд╕рдВрдЦреНрдпрд╛ рдХреА рд╕реАрдорд╛ рдкрд╛рд░ рдХрд░ (overflow) рдЬрд╛рддрд╛ рд╣реИ; рдпрд╣ рдЕрдкреЗрдХреНрд╖рд┐рдд рд╡реНрдпрд╡рд╣рд╛рд░ рд╣реИ, рдХреЛрдИ рддреНрд░реБрдЯрд┐ рдирд╣реАрдВред

рдХреНрдпрд╛ рдореИрдВ рдЛрдгрд╛рддреНрдордХ x рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░ рд╕рдХрддрд╛ рд╣реВрдБ? рд╣рд╛рдБред рдмрдбрд╝реЗ рдЛрдгрд╛рддреНрдордХ x рдХреЗ рд▓рд┐рдП рдпреЗ рдлрд╝рдВрдХреНрд╢рди рджреЛрд▓рди рдХрд░рддреЗ рд╣реИрдВ, рдФрд░ рд╕рдЯреАрдХрддрд╛ рдХреЗ рд▓рд┐рдП рдХреИрд▓рдХреБрд▓реЗрдЯрд░ oscillatory asymptotic рд░реВрдкреЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред

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