рдкреЙрд▓реАрдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреНрдпрд╛ рд╣реИ?
рдХреНрд░рдо n рдХрд╛ рдкреЙрд▓реАрдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди, рдЬрд┐рд╕реЗ \(\psi^{(n)}(a)\) рд▓рд┐рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреЗ рдкреНрд░рд╛рдХреГрддрд┐рдХ рд▓рдШреБрдЧрдгрдХ (рдиреИрдЪреБрд░рд▓ рд▓реЙрдЧрд░рд┐рджрдо) рдХрд╛ (n+1)-рд╡рд╛рдБ рдЕрд╡рдХрд▓рдЬ (рдбреЗрд░рд┐рд╡реЗрдЯрд┐рд╡) рд╣реЛрддрд╛ рд╣реИред рдпрд╣рд╛рдБ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд┐рдП рдЧрдП рдЗрдВрдбреЗрдХреНрд╕рд┐рдВрдЧ рдХреЗ рдЕрдиреБрд╕рд╛рд░, \(n = -1\) рд╕реНрд╡рдпрдВ рд▓реЙрдЧ-рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди \(\ln\Gamma(a)\) рджреЗрддрд╛ рд╣реИ, \(n = 0\) рдбрд╛рдЗрдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди \(\psi(a)\) рджреЗрддрд╛ рд╣реИ, \(n = 1\) рдЯреНрд░рд╛рдЗрдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди, \(n = 2\) рдЯреЗрдЯреНрд░рд╛рдЧрд╛рдорд╛, рдФрд░ рдЗрд╕реА рддрд░рд╣ рдЖрдЧреЗред рдпреЗ рд╡рд┐рд╢реЗрд╖ рдлрд╝рдВрдХреНрд╢рди рдкреНрд░рд╛рдпрд┐рдХрддрд╛ (рдкреНрд░реЛрдмреЗрдмрд┐рд▓рд┐рдЯреА), рд╕рд╛рдВрдЦреНрдпрд┐рдХреА (рдЧрд╛рдорд╛ рдФрд░ рдмреАрдЯрд╛ рд╡рд┐рддрд░рдгреЛрдВ рдХреЗ рд▓рд┐рдП рдореИрдХреНрд╕рд┐рдордо-рд▓рд╛рдЗрдХрд▓реАрд╣реБрдб рдЕрдиреБрдорд╛рди), рд╕рдВрдЦреНрдпрд╛ рд╕рд┐рджреНрдзрд╛рдВрдд рдФрд░ рдПрд╕рд┐рдореНрдкреНрдЯреЛрдЯрд┐рдХ рд╡рд┐рд╢реНрд▓реЗрд╖рдг рдореЗрдВ рд╣рд░ рдЬрдЧрд╣ рджрд┐рдЦрд╛рдИ рджреЗрддреЗ рд╣реИрдВред рдпрд╣ рд╢реБрджреНрдз рдЧрдгрд┐рдд рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рдПрдХ рд╕рдорд╛рди рд░реВрдк рд╕реЗ рд▓рд╛рдЧреВ рд╣реЛрддрд╛ рд╣реИред
$$\psi^{(n)}(a) = \frac{d^{n+1}}{da^{n+1}} \ln\Gamma\!\left(a\right), \quad n = \text{order}$$
рдЗрд╕ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рдбреНрд░реЙрдкрдбрд╛рдЙрди рд╕реЗ рдЕрд╡рдХрд▓рдЬ рдХреНрд░рдо \(n\) рдЪреБрдиреЗрдВ (-1, 0, 1, 2, 3, рдпрд╛ 4) рдФрд░ рддрд░реНрдХ (рдЖрд░реНрдЧреНрдпреБрдореЗрдВрдЯ) \(a\) рдХреЛ рдПрдХ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рдХреЗ рд░реВрдк рдореЗрдВ рдЯрд╛рдЗрдк рдХрд░реЗрдВред calculate рджрдмрд╛рдПрдБ рдФрд░ рд▓рдЧрднрдЧ 12 рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ рддрдХ \(\psi^{(n)}(a)\) рдХрд╛ рдорд╛рди рдкрд╛рдПрдБред рдпрд╣ рдЧрдгрдирд╛ \(a > 0\) рдХреЗ рд▓рд┐рдП рдкреВрд░реА рддрд░рд╣ рдорд╛рдиреНрдп рд╣реИ; рдЬрдм \(a\), 0 рд╕реЗ рдХрдо рдпрд╛ рдЙрд╕рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛ рддреЛ рдЗрди рдлрд╝рдВрдХреНрд╢рдиреЛрдВ рдХреЗ рд▓рд┐рдП рд░рд┐рдлреНрд▓реЗрдХреНрд╢рди рдХреА рдЬрд╝рд░реВрд░рдд рдкрдбрд╝рддреА рд╣реИ, рдФрд░ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдЕрдЛрдгрд╛рддреНрдордХ рди рд╣реЛрдиреЗ рд╡рд╛рд▓реЗ рдкреВрд░реНрдгрд╛рдВрдХреЛрдВ (рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреЗ рдкреЛрд▓/рдзреНрд░реБрд╡реЛрдВ) рдХреЛ рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд рдмрддрд╛рддрд╛ рд╣реИред
рд╕реВрддреНрд░ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛
\(n = -1\) рдХреЗ рд▓рд┐рдП рдорд╛рди, \(\ln\Gamma\) рд╣реЗрддреБ рд▓реИрдВрдХреНрдЬрд╝реЛрд╕ рд╕рдиреНрдирд┐рдХрдЯрди (Lanczos approximation) рд╕реЗ рдирд┐рдХрд╛рд▓рд╛ рдЬрд╛рддрд╛ рд╣реИред \(n = 0\) рдХреЗ рд▓рд┐рдП рдбрд╛рдЗрдЧрд╛рдорд╛ рдХреЛ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕реВрддреНрд░ \(\psi(x) = \psi(x+1) - \frac{1}{x}\) рд╕реЗ рдЬреНрдЮрд╛рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ рддрд╛рдХрд┐ рддрд░реНрдХ рдХреЛ 6 рд╕реЗ рдКрдкрд░ рд▓реЗ рдЬрд╛рдпрд╛ рдЬрд╛ рд╕рдХреЗ, рдЙрд╕рдХреЗ рдмрд╛рдж рдПрд╕рд┐рдореНрдкреНрдЯреЛрдЯрд┐рдХ рд╢реНрд░реЗрдгреА \(\psi(x) \sim \ln x - \frac{1}{2x} - \frac{1}{12 x^2} + \cdots\) рд▓рдЧрд╛рдИ рдЬрд╛рддреА рд╣реИред \(n \geq 1\) рдХреЗ рд▓рд┐рдП рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕реВрддреНрд░ \(\psi^{(n)}(x) = \psi^{(n)}(x+1) + \frac{(-1)^{n+1}\, n!}{x^{n+1}}\) рддрд░реНрдХ рдХреЛ 10 рд╕реЗ рдЖрдЧреЗ рдЦрд┐рд╕рдХрд╛рддрд╛ рд╣реИ, рдлрд┐рд░ рдмрд░реНрдиреВрд▓реА-рд╕рдВрдЦреНрдпрд╛ рдЖрдзрд╛рд░рд┐рдд рдПрд╕рд┐рдореНрдкреНрдЯреЛрдЯрд┐рдХ рд╢реНрд░реЗрдгреА рд╕реЗ рдЧрдгрдирд╛ рдкреВрд░реА рд╣реЛрддреА рд╣реИред
рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг
рдорд╛рди рд▓реАрдЬрд┐рдП \(n = 1\), \(a = 5\) (рдЕрд░реНрдерд╛рддреН 5 рдкрд░ рдЯреНрд░рд╛рдЗрдЧрд╛рдорд╛)ред рд╕рд░реНрд╡рд╕рдорд┐рдХрд╛ \(\psi^{(1)}(m) = \frac{\pi^2}{6} - \sum_{k=1}^{m-1} \frac{1}{k^2}\) рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП, рд╣рдореЗрдВ рдорд┐рд▓рддрд╛ рд╣реИ $$1.6449340668 - \left(1 + \frac{1}{4} + \frac{1}{9} + \frac{1}{16}\right) = 1.6449340668 - 1.4236111111 = 0.2213229557$$ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рд▓рдЧрднрдЧ \(0.221322955737\) рд▓реМрдЯрд╛рддрд╛ рд╣реИред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди
\(a = 0\) рдпрд╛ рдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХ рдкрд░ рдорд╛рди рдЕрдкрд░рд┐рднрд╛рд╖рд┐рдд рдХреНрдпреЛрдВ рд╣реЛрддрд╛ рд╣реИ? рдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреЗ \(0, -1, -2, \ldots\) рдкрд░ рдкреЛрд▓ (рдзреНрд░реБрд╡) рд╣реЛрддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП рд╡рд╣рд╛рдБ \(\ln\Gamma\) рдФрд░ рд╣рд░ рдкреЙрд▓реАрдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдЕрдкрд╕рд░рдг (diverge) рдХрд░рддреЗ рд╣реИрдВред
рдбрд╛рдЗрдЧрд╛рдорд╛ рдФрд░ рдЯреНрд░рд╛рдЗрдЧрд╛рдорд╛ рдореЗрдВ рдХреНрдпрд╛ рдЕрдВрддрд░ рд╣реИ? рдбрд╛рдЗрдЧрд╛рдорд╛ (\(n = 0\)), \(\ln\Gamma\) рдХрд╛ рдкрд╣рд▓рд╛ рдЕрд╡рдХрд▓рдЬ рд╣реИ; рдЯреНрд░рд╛рдЗрдЧрд╛рдорд╛ (\(n = 1\)) рджреВрд╕рд░рд╛ рдЕрд╡рдХрд▓рдЬ рд╣реИ, рдЬреЛ рдбрд╛рдЗрдЧрд╛рдорд╛ рдХреЗ рдЕрд╡рдХрд▓рдЬ рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛрддрд╛ рд╣реИред
рдкрд░рд┐рдгрд╛рдо рдХрд┐рддрдирд╛ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИ? рд╢рд┐рдлреНрдЯ-рдПрдВрдб-рдПрд╕рд┐рдореНрдкреНрдЯреЛрдЯрд┐рдХ рд╡рд┐рдзрд┐ рдХреЗ рд╕рд╛рде рдбрдмрд▓-рдкреНрд░рд┐рд╕рд┐рдЬрд╝рди рдЕрдВрдХрдЧрдгрд┐рдд, \(a > 0\) рдХреЗ рд▓рд┐рдП рд▓рдЧрднрдЧ 12-15 рд╕рд╣реА рд╕рд╛рд░реНрдердХ рдЕрдВрдХ рджреЗрддрд╛ рд╣реИред