Numeri di Bernoulli

Jakob Bernoulli. *Summae Potestatum* in Ars Conjectandi, 1713

In matematica, i numeri di Bernoulli $B_{n}$ ^[1] costituiscono una successione di numeri razionali che gioca un ruolo importante in vari problemi. Accanto a essi conviene prendere in considerazione i polinomi di Bernoulli che si possono considerare una loro generalizzazione.

Genesi storica

Questi numeri furono individuati quasi contemporaneamente ma indipendentemente da Kōwa Seki nel 1712 e da Jakob Bernoulli nel 1713^[2]. Bernoulli li tratta nella sua opera Ars Conjectandi, in relazione con le forme chiuse per le somme di potenze di interi successivi

\sum_{k = 1}^{n} k^{m} = 1^{m} + 2^{m} + \dots + n^{m}

per valori interi positivi fissati di $m .$

Queste forme chiuse erano state individuate già nel 1631 da Johann Faulhaber^[3] cui Bernoulli fa riferimento. Dopo la morte di Bernoulli, nel 1721 Abraham de Moivre ed Eulero diedero ai numeri $B_{n}$ il nome con il quale sono tuttora conosciuti^[2].

Le precedenti somme sono esprimibili per ogni $m$ come polinomi in $n$ di grado $m + 1 .$ La formula rivelata e forse scoperta^[4] ma non dimostrata da Jakob Bernoulli, scritta in notazione moderna utilizzando la notazione del fattoriale decrescente, è:

\sum_{k = 1}^{n} k^{m} = \frac{1}{m + 1} n^{m + 1} + \frac{1}{2} n^{m} + \sum_{k = 2}^{n} \frac{m^{\underline{k - 1}}}{k!} B_{k} n^{m - k + 1} .

Bernoulli però nel suo trattato non considerava quelli che oggi per noi sono invece i primi due numeri della sequenza numerica che porta il suo nome per cui oggi preferiamo esprimere quella stessa formula nel seguente modo

\sum_{k = 1}^{n} k^{m} = \frac{1}{m + 1} \sum_{k = 0}^{m} (\binom{m + 1}{k}) B_{k} n^{m - k + 1},

dove si è utilizzata la variante $B_{1} = 1 / 2$ dei numeri di Bernoulli.

Esempi

(B_{0} = 1, B_{1} = 1 / 2, B_{2} = 1 / 6, B_{3} = 0, B_{4} = - 1 / 30) .

Nel caso $m = 1$ abbiamo

1 + 2 + \dots + n = \frac{1}{2} (B_{0} n^{2} + 2 B_{1} n^{1}) = \frac{n (n + 1)}{2} = \frac{1}{2} n^{2} + \frac{1}{2} n .

Nel caso $m = 4$ abbiamo

1^{4} + 2^{4} + 3^{4} + \dots + n^{4} = \frac{1}{5} \sum_{j = 0}^{4} (\binom{5}{j}) B_{j} n^{5 - j} =

= \frac{1}{5} (B_{0} n^{5} + 5 B_{1} n^{4} + 10 B_{2} n^{3} + 10 B_{3} n^{2} + 5 B_{4} n) =

= \frac{1}{5} n^{5} + \frac{1}{2} n^{4} + \frac{1}{3} n^{3} - \frac{1}{30} n .

Definizione ricorsiva

I numeri di Bernoulli, nella variante $B_{1} = - 1 / 2$ , possono essere calcolati usando la seguente formula di ricorrenza:

\sum_{j = 0}^{m} (\binom{m + 1}{j}) B_{j} = 0, con la condizione iniziale B_{0} = 1,

che equivale a:

\sum_{j = 0}^{m - 1} (\binom{m + 1}{j}) B_{j} + (m + 1) B_{m} = 0,

da cui la forma esplicita:

B_{m} = - \frac{1}{m + 1} \sum_{j = 0}^{m - 1} (\binom{m + 1}{j}) B_{j} .

L'algoritmo di Ada Lovelace

Nella nota G delle Note di Ada Lovelace sull'analytical engine del 1842^[5] è stato descritto per la prima volta un algoritmo per la costruzione dei numeri di Bernoulli con una macchina in grado di eseguire calcoli automatici. L'algoritmo di Ada Lovelace per i numeri di Bernoulli si basa sulla formula ricorsiva che abbiamo visto anche se per riconoscerla si deve tener conto che anche ai suoi tempi, come in quelli di Bernoulli, non erano considerati i primi due numeri della sequenza^[6].

Tavola dei numeri di Bernoulli

Si può dimostrare che $B_{n} = 0$ per tutti gli $n$ dispari maggiori di 1.

I primi valori diversi da 0 sono i seguenti:

n	0	1	2	4	6	8	10	12	14	16	18	20
B_n	1	±1/2^[7]	1/6	-1/30	1/42	-1/30	5/66	-691/2730	7/6	-3617/510	43867/798	-174611/330

I numeri di Bernoulli compaiono anche negli sviluppi in serie di Taylor della tangente e della tangente iperbolica, nella formula di Eulero-Maclaurin e nelle espressioni di certi valori della funzione zeta di Riemann.

Funzioni generatrici

I numeri di Bernoulli possono anche essere definiti usando delle funzioni generatrici esponenziali sviluppate in serie di Maclaurin. Per la variante della sequenza bernoulliana con $B_{1} = - 1 / 2$ abbiamo:

\frac{x}{e^{x} - 1} = : \sum_{n = 0}^{\infty} B_{n} \frac{x^{n}}{n!},

mentre sommando $x$ ai due membri della precedente^[8], troviamo la funzione generatrice della variante con $B_{1} = 1 / 2$

\frac{- x}{e^{- x} - 1} = : \sum_{n = 0}^{\infty} B_{n} \frac{x^{n}}{n!} .

Queste possono considerarsi uguaglianze fra serie formali di potenze; in questo caso per la convergenza della serie si chiede che $x$ abbia valore assoluto minore di $2 π$ (il raggio di convergenza della serie stessa).