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

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

n = -1, 0, 1, 2, 3, 4
argument (real number, valid for a > 0)

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

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

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

Polygamma value ╬иn(a)
0.221322955737
at n = 1, a = 5.0
рдЕрд╡рдХрд▓рдЬ рдХреНрд░рдо n 1
рддрд░реНрдХ a 5.0
рдлрд╝рдВрдХреНрд╢рди Trigamma ╬и1(a)

рдкреЙрд▓реАрдЧрд╛рдорд╛ рдлрд╝рдВрдХреНрд╢рди рдХреНрдпрд╛ рд╣реИ?

рдХреНрд░рдо 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 рдХреЛ ln ╬У(a) рдХрд╛ n+1 рдмрд╛рд░ рдФрд░ рдЕрд╡рдХрд▓рди рдХрд░рдХреЗ рдкреНрд░рд╛рдкреНрдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред

рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг

рдорд╛рди рд▓реАрдЬрд┐рдП \(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 рд╕рд╣реА рд╕рд╛рд░реНрдердХ рдЕрдВрдХ рджреЗрддрд╛ рд╣реИред

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