On Bernoulli polynomials. Extract from a correspondence between Mr.~Sonin and Mr.~Hermite. (Q1521847)

In dem ersten Briefe Sonin's wird die \(i^{\text{te}}\) Bernoulli'sche Function \(\varphi_i (x)\) durch die Gleichungen \(\varphi_i (x+1) -\varphi_i (x) = ix^{i-1},\) \(\varphi_i (0) = 0\) definirt und in der Form \(\varphi_i (x) = p^{i-1} \left\{ \sum_{k=0}^{p-1} \varphi_i \left( \frac {x+k}{p} \right) -A_i (p) \right\}\) dargestellt, wo \(A_i (p) = \sum_{k=0}^{p-1} \varphi_i \left( \frac {k}{p} \right) ,\) \(A_{2n+1} (p) = 0\) und \(A_ {2n} (p) = \sum_{k=0}^{p-1} \varphi_{2n} \left( \frac {x+k}{p} \right) - p^{1-2n} \varphi_{2n} (x) = p (1-p^{-2n}) (-1)^n B_n,\) \((-1)^n B_n = \int_0^1 \varphi_{2n} (x) dx\) ist. Eine andere Darstellung ist \[ \varphi_i (x) = x^i - \tfrac{i}{2}\, x^{i-1} + {i \choose 2} B_1 x^{i-2} - {i \choose 4} B_2 x^{i-4} + {i \choose 6} B_3 x^{i-6} - \cdots; \] das letzte Glied der rechten Seite ist \[ -(-1)^{\frac {i-1}{2}} {i\choose i-1} B_{\frac {i-1}{2}}\;x \quad \text{oder} \quad (-1)^n {2n \choose 2} B_{n-1} x^2, \] je nachdem \(i\) eine ungerade oder gerade Zahl ist. Für eine ganze Function \(f (z)\) vom \(n^{\text{ten}}\) Grade gilt die Entwickelung \[ f(y+xh) = f(y) + \sum_{i=1}^m \left\{f^{i-1}(y+h) - f^{i-1}(y)\right\} h^{i-1}\;\frac {\varphi_i (x)}{i!}\,, \] \[ \frac{f(y+xh) + f(y+h-xh)}{2} = \frac{f(y+h) + f(y)}{2}+\sum_{m=1}^{[n/2]} \left\{f^{2m-1} (y+h) - f^{2m-1}(y)\right\} h^{2m-1}\;\frac{\varphi_{2m} (x)}{(2m)!}, \] wo \([n/2]\) die grösste ganze Zahl bedeutet, welche in \(n/2\) enthalten ist. Diese Formeln werden verallgemeinert und dann auf \(f(x) = \log \varGamma (x)\) angewandt: \[ (1)\qquad\log \frac {\varGamma (y+x)}{\varGamma (y)} = x \log y+\sum_{i=2}^n (-1)^i\;\frac{\varphi_i(x)}{(i-1)i}\;y^{1-i} + R_n(x,y); \] für \(0<x\leqq \frac{1}{2}\) hat \(R_n\) die Form \((0<\vartheta<1):\) \[ (2)\qquad R_{2k-1} = \frac{\vartheta}{2k} \left\{\frac{\varphi_{2k}(x)}{2k-1}\;y- \frac{\varphi_{2k+1}(x)}{2k+1}\right\} y^{-2k}. \] Für \(x=\frac{1}{2}\) ist: \[ \log\;\frac {\varGamma(y+\frac {1}{2})}{\varGamma (y)} = \tfrac{1}{2} \log y +2 \sum_{m=1}^{k-1} (-1)^m (1-2^{-2m}) \frac {B_m}{(2m-1) 2m}\;y^{1-2m} +2 \vartheta (-1)^k (1-2^{-2k}) \frac{B_k}{(2k-1)2k}\;y^{1-2k}. \] Mit Hülfe der Formel (wo \(k=(n-1/2):)\) \[ \tfrac{h}{2}(e^h+1) = \left\{1+\frac{B_1}{2!}\;h^2 - \frac{B_2}{4!}\;h^4 + \cdots \pm \frac {B_k}{(2k)!}\;h^{2k} \right\} (e^h -1)+\frac{(-h)^{n+1}}{n!} \;\int_0^1 \varphi_n(x) e^{hx} dx \] werden die Resultate: \[ (-1)^m \int_0^1 \varphi_{2m} (x) dx = B_m = 2\;\frac {(2m)!}{(2\pi)^{2m}}\;\sum_{r=1}^\infty r^{-2m}, \] \[ \int_0^1 \varphi_{2m+1} (x)\text{cotg}\,\pi xdx = -2(-1)^m\;\frac {(2m+1)!}{(2\pi)^{2m+1}}\;\sum_{r=1}^\infty\;r^{-2m-1} \] \[ = -\frac {2m+1}{\pi}\;\int_0^1 \varphi_{2m} (x) \log \sin \pi xdx - (-1)^m B_m (2m+1)\;\frac {\log 2}{\pi} \] abgeleitet. Das Wesentliche dieser Arbeit hatte der Verf. in russischer Sprache bereits im Jahre 1888 in den Nachrichten der Universität Warschau veröffentlicht (F. d. M. 20, 424, 1888, JFM 20.0424.01). Daher gebührt Sonin, wie Hermite in seiner Antwort bemerkt, die Priorität für die auch von ihm in seiner Arbeit: Sur la fonction \(\log \varGamma (a)\) (J. für Math. 115, F. d. M. 26, 474, 1895, JFM 26.0474.01) entwickeln Formeln (1) und (2). Hermite definirt sodann die Function \(S_n (x) = \frac {\varphi_{n+1} (x)}{n+1}\) durch die Gleichung \(\frac {e^{xy} -1}{e^{y-1}} = \sum \frac {S_n (x)}{n!}\;y^n\) und leitet daraus \((c=0,1,2,\dots,n)\) die Entwickelungen ab: \[ \frac {S_n (x+z) -S_n (z)}{n!} = \sum\;\frac {S_c (x) z^{n-c}}{c! (n-c)!}, \quad \frac {x^n}{n!} = \sum\;\frac {S_c (x)}{n! (n-c)!} \] und, wenn \(f(y)\) eine beliebige Function bedeutet, \[ f(y+x) - f(y) = \sum\;\frac {f^c (y+1) -f^c (y)}{c!}\;S_c(x) - \frac{R_n} {n!}\,; \] für \(R_n\) werden zwei Darstellungen gegeben, deren eine ist: \[ R_n = \int_0^x (x-y)^n f^{n+1} (y+z) dz - \int_0^1 \left[S_n (x+1-z) -S_n (z) \right]f^{n+1} (y+z) dz. \] Weiter giebt er die Formel \((0<\theta<1):\) \[ \log \varGamma (a) = (a- \tfrac {1}{2}) \log a - a + \log \sqrt 2 \pi + \frac {\theta}{12a} \] und entwickelt eine beliebige Function \(f(y+x)\) nach Functionen \(\chi_n (x),\) \(\overline {\varphi}_n (x)\) (und ähnlichen), die durch die Gleichungen definirt sind: \[ \frac {e^{xy}}{e^y+1} = \sum\;\frac {\chi_n (x)}{n!}\;y^n, \quad \frac{y^s e^{xy}}{\varPhi(y)} = \sum\;\frac {\overline {\varphi}_n (x)} {n!}\;y^n, \quad (n = 0,1,2,\dots), \] wo \(\varPhi(x) = Ae^{ax} + Be^{bx} + \cdots + Le^{lx}\) ist. In seinem zweiten Briefe leitet Sonin für die Bernoulli'schen Zahlen verschiedene Sätze ab, z. B. \[ (-1)^n B_n = \frac {1}{2n+1} -\tfrac {1}{2} + {2n \choose 1} \frac {B_1}{2} - {2n \choose 3} \frac {B_2}{4} + \cdots + (-1)^n {2n \choose 2n-3} \frac {B_{n-1}}{2n-2}\,, \] \[ \frac{(-1)^nB_np(p^{2n}-1)}{n}=\sum^{p-1}_{k=0}(p-1-2k)(x+k)^{2n-1}+{2n-1\choose 1}\;\frac{B_1(p^2-1)}{1}\sum^{p-1}_{k=0}(x+k)^{2n-2}-\cdots \] \[ +(-1)^n{2n-1\choose 2}\frac{B_{n-1}(p^{2n-2}-1)}{n-1}\;\sum^{p-1}_{k=0}\;(x+k)^2, \] woraus sich ergiebt, dass die linke Seite der letzten Formel eine ganze Zahl ist. Zum Schlusse behandelt er das Integral \[ \Omega (z) = \int_0^\infty F (x)\;\frac {zdx}{x^2 + z^2}\,, \] wo \(F(x)\) zwischen den Grenzen nicht negativ wird, ferner \(\log\,\frac {\varGamma (a+\xi)}{\varGamma (a)}\) und berechnet die Werte von \(\Omega_k = (k + \frac {1}{2}) (k+1)\;\frac {B_k}{B_{k+1}}\) für \(k = 1, 2, 9, \infty.\)