рдПрдпрд░реА рдлрд╝рдВрдХреНрд╢рди рдХреНрдпрд╛ рд╣реИрдВ?
рдПрдпрд░реА рдлрд╝рдВрдХреНрд╢рди 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) рдХреЗ рдорд╛рди рджреЗрдЦреЗрдВред рд╢реБрд░реБрдЖрдд рдХреЗ рд▓рд┐рдП рдбрд┐рдлрд╝реЙрд▓реНрдЯ рдорд╛рди 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 рд░реВрдкреЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред