рд▓реЙрдЧ рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреНрдпрд╛ рд╣реИ?
рд▓реЙрдЧ рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди, рдЬрд┐рд╕реЗ \(\ln\Gamma(a)\) рдпрд╛ \(\ln(\Gamma(a))\) рд▓рд┐рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИ, рджрд░рдЕрд╕рд▓ рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди \(\Gamma(a)\) рдХрд╛ рдкреНрд░рд╛рдХреГрддрд┐рдХ рд▓рдШреБрдЧрдгрдХ (рдиреИрдЪреБрд░рд▓ рд▓реЙрдЧ) рд╣реЛрддрд╛ рд╣реИред рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдлрд╝реИрдХреНрдЯреЛрд░рд┐рдпрд▓ рдХреЛ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдФрд░ рд╕рдореНрдорд┐рд╢реНрд░ (complex) рд╕рдВрдЦреНрдпрд╛рдУрдВ рддрдХ рд╡рд┐рд╕реНрддреГрдд рдХрд░рддрд╛ рд╣реИ: рдХрд┐рд╕реА рдзрди рдкреВрд░реНрдгрд╛рдВрдХ \(n\) рдХреЗ рд▓рд┐рдП \(\Gamma(n) = (n-1)!\) рд╣реЛрддрд╛ рд╣реИред рдЪреВрдБрдХрд┐ \(\Gamma(a)\) рдмрд╣реБрдд рддреЗрдЬрд╝реА рд╕реЗ рдмрдврд╝рддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╡реИрдЬреНрдЮрд╛рдирд┐рдХ рдФрд░ рд╕рд╛рдВрдЦреНрдпрд┐рдХреАрд╡рд┐рджреН рд▓рдЧрднрдЧ рд╣рдореЗрд╢рд╛ рдЗрд╕рдХреЗ рд▓рдШреБрдЧрдгрдХ рдХреЗ рд╕рд╛рде рд╣реА рдХрд╛рдо рдХрд░рддреЗ рд╣реИрдВ рддрд╛рдХрд┐ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рдУрд╡рд░рдлрд╝реНрд▓реЛ рд╕реЗ рдмрдЪрд╛ рдЬрд╛ рд╕рдХреЗред \(\ln\Gamma\) рдкреНрд░рд╛рдпрд┐рдХрддрд╛ (рдмреАрдЯрд╛ рдФрд░ рдЧрд╛рдорд╛ рдмрдВрдЯрди, рдХрд╛рдИ-рд╡рд░реНрдЧ рдкрд░реАрдХреНрд╖рдг), рд╕рдВрдпреЛрдЬрд┐рдХреА (combinatorics) рдФрд░ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╡рд┐рд╢реНрд▓реЗрд╖рдг рдореЗрдВ рд╣рд░ рдЬрдЧрд╣ рджрд┐рдЦрд╛рдИ рджреЗрддрд╛ рд╣реИред
рдЗрд╕ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рддрд░реНрдХ (argument) \(a\) рдбрд╛рд▓реЗрдВ тАФ рдпрд╣ рдХреЛрдИ рднреА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реЛ рд╕рдХрддреА рд╣реИ тАФ рдФрд░ 'рдХреИрд▓рдХреБрд▓реЗрдЯ' рджрдмрд╛рдПрдБред рдпрд╣ рдЯреВрд▓ \(\ln\Gamma(a)\) рдХреЛ рдкреНрд░рд╛рдХреГрддрд┐рдХ рд▓реЙрдЧ рд░реВрдк рдореЗрдВ рджреЗрддрд╛ рд╣реИ, рд╕рд╛рде рд╣реА рдЖрдзрд╛рд░-10 рд╡рд╛рд▓рд╛ рд░реВрдк (Log10 Gamma) рднреАред рдЬрдм \(a\), 0 рд╕реЗ рдмрдбрд╝рд╛ рд╣реЛ рддреЛ рдкрд░рд┐рдгрд╛рдо рдПрдХ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реЛрддреА рд╣реИред \(a = 0, -1, -2, -3, \ldots\) рдкрд░ рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдореЗрдВ рдзреНрд░реБрд╡ (poles) рд╣реЛрддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП рд╡рд╣рд╛рдБ \(\ln\Gamma\) рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд рд░рд╣рддрд╛ рд╣реИ рдФрд░ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдпрд╣ рдмрддрд╛ рджреЗрддрд╛ рд╣реИред рдЛрдгрд╛рддреНрдордХ рдЧреИрд░-рдкреВрд░реНрдгрд╛рдВрдХ \(a\) рдХреЗ рд▓рд┐рдП \(\Gamma\) рдЛрдгрд╛рддреНрдордХ рд╣реЛ рд╕рдХрддрд╛ рд╣реИ; рдРрд╕реЗ рдорд╛рдорд▓реЛрдВ рдореЗрдВ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдкрд░рд╛рд╡рд░реНрддрди рд╕реВрддреНрд░ (reflection formula) рдХреА рдорджрдж рд╕реЗ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдореБрдЦреНрдп рдорд╛рди (principal value) рджреЗрддрд╛ рд╣реИред
рд╕реВрддреНрд░ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛
рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ Lanczos рдЕрдиреБрдорд╛рди (g = 7, рдиреМ рдЧреБрдгрд╛рдВрдХ) рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИ, рдЬреЛ рдПрдХ рддреЗрдЬрд╝ рдФрд░ рдмреЗрд╣рдж рд╕рдЯреАрдХ рдмрдВрдж-рд░реВрдк (closed-form) рд╡рд┐рдзрд┐ рд╣реИред \(x = a - 1\) рдФрд░ \(t = x + 7.5\) рд▓реЗрдХрд░, рдпрд╣ Lanczos рдЧреБрдгрд╛рдВрдХреЛрдВ рдХрд╛ рднрд╛рд░рд┐рдд рдпреЛрдЧ \(A\) рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИ рдФрд░ рдлрд┐рд░
$$\ln\Gamma(a) = 0.5\,\ln(2\pi) + \left(x + 0.5\right)\ln(t) - t + \ln(A)$$рдХрд╛ рдорд╛рди рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИред рдЬрдм \(a\), 0.5 рд╕реЗ рдХрдо рдпрд╛ рдЙрд╕рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛ, рддреЛ рдкрд░рд╛рд╡рд░реНрддрди рд╕реВрддреНрд░
$$\Gamma(a)\,\Gamma(1-a) = \frac{\pi}{\sin(\pi a)}$$рд╕рдорд╕реНрдпрд╛ рдХреЛ рдРрд╕реЗ рдорд╛рди рдкрд░ рд▓реЗ рдЬрд╛рддрд╛ рд╣реИ рдЬрд╣рд╛рдБ Lanczos рд╢реНрд░реЗрдгреА рдЕрдЪреНрдЫреА рддрд░рд╣ рдЕрднрд┐рд╕рд░рд┐рдд (converge) рд╣реЛрддреА рд╣реИред
рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг
рдорд╛рди рд▓реАрдЬрд┐рдП \(a = 3.5\):
$$\Gamma(3.5) = 2.5 \times 1.5 \times 0.5 \times \sqrt{\pi} = 1.875 \times 1.7724538509 = 3.32335097045$$рдЗрд╕рд▓рд┐рдП
$$\ln\Gamma(3.5) = \ln(3.32335097045) = 1.20097360234707$$рдПрдХ рдФрд░ рдЬрд╛рдБрдЪ рдХреЗ рддреМрд░ рдкрд░, \(a = 5\) рдХреЗ рд▓рд┐рдП \(\Gamma(5) = 4! = 24\) рд╣реЛрддрд╛ рд╣реИ, рддреЛ
$$\ln\Gamma(5) = \ln(24) = 3.17805383034795$$рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рд╕рд╡рд╛рд▓
рдЧрд╛рдорд╛ рдХреЗ рдмрдЬрд╛рдп рдЧрд╛рдорд╛ рдХрд╛ рд▓реЙрдЧ рдХреНрдпреЛрдВ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд░реЗрдВ? рд╕рд╛рдорд╛рдиреНрдп рдлрд╝реНрд▓реЛрдЯрд┐рдВрдЧ-рдкреЙрдЗрдВрдЯ рд╕рдВрдЦреНрдпрд╛рдУрдВ рдореЗрдВ \(\Gamma(a)\) рдордзреНрдпрдо рдорд╛рдиреЛрдВ рдкрд░ рд╣реА рдУрд╡рд░рдлрд╝реНрд▓реЛ рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ (\(\Gamma(171)\) рд╣реА рдбрдмрд▓ рд░реЗрдВрдЬ рд╕реЗ рдмрд╛рд╣рд░ рдирд┐рдХрд▓ рдЬрд╛рддрд╛ рд╣реИ), рдЬрдмрдХрд┐ \(\ln\Gamma\) рд╕рдВрднрд╛рд▓рдиреЗ рдпреЛрдЧреНрдп рдмрдирд╛ рд░рд╣рддрд╛ рд╣реИред рдЗрд╕реАрд▓рд┐рдП рд▓реЙрдЧ рд░реВрдк рд╡реНрдпрд╛рд╡рд╣рд╛рд░рд┐рдХ рд╡рд┐рдХрд▓реНрдк рд╣реИред
\(\ln\Gamma(1)\) рдФрд░ \(\ln\Gamma(2)\) рдХреНрдпрд╛ рд╣реЛрддреЗ рд╣реИрдВ? рджреЛрдиреЛрдВ 0 рд╣реЛрддреЗ рд╣реИрдВ, рдХреНрдпреЛрдВрдХрд┐ \(\Gamma(1) = \Gamma(2) = 1\) рдФрд░ \(\ln(1) = 0\)ред рдЗрд╕ рдлрд╝рдВрдХреНрд╢рди рдХрд╛ рдиреНрдпреВрдирддрдо рдорд╛рди рд▓рдЧрднрдЧ \(a = 1.4616\) рдХреЗ рдкрд╛рд╕ рдЖрддрд╛ рд╣реИред
рдХреНрдпрд╛ рдкрд░рд┐рдгрд╛рдо рдмрд┐рд▓реНрдХреБрд▓ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИ? Lanczos рдЕрдиреБрдорд╛рди рд╕рд╛рдорд╛рдиреНрдп рддрд░реНрдХреЛрдВ рдХреЗ рд▓рд┐рдП рд▓рдЧрднрдЧ 15 рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ рддрдХ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рдбрдмрд▓-рдкреНрд░рд┐рд╕рд┐рдЬрд╝рди рдлрд╝реНрд▓реЛрдЯрд┐рдВрдЧ рдкреЙрдЗрдВрдЯ рдХреЗ рдмрд░рд╛рдмрд░ рд╣реИред