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

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

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

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

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

рдкрд╣рд▓реЗ рддрд░реНрдХ рдкрд░ рдЧрд╛рдорд╛ рдХрд╛ рдорд╛рди
1.772454
over 5 rows (5 finite)
рдмрдирд╛рдИ рдЧрдИ рдкрдВрдХреНрддрд┐рдпрд╛рдБ 5
рдкрд░рд┐рдорд┐рдд рдорд╛рди 5
рдиреНрдпреВрдирддрдо рдкрд░рд┐рдорд┐рдд рдорд╛рди 0.886227
рдЕрдзрд┐рдХрддрдо рдкрд░рд┐рдорд┐рдд рдорд╛рди 1.772454
a Gamma(a)
0.500000 1.772454
1.000000 1.000000
1.500000 0.886227
2.000000 1.000000
2.500000 1.329340

рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рд╣реИ?

рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди, рдЬрд┐рд╕реЗ \(\Gamma(a)\) рд▓рд┐рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЧрдгрд┐рдд рдХреЗ рд╕рдмрд╕реЗ рдорд╣рддреНрд╡рдкреВрд░реНрдг рд╡рд┐рд╢реЗрд╖ рдлрд╝рдВрдХреНрд╢рдиреЛрдВ рдореЗрдВ рд╕реЗ рдПрдХ рд╣реИред рдпрд╣ рджрд░рдЕрд╕рд▓ рдлрд╝реИрдХреНрдЯреЛрд░рд┐рдпрд▓ рдХрд╛ рд╡рд┐рд╢реНрд▓реЗрд╖рдгрд╛рддреНрдордХ рд╡рд┐рд╕реНрддрд╛рд░ (analytic continuation) рд╣реИ: рдХрд┐рд╕реА рднреА рдЕрдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХ \(n\) рдХреЗ рд▓рд┐рдП \(\Gamma(n+1) = n!\) рд╣реЛрддрд╛ рд╣реИред рдлрд╝реИрдХреНрдЯреЛрд░рд┐рдпрд▓ рдХреЗ рд╡рд┐рдкрд░реАрдд, рдЧрд╛рдорд╛ рд▓рдЧрднрдЧ рд╣рд░ рд╡рд╛рд╕реНрддрд╡рд┐рдХ (рдФрд░ рд╕рдореНрдорд┐рд╢реНрд░) рд╕рдВрдЦреНрдпрд╛ рдХреЗ рд▓рд┐рдП рдкрд░рд┐рднрд╛рд╖рд┐рдд рд╣реИ тАФ рднрд┐рдиреНрдиреЛрдВ рдФрд░ рдЛрдгрд╛рддреНрдордХ рд╕рдВрдЦреНрдпрд╛рдУрдВ рдХреЗ рд▓рд┐рдП рднреАред рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рд╕рдорд╛рди рдЕрдВрддрд░рд╛рд▓ рд╡рд╛рд▓реЗ рдорд╛рдиреЛрдВ рдХреА рдПрдХ рд╢реГрдВрдЦрд▓рд╛ рдкрд░ \(\Gamma(a)\) рдХреА рд╕рд╛рд░рдгреА рддреИрдпрд╛рд░ рдХрд░рддрд╛ рд╣реИ рдФрд░ рдкрд░рд┐рдгрд╛рдореА рд╡рдХреНрд░ рдХрд╛ рдЧреНрд░рд╛рдлрд╝ рдЦреАрдВрдЪрддрд╛ рд╣реИред рдпрд╣ рд╕рд╛рд░реНрд╡рднреМрдорд┐рдХ рдЧрдгрд┐рдд рд╣реИ тАФ рд╣рд░ рдЬрдЧрд╣ рд▓рд╛рдЧреВ, рдХрд┐рд╕реА рджреЗрд╢-рд╡рд┐рд╢реЗрд╖ рдирд┐рдпрдо рдХреЗ рдмрд┐рдирд╛ред

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

рддреАрди рдорд╛рди рджрд░реНрдЬ рдХрд░реЗрдВ: a рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди (рдкрд╣рд▓рд╛ рддрд░реНрдХ/argument), рд╕реНрдЯреЗрдк (рд╣рд░ рдкрдВрдХреНрддрд┐ рдореЗрдВ рдЬреБрдбрд╝рдиреЗ рд╡рд╛рд▓реА рдирд┐рдпрдд рд╡реГрджреНрдзрд┐), рдФрд░ рдкрдВрдХреНрддрд┐рдпреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ред рдкрдВрдХреНрддрд┐ \(k\) рдХреЗ рд▓рд┐рдП рддрд░реНрдХ $$a_k = \text{Initial }a + k \cdot \text{Step}, \quad k = 0,1,\dots,\text{Rows}-1$$ рд╣реЛрддрд╛ рд╣реИред рдЯреВрд▓ рд╣рд░ рддрд░реНрдХ рдкрд░ рдЧрд╛рдорд╛ рдХрд╛ рдорд╛рди рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИ, рдЗрди рдЬреЛрдбрд╝реЛрдВ рдХреЛ рдПрдХ рд╕рд╛рд░рдгреА рдореЗрдВ рджрд┐рдЦрд╛рддрд╛ рд╣реИ, рдФрд░ рдиреНрдпреВрдирддрдо рддрдерд╛ рдЕрдзрд┐рдХрддрдо рдкрд░рд┐рдорд┐рдд рдорд╛рди рдмрддрд╛рддрд╛ рд╣реИред \(a = 0, -1, -2, \dots\) рдкрд░ рдореМрдЬреВрдж рдзреНрд░реБрд╡реЛрдВ (poles) рдХреЛ рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд (undefined) рдХреЗ рд░реВрдк рдореЗрдВ рджрд░реНрд╢рд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред

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

рдкрд░рд┐рднрд╛рд╖рд┐рдд рд╕рдорд╛рдХрд▓рди рд╣реИ $$\Gamma(a) = \int_{0}^{\infty} t^{\,a - 1}\, e^{-t}\, dt$$ рдЬрд╣рд╛рдБ \(\operatorname{Re}(a) > 0\)ред рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╣рдо рд▓реИрдВрдХреНрдЬрд╝реЛрд╕ рд╕рдиреНрдирд┐рдХрдЯрди (Lanczos approximation) рдХрд╛ рдкреНрд░рдпреЛрдЧ рдХрд░рддреЗ рд╣реИрдВ, рдЬреЛ рд▓рдЧрднрдЧ 15 рдЕрдВрдХреЛрдВ рддрдХ рдХреА рд╕рдЯреАрдХрддрд╛ рджреЗрддрд╛ рд╣реИред \(a \le 0.5\) рдХреЗ рд▓рд┐рдП рдкрд░рд╛рд╡рд░реНрддрди рд╕реВрддреНрд░ (reflection formula) $$\Gamma(a) = \frac{\pi}{\sin(\pi a) \cdot \Gamma(1-a)}$$ рдЛрдгрд╛рддреНрдордХ рддрдерд╛ рдЫреЛрдЯреЗ рддрд░реНрдХреЛрдВ рдХреЛ рд╕рдБрднрд╛рд▓рддрд╛ рд╣реИ рдФрд░ рдЕрдкрд╕рд░рдг (divergence) рд╕реЗ рдмрдЪрд╛рддрд╛ рд╣реИред рдкреНрд░рдореБрдЦ рд╡рд┐рд╢реЗрд╖ рдорд╛рди: \(\Gamma(1/2) = \sqrt{\pi} = 1.77245\), \(\Gamma(1) = 1\), \(\Gamma(n+1) = n!\)ред

рдЕрдирдзрдирд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХреЛрдВ рдкрд░ рдзреНрд░реБрд╡реЛрдВ рд╡рд╛рд▓реЗ рдЧрд╛рдорд╛ рдлрд▓рди рдХрд╛ рдЧреНрд░рд╛рдл
рдЧрд╛рдорд╛ рдлрд▓рди рдХрд╛ рд╡рдХреНрд░, рдЬреЛ рдЗрд╕рдХреА рддреАрд╡реНрд░ рд╡реГрджреНрдзрд┐ рдФрд░ рд╢реВрдиреНрдп рд╡ рдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХреЛрдВ рдкрд░ рдКрд░реНрдзреНрд╡рд╛рдзрд░ рдЕрдирдВрддрд╕реНрдкрд░реНрд╢реА (рдзреНрд░реБрд╡) рджрд┐рдЦрд╛рддрд╛ рд╣реИред

рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг

рдорд╛рди рд▓реЗрдВ \(\text{initialA} = 0.5\), \(\text{step} = 0.5\), \(\text{count} = 5\), рддреЛ рддрд░реНрдХ рд╣реЛрдВрдЧреЗ \(0.5, 1.0, 1.5, 2.0, 2.5\)ред рдкрд░рд┐рдгрд╛рдо рд╣реИрдВ \(\Gamma(0.5) = 1.77245\) (\(= \sqrt{\pi}\)), \(\Gamma(1.0) = 1.0\), \(\Gamma(1.5) = 0.88623\), \(\Gamma(2.0) = 1.0\), рдФрд░ \(\Gamma(2.5) = 1.32934\)ред рдзрдирд╛рддреНрдордХ \(a\) рдХреЗ рд▓рд┐рдП рд╡рдХреНрд░ рдЕрдкрдиреЗ рдиреНрдпреВрдирддрдо (рд▓рдЧрднрдЧ \(0.8856\), \(a = 1.4616\) рдХреЗ рдЖрд╕рдкрд╛рд╕) рддрдХ рдЧрд┐рд░рддрд╛ рд╣реИ, рдлрд┐рд░ рджреЛрдмрд╛рд░рд╛ рдКрдкрд░ рдЙрдарддрд╛ рд╣реИред

рдЧрд╛рдорд╛ рд╕рдорд╛рдХрд▓ рдХреЛ рджрд░реНрд╢рд╛рддреЗ t^(a-1) e^(-t) рд╡рдХреНрд░ рдХреЗ рдиреАрдЪреЗ рдХрд╛ рдХреНрд╖реЗрддреНрд░рдлрд▓
рд╕рдорд╛рдХрд▓ рдкрд░рд┐рднрд╛рд╖рд╛ рдХрд╛ рдЬреНрдпрд╛рдорд┐рддреАрдп рдЕрд░реНрде: \(\Gamma(a)\), \(t^{a-1} \cdot e^{-t}\) рдХреЗ рдиреАрдЪреЗ рдЫрд╛рдпрд╛рдВрдХрд┐рдд рдХреНрд╖реЗрддреНрд░рдлрд▓ рд╣реИред

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

\(\Gamma(0)\) рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд рдХреНрдпреЛрдВ рд╣реИ? рдЕрдзрдирд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХ (\(0, -1, -2, \dots\)) рд╕рд░рд▓ рдзреНрд░реБрд╡ (simple poles) рд╣реИрдВ, рдЬрд╣рд╛рдБ рдЧрд╛рдорд╛ рдзрди рдпрд╛ рдЛрдг рдЕрдирдВрдд рдХреА рдУрд░ рдЕрдкрд╕рд░рд┐рдд рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╡рд╣рд╛рдБ рдХреЛрдИ рдкрд░рд┐рдорд┐рдд рдорд╛рди рдореМрдЬреВрдж рдирд╣реАрдВ рд╣реЛрддрд╛ред

рдХреНрдпрд╛ \(a\) рдЛрдгрд╛рддреНрдордХ рд╣реЛ рд╕рдХрддрд╛ рд╣реИ? рд╣рд╛рдБред рдЛрдгрд╛рддреНрдордХ рдЕрдкреВрд░реНрдгрд╛рдВрдХ (non-integer) рдорд╛рди рдорд╛рдиреНрдп рд╣реИрдВ; рд▓рдЧрд╛рддрд╛рд░ рдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХреЛрдВ рдХреЗ рдмреАрдЪ рдЗрди рдорд╛рдиреЛрдВ рдХрд╛ рдЪрд┐рд╣реНрди рдмрджрд▓рддрд╛ рд░рд╣рддрд╛ рд╣реИ рдФрд░ рдкрд░рд┐рдорд╛рдг рдмрдбрд╝рд╛ рд╣реЛрддрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬреИрд╕реЗ \(\Gamma(-0.5) = -3.5449\)ред

рдкрд░рд┐рдгрд╛рдо рдХрд┐рддрдирд╛ рд╕рдЯреАрдХ рд╣реИ? рд▓реИрдВрдХреНрдЬрд╝реЛрд╕ рд╕рдиреНрдирд┐рдХрдЯрди рд▓рдЧрднрдЧ 15 рд╕рд╛рд░реНрдердХ рдЕрдВрдХ рджреЗрддрд╛ рд╣реИ, рдЬреЛ рд▓рдЧрднрдЧ рд╕рднреА рд╡реНрдпрд╛рд╡рд╣рд╛рд░рд┐рдХ рдФрд░ рд╢реИрдХреНрд╖рд┐рдХ рдЙрдкрдпреЛрдЧреЛрдВ рдХреЗ рд▓рд┐рдП рдкрд░реНрдпрд╛рдкреНрдд рд╣реИред

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