Extracting Bernoulli polynomials from their generating function

The generating function for Bernoulli polynomials is

$$\frac{te^{tx}}{e^t-1} = \sum_{n=0}^\infty B_n(x) \frac{t^n}{n!}$$

The only way that I know of to get the coefficients out of this is to use Taylor's theorem. But to use Taylor's theorem, it needs to have derivatives at zero, and this doesn't, because there is always a $e^t-1$ in some denominator, which at zero is undefined.

However, if instead of evaluating things at zero, I take the limit as $t \to 0$, then it looks like it works. As an example, for the second "term",

$$\lim_{t\to 0} \left(\frac{te^{tx}}{e^t-1} \right)^{(1)}\frac{t^0}{n!} = \frac{x-1/2}{n!}$$

Which works, since $B_1(x) = x - 1/2$.

My questions are,

1) Does this actually work, and if so, why? Are there other cases where we have formal power series which sum to something that doesn't have a Taylor series at zero?

2) Is there a better way to extract Bernoulli polynomials from that generating function (not using some other formula)?

• Here is a practical way to do it: First derive the relation $\frac{dB_n(x)}{dx} = n B_{n-1}(x)$ (from the generating series) and use this to generate the polynomials recursively as $B_n(x) = B_n(0) + \int_0^x nB_{n-1}(x)dx$. You need to already know $B_n(0)$ (which are just the bernouilli numbers) for this method though. Commented Apr 17, 2015 at 21:28

As we are not permitted to use anything else we start from the definition $$\frac{te^{tx}}{e^t-1} = \sum_{n\ge 0} B_n(x) \frac{t^n}{n!}$$ and obtain $$B_n(x) = \frac{n!}{2\pi i} \int_{|t|=\epsilon} \frac{1}{t^{n+1}} \frac{te^{tx}}{e^t-1} \; dt$$ and hence $$[x^q] B_n(x) = \frac{n!}{q!\times 2\pi i} \int_{|t|=\epsilon} \frac{1}{t^{n+1}} \frac{t^{q+1}}{e^t-1} \; dt \\ = \frac{n!}{q!\times 2\pi i} \int_{|t|=\epsilon} \frac{1}{t^{n-q+1}} \frac{t}{e^t-1} \; dt.$$ We see that this is zero for $q>n$ and moreover, $$[x^q] B_n(x) = \frac{n!}{q!} \frac{B_{n-q}}{(n-q)!} = {n\choose q} B_{n-q}.$$

This means we have reduced the problem of computing Bernoulli polynomials to the problem of computing Bernoulli numbers, since as we now know, $$B_n(x) = \sum_{q=0}^n {n\choose q} B_{n-q} x^q.$$

This can be done by differentiating the generating function $$\frac{t}{e^t-1} = \sum_{n\ge 0} B_n \frac{t^n}{n!}$$ to get $$\frac{1}{e^t-1} - \frac{t}{(e^t-1)^2} e^t = \sum_{n\ge 1} B_n \frac{t^{n-1}}{(n-1)!}.$$

The left is $$\frac{1}{e^t-1} - \frac{1}{e^t-1} \frac{te^t}{e^t-1} \\ = \frac{1}{t} + \left(-\frac{1}{t} + \frac{1}{e^t-1}\right) - \frac{1}{t} \frac{te^t}{e^t-1} - \left(-\frac{1}{t} + \frac{1}{e^t-1}\right) \frac{te^t}{e^t-1}.$$

Extracting coefficients from these exponential generating functions we obtain for $n\ge 0$

$$B_{n+1} = \frac{B_{n+1}}{n+1} - \frac{B_{n+1}(1)}{n+1} - \sum_{q=0}^n {n\choose q} \frac{B_{q+1}}{q+1} B_{n-q}(1).$$

Observe that $$\frac{te^t}{e^t-1} = t + \frac{t}{e^t-1}$$

and hence $$B_{n+1}(1) = B_{n+1} \quad\text{except for}\quad B_1(1) = 1 + B_1.$$

This gives for $n\ge 1$ $$B_{n+1} = - {n\choose n-1} \frac{B_n}{n} - \sum_{q=0}^n {n\choose q} \frac{B_{q+1}}{q+1} B_{n-q} \\ = - B_n - \frac{B_{n+1}}{n+1} - \sum_{q=0}^{n-1} {n\choose q} \frac{B_{q+1}}{q+1} B_{n-q}$$ or $$B_{n+1}\frac{n+2}{n+1} = - B_n - \sum_{q=0}^{n-1} {n\choose q} \frac{B_{q+1}}{q+1} B_{n-q}$$

and finally $$B_{n+1} = - \frac{n+1}{n+2} B_n - \frac{n+1}{n+2} \sum_{q=0}^{n-1} {n\choose q} \frac{B_{q+1}}{q+1} B_{n-q}.$$

This recurrence lets us calculate the $B_n$ starting from $B_0=1$ and $B_1=-\frac{1}{2}.$

Addendum Wed Apr 22 22:47:07 CEST 2015. Since I did not have a reference at hand when I wrote the above I failed to see that there is a much simpler identity which results from writing $$t = (e^t-1) \sum_{n\ge 0} B_n \frac{t^n}{n!}.$$

This is $$t = e^t \sum_{n\ge 0} B_n \frac{t^n}{n!} - \sum_{n\ge 0} B_n \frac{t^n}{n!}.$$

Extracting coefficients for $n\ge 2$ and observing the convolution of two exponential generating functions in the first term we have $$0 = \sum_{q=0}^n {n\choose q} B_q - B_n = \sum_{q=0}^{n-1} {n\choose q} B_q.$$ This yields $$B_{n-1} {n\choose n-1} = -\sum_{q=0}^{n-2} {n\choose q} B_q$$ or $$B_{n-1} = -\frac{1}{n} \sum_{q=0}^{n-2} {n\choose q} B_q \quad\text{or}\quad B_m = -\frac{1}{m+1} \sum_{q=0}^{m-1} {m+1\choose q} B_q.$$

• The first step is a bit of a leap for me. What is $\int _{|t|=\epsilon}$? I haven't come across notation like that before. Commented Apr 18, 2015 at 1:42
• @dannynoest: It is a contour integration from complex variables. Commented Apr 18, 2015 at 1:48
• Start additional readings at Wikipedia I and Wikipedia II. Commented Apr 18, 2015 at 1:49