Transformation du boustrophédon

La transformation du boustrophédon, ou algorithme boustrophédon, est une méthode mathématique permettant d'obtenir les coefficients du développement en série de Taylor des fonctions tangente, et sécante. Son nom fait référence au boustrophédon, une écriture dont le sens de lecture alterne d'une ligne à l'autre.

Construction du triangle boustrophédon

Schéma de principe des additions

Chaque ligne s'écrit dans le sens contraire de la précédente, en commençant par zéro ; et chaque terme se calcule en effectuant la somme du terme écrit précédemment et du terme écrit au-dessus, entre le précédent et lui. En partant de 1, on obtient le tableau triangulaire:

								1
						→	0		1
						1		1		0	←
				→	0		1		2		2
				5		5		4		2		0	←
		→	0		5		10		14		16		16
		61		61		56		46		32		16		0	←
→	0		61		122		178		224		256		272		272, etc.

Par exemple, le terme 14 est obtenu en faisant 10 + 4, ou 5 + 5 + 4 tandis que le terme 56 est obtenu en faisant 46 + 10 ou 16 + 16 + 14 + 10.

Plus formellement, le triangle des $B(n,k)$ est défini pour $0\leqslant k\leqslant n$ par :

$B(0,0)=1$

et pour $n\geqslant 1$ :

• Si $n$ est pair : $B(n,n)=0,B(n,k)=B(n,k+1)+B(n-1,k)$ pour $0\leqslant k\leqslant n-1$
• Si $n$ est impair : $B(n,0)=0,B(n,k)=B(n,k-1)+B(n-1,k-1)$ pour $1\leqslant k\leqslant n$.

Ce triangle est nommé triangle d'Euler-Bernoulli par Vladimir Arnold en 1992 ("parce que Pascal ne l'a pas considéré, et parce qu'Euler et Bernoulli ne l'ont pas considéré non plus") mais surtout car les nombres non nuls de gauche sont les nombres d'Euler et ceux de droite sont liés aux nombres de Bernoulli. L'appellation "boustrophédon" apparait sous la plume de Millar, Sloane, Young en 1996, reprise par Conway et Guy.

Il est répertorié comme suite A008280 de l'OEIS.

Développement de la fonction tangente

La transformation du boustrophédon permet d'obtenir le développement limité de la fonction tangente en 0.

La suite $(B(n,n))_{n\geqslant 1}$ des nombres formant le côté droit de ce triangle (sans le premier chiffre), soit 1, 0, 2, 0, 16, 0, 272, etc., donne la suite des coefficients du développement limité de la fonction tangente en 0 (en commençant par celui de $x^{1}$) :

$\tan(x)=1\times {\frac {x^{1}}{1!}}+0\times {\frac {x^{2}}{2!}}+2\times {\frac {x^{3}}{3!}}+0\times {\frac {x^{4}}{4!}}+16\times {\frac {x^{5}}{5!}}+0\times {\frac {x^{6}}{6!}}+272\times {\frac {x^{7}}{7!}}+o(x^{7})$

ce qui donne, après simplification :

$\tan(x)=x+{\frac {1}{3}}x^{3}+{\frac {2}{15}}x^{5}+{\frac {17}{315}}x^{7}+o(x^{7})$.

En poursuivant à l'infini, on obtient le développement en série de Taylor de la fonction tangente en 0 : $\operatorname {tan} x=\sum _{n=1}^{\infty }T_{n}{\frac {x^{2n-1}}{(2n-1)!}},|x|<\pi /2.$

La suite $(T_{n})_{n\geqslant 1}=(B(2n-1,2n-1))_{n\geqslant 1}=(\tan ^{(2n-1)}(0))=(1,2,16,272,...)$ est appelée suite des nombres tangents, ou parfois des nombres d'Euler de deuxième espèce.

Elle est répertoriée comme suite A000182 de l'OEIS, et avec les zéros intercalés, comme suite A350972 de l'OEIS.

Une définition par récurrence forte est $T_{1}=1,T_{n+1}=\sum _{k=1}^{n}{\binom {2n}{2k-1}}T_{k}T_{n+1-k}$ (application par la formule de Leibniz de $\tan '=1+\tan ^{2}$).

Le nombre $T_{n}$ est le nombre de permutations alternées ascendantes de longueur $2n-1$ ; il s'exprime en fonction des nombres de Bernoulli.

Calcul du développement de la fonction sécante

La suite $(B(n,0))_{n\geqslant 0}$ formant le côté gauche du triangle (avec le premier chiffre), soit 1, 0, 1, 0, 5, 0, 61, 0, etc., donne la suite des coefficients du développement limité de la fonction sécante en 0 (en commençant par celui de $x^{0}$, c'est-à-dire le terme constant).

$\sec(x)={\frac {1}{\cos x}}=1+0\times {\frac {x^{1}}{1!}}+1\times {\frac {x^{2}}{2!}}+0\times {\frac {x^{3}}{3!}}+5\times {\frac {x^{4}}{4!}}+0\times {\frac {x^{5}}{5!}}+61\times {\frac {x^{6}}{6!}}+0\times {\frac {x^{7}}{7!}}+o(x^{7})$.

En poursuivant à l'infini, on obtient le développement en série de Taylor de la fonction sécante : $\operatorname {sec} x=\sum _{n=0}^{\infty }E_{n}{\frac {x^{2n}}{(2n)!}},|x|<\pi /2.$

La suite $(E_{n})_{n\geqslant 0}=(B(2n,0))_{n\geqslant 0}=(\sec ^{(2n)}(0))=(1,1,5,61,1385,...)$ est appelée suite des nombres sécants, ou des nombres d'Euler. Elle est répertoriée comme suite A000364 de l'OEIS, et avec les zéros intercalés et alternance de signes, comme suite A122045 de l'OEIS.

Une définition par récurrence forte faisant intervenir les nombres tangents est $E_{0}=1,E_{n+1}=\sum _{k=0}^{n}{\binom {2n+1}{2k}}E_{k}T_{n+1-k}$ (application par la formule de Leibniz de $\sec '=\sec .\tan$).

Le nombre $E_{n}$ est le nombre de permutations alternées ascendantes de longueur $2n$ .

Autre présentation du triangle

Si l'on définit le triangle des nombres $E(n,k)$ pour $0\leqslant k\leqslant n$ par  :

$E(0,0)=1$

et pour $n\geqslant 1$ :

$E(n,0)=0$

$E(n,k)=E(n-1,n-1)+E(n-1,n-2)+\cdots +E(n-1,n-k)$ pour $1\leqslant k\leqslant n$, (ou $E(n,k)=E(n,k-1)+E(n-1,n-k)$),

autrement dit $E(n,k)$ est la somme des $k$ derniers termes de la ligne précédente,

on obtient le même triangle que précédemment sauf qu'une ligne sur deux a son sens inversé :

k n	0	1	2	3	4	5	6
0	1
1	0	1
2	0	1	1
3	0	1	2	2
4	0	2	4	5	5
5	0	5	10	14	16	16
6	0	16	32	46	56	61	61

Par exemple, $E(5,4)=5+5+4+2=14+2=16$.

Sous cette forme, le triangle est nommé triangle des nombres d'Entringer, ce dernier l'ayant étudié en 1966, et répertorié  A008282.

On a alors sur la diagonale alternativement les nombres sécants et les nombres tangents : $E(2n,2n)=E_{n}$, $E(2n-1,2n-1)=T_{n}$.

La suite $(E(n,n))_{n\geqslant 0}=(A_{n})=(1,1,1,2,5,16,\cdots )$ est la suite des nombres de permutations alternées ascendantes, répertoriée comme suite A000111 de l'OEIS.

Elle peut être définie par sa fonction génératrice exponentielle ${\frac {1}{\cos x}}+\tan x=\tan \left({\frac {x}{2}}+{\frac {\pi }{4}}\right)=\sum _{n=0}^{+\infty }A_{n}{\frac {x^{n}}{n!}},|x|<\pi /2$,

ou par récurrence forte par $A_{0}=1,2A_{n+1}=\sum _{k=0}^{n}{\binom {n}{k}}A_{k}A_{n-k}$.

Application à des valeurs approchées du nombre Pi

Le rayon de convergence de la série $\Sigma a_{n}x^{n}$ où $a_{n}={\frac {E(n,n)}{n!}}$ étant égal à $\pi /2$, $\lim {\frac {a_{n-1}}{a_{n}}}=\pi /2$, donc $\pi =\lim {\frac {2n\times E(n-1,n-1)}{E(n,n)}}$.

Par exemple, pour $n=9$, on obtient $\pi \approx 18\times 1385/7936\approx 3,1414$. L’intérêt de cette méthode est de donner des valeurs approchées de $\pi$ uniquement à partir d'additions d'entiers, d'une multiplication et d'une division.

Voir aussi

• Nombre d'Euler
• Permutation alternée

• Tableau triangulaire
• Triangle de Catalan, où $T(n,k)$ est la somme des $k+1$ premiers termes de la ligne précédente.