Cuvée 3.14

Vin du nouvel an

Bonne année en retard ! La cuvée 3.14 de Foillard est un petit délice que j’ai eu la chance de boire entre amis le soir du nouvel an (avec en prime un calembour mathématique).

\(\text{Côte du Py }\approx\text{ Côte du }\pi\)

Peut-être le plus grand vin du Beaujolais ? La remarquable cuvée 3.14 de Foillard est issue de vignes de plus de 100 ans de la colline volcanique de la Côte du Py et uniquement dans les meilleurs millésimes. Ce vin éthéré et élégant, produit en quantités minuscules, est exceptionnellement rare. Dans une leçon magistrale d’équilibre, Foillard parvient à faire ressortir des arômes incroyablement complexes et toute la profondeur du fruit, tout en conférant à son vin une texture légère, séduisante et soyeuse.

Traduit et adapté de : https://cellanextdoor.com

\(\pi\)-ètre formulation

Plus tôt cette année, un peu avant la journée de \(\pi\), j’ai vu dans une infolettre mathématique hebdomadaire cette rubrique :

Je dois avouer que l’infolettre est d’ordinaire assez bien ficelée et que je trouve les contenus souvent intéressants et pertinents, surtout pour un enseignant au secondaire. Sans surprise j’étais donc enthousiasmé par la suite présentée, une suite simple, nouvelle en ce qui me concerne, et qui tend vers une constante célébrée. Je l’ai donc soumise avec entrain à certains élèves qui étaient encore dans la classe au moment où j’ai consulté l’infolettre. J’ai même osé dire à un élève dont je savais les parents professeurs de mathématiques : « Tiens, tu montreras ça à ton père ! Je ne sais pas encore pourquoi ça tend vers \(\pi\), mais c’est intriguant ! »

Curieux de connaître la rapidité avec laquelle la suite converge vers \(\pi\), en arrivant à la maison j’ai concocté un petit script en python :

Aïe aïe aïe ! J’ai été surpris de constater qu’en réalité la suite ne tend pas vers \(\pi\), malgré ce qu’insinuait l’infolettre. Elle semble tendre vers un autre nombre pas loin, à près de \(0,\!00005\) de \(\pi\). Close, but no cigar. Oups ! Un petit courriel plus tard pour expliquer mon trop plein d’enthousiasme (et la bévue) au père de mon élève, et je me mettais au travail. « … numbers that converge to a recognisable number ». Ah ? Lequel ? Je ne vois pas.

Feu l’Inverseur de Plouffe

Étant incapable d’utiliser à bon escient la référence du Journal of Recreational Mathematics, mon premier réflexe a été de tenter de me rendre sur le site de l’Inverseur de Plouffe, pour voir si ce nombre mystère s’y trouvait. Je dis « tenter » car c’est à ce moment que j’ai découvert que l’Inverseur de Plouffe n’est plus en ligne… depuis… 12 ans. Re-oups. Ça ne nous rajeunit pas, disons.

La fraction médiane ? (non)

Je dois avouer que j’ai ensuite passé un peu (trop) de temps à investiguer du côté des fractions médianes. La fraction médiane de \(\frac{a}{b}\) et \(\frac{c}{d}\) est \(\frac{a+c}{b+d}\). Et dans la suite, on somme les dividendes (numérateurs) et les diviseurs (dénominateurs) des deux termes précédents. On note parfois \(\oplus\) l’opération de faire la somme des numérateurs et des dénominateurs \[\frac{a}{b}\oplus \frac{c}{d} = \frac{a+c}{b+d}\]que les anglophones nomment « the freshman’s sum » parce qu’il s’agit d’une façon erronée d’effectuer l’addition (traditionnelle) de fractions, erreur que répètent parfois certains élèves novices. L’intérêt de la fraction médiane est que si \[\frac{a}{b}<\frac{c}{d}\] et que \(a\), \(b\), \(c\) et \(d\) sont positifs, alors \[\frac{a}{b}<\frac{a+c}{b+d}<\frac{c}{d}\]

Les pentes des côtés du parallélogramme sont \(\frac{a}{b}\) et \(\frac{c}{d}\) et \(\frac{a}{b}<\frac{c}{d}\).
La pente de la diagonale du parallélogramme est \(\frac{a+c}{b+d}\) et \(\frac{a}{b}<\frac{a+c}{b+d}<\frac{c}{d}\).

On peut s’en convaincre algébriquement en calculant les différences strictement positives entre \(\frac{a+c}{b+d}\) et \(\frac{a}{b}\) \begin{align*}\frac{a+c}{b+d}-\frac{a}{b} &=\frac{b(a+c)}{b(b+d)}-\frac{a(b+d)}{b(b+d)} \\ \\ &=\frac{ab+bc}{b^{2}+bd}-\frac{ab+ad}{b^{2}+bd} \\ \\ &= \frac{ab + bc-ab-ad}{b^{2}+bd} \\ \\ &=\frac{bc-ad}{b^{2}+bd}\\ \\ &=\frac{\frac{bc}{bd}-\frac{ad}{bd}}{\frac{b^{2}}{bd}+\frac{bd}{bd}} \\ \\ &= \frac{\frac{c}{d}-\frac{a}{b}}{\frac{b}{d}+1} > 0 \end{align*}car \begin{align*}\frac{c}{d}&>\frac{a}{b} \\ \\ \frac{c}{d}-\frac{a}{b}&>0\end{align*}et \[\frac{b}{d}+1>0\]et entre \(\frac{c}{d}\) et \(\frac{a+c}{b+d}\) \begin{align*}\frac{c}{d}-\frac{a+c}{b+d}&=\frac{c(b+d)}{d(b+d)}-\frac{d(a+c)}{d(b+d)}\\ \\ &=\frac{bc+cd}{bd+d^{2}}-\frac{ad+cd}{bd+d^{2}}\\ \\ &=\frac{bc+cd-ad-cd}{bd+d^{2}} \\ \\ &=\frac{bc-ad}{bd+d^{2}} \\ \\ &=\frac{\frac{bc}{bd}-\frac{ad}{bd}}{\frac{bd}{bd}+\frac{d^{2}}{bd}} \\ \\ &=\frac{\frac{c}{d}-\frac{a}{b}}{1+\frac{d}{b}} > 0\end{align*} car \[\frac{c}{d}-\frac{a}{b}>0\] et \[1+\frac{b}{d}>0\]Les fractions médianes peuvent donc nous servir à trouver des approximations rationnelles à des nombres irrationnels sans trop de peine. Par exemple, pour \(\pi\), on peut utiliser les fractions \(\frac{3}{1}\) et \(\frac{4}{1}\) puisque \[\frac{3}{1}<\pi< \frac{4}{1}\]La première fraction médiane est \[\frac{3}{1}\oplus \frac{4}{1}=\frac{3+4}{1+1} = \frac{7}{2}\]Puisque \[\frac{3}{1}< \pi <\frac{7}{2}\]la deuxième fraction médiane serait \[\frac{3}{1}\oplus\frac{7}{2} = \frac{3+7}{1+3} = \frac{10}{3}\]Puisque \[\frac{3}{1}<\pi <\frac{10}{3}\]la suivante serait \[\frac{3}{1}\oplus \frac{10}{3} = \frac{10+3}{1+3} = \frac{13}{4}\]Puis comme \[\frac{3}{1}<\pi <\frac{13}{4}\]la suivante serait \[\frac{3}{1}\oplus \frac{13}{4} = \frac{3+13}{1+4} = \frac{16}{5}\]Puisque \(\pi\) était beaucoup plus près de la fraction de départ \(\frac{3}{1}\) que de la fraction de départ \(\frac{4}{1}\), on semble s’approcher de \(\pi\) par la droite, en d’autres mots, la fraction médiane est pour l’instant toujours supérieure à \(\pi\). Bien sûr, cela ne sera pas toujours le cas. Poursuivons encore quelques étapes,\[\frac{3}{1}<\pi <\frac{16}{5}\]nous donne \[\frac{3}{1}\oplus \frac{16}{5} = \frac{3+16}{1+5} =\frac{19}{6}\]puis \[\frac{3}{1}<\pi <\frac{19}{6}\]nous donne \[\frac{3}{1}\oplus \frac{19}{6} = \frac{22}{7}\]approximation célébrée de \(\pi\). \[\frac{3}{1}<\pi <\frac{22}{7}\]nous donne une enième approximation par la droite \[\frac{3}{1}\oplus \frac{22}{7} = \frac{3+22}{1+7} = \frac{25}{8}\]avant que \[\frac{22}{7}<\pi <\frac{25}{8}\]nous donne une première approximation par la gauche \[\frac{22}{7}\oplus\frac{25}{8} = \frac{25+22}{7+8} = \frac{47}{15}\]Les approximations obtenues avec les fractions médianes sont sensibles au choix des fractions de départ et convergent lentement par rapport à d’autres méthodes. Les meilleures approximations rationnelles sont celles obtenues à l’aide des fractions continues, et en comparaison, celles des fractions médianes font pâle figure. Il nous a fallu de nombreuses étapes pour trouver \[\frac{47}{15}\approx 3,\!1333\]une approximation qui n’est pas, en soi, exceptionnelle.

En outre, dans la suite,

en considérant les divisions comme des fractions, on approche le nombre mystère en trouvant la fraction médiane des deux dernier termes. Or, sans plan clair et agacé par le fait qu’on utilise uniquement les deux dernières fractions et jamais celle qui précède les deux dernières, comme c’est arrivé dans mon exemple pour \(\pi\), j’ai eu l’impression d’emprunter un cul-de-sac. Or, ce n’était pas en vain : à quelle notion pourrait-on faire référence lorsqu’on considère les sommes des deux derniers termes ?

La suite de Fibonacci (oui !)

En effet, si on considère la suite de Fibonacci, \[1, \ 1, \ 2, \ 3, \ 5, \ 8, \ 13, \ 21, \ 34, \ 55, \ 89, \ 144, \ 233, \ 377, \ 610\, \ \dots\]on peut observer que les diviseurs correspondent aux quintuple des termes de de la suite.\[\textcolor{Red}{5}, \ \textcolor{Red}{5}, \ 10, \ 15, \ 25, \ 40, \ 65, \ 105, \ 170, \ 275, \ 445, \ 720, \ , 1\,165,\ 1\,885, \ 3\,050\, \ \dots\]Pour les dividendes, c’est moins évident, mais ce sont les sextuple des termes de la suite, avec un déphasage \[6, \ 6, \ \textcolor{Blue}{12}, \ \textcolor{Blue}{18}, \ 30, \ 48, \ 78, \ 126, \ 204, \ 330, \ 534, \ 864, \ 1\,398, \ 2\,262, \ 3\,660\, \ \dots\]Si \(F_{n}\) représente le \(n^{\text{e}}\) terme de la suite de Fibonacci, alors les termes de la suite de l’infolettre sont \[\displaystyle \frac{6F_{n+2}}{5F_{n}}\]En se rappelant que \(F_{n+2} = F_{n+1}+F_{n}\), o peut évaluer la limite \begin{align*}\lim_{n\to \infty}\frac{6F_{n+2}}{5F_{n}}&=\frac{6}{5}\cdot \lim_{n\to \infty}\frac{F_{n+2}}{F_{n}} \\ \\ &=\frac{6}{5}\cdot \lim_{n\to \infty}\frac{F_{n+1}+F_{n}}{F_{n}} \\ \\ &= \frac{6}{5}\left( \lim_{n\to \infty} 1 + \frac{F_{n+1}}{F_{n}}\right) \\ \\ &=\frac{6}{5}\left(1+\lim_{n \to \infty}\frac{F_{n+1}}{F_{n}}\right)\\ \\ &= \frac{6}{5}\left(1 + \varphi\right)\end{align*}car tel que vu précédemment sur ce blogue, le rapport des termes consécutifs de la suite tendent vers \(\varphi = \frac{1+\sqrt{5}}{2}\).  Sachant en plus que \begin{align*}1 + \varphi &= 1+ \frac{1+\sqrt{5}}{2} \\ \\ &=\frac{3+\sqrt{5}}{2} \\ \\ &= \frac{6+2\sqrt{5}}{4} \\ \\ &= \frac{1+2\sqrt{5}+5}{4} \\ \\ &=\frac{(1+\sqrt{5})^{2}}{2^{2}} \\ \\ &=\left(\frac{1+\sqrt{5}}{2}\right)^{2} \\ \\ &=\varphi^{2}\end{align*}il est possible d’écrire \[\frac{6}{5}\left(1+\varphi\right) = \frac{6}{5}\varphi^{2}\]La suite tend donc vers \[\frac{6}{5}\varphi^{2} = \frac{9+3\sqrt{5}}{5} \approx 3,\!141640786499873817845504201238\]un nombre qui est à moins de \(0,\!000048132910080579383\) de \(\pi\). Grâce à cette valeur, je suis ensuite tombé sur cette page chez Futility Closet, contenant la même référence au Journal of Recreational Mathematics, mais hélas sans plus de détails.

 

PS. Au moment où j’allais publier le billet, j’ai découvert une version de l’Inverseur disponible ici : http://wayback.cecm.sfu.ca/projects/ISC/ISCmain.html

J’avais foi en l’Inverseur et j’aurai eu raison ! En effet, on y trouve :
\[18 + 6\sqrt{5} \approx 31,\!41640786499873817\]ce qui nous aurait bien sûr mis sur la bonne piste.

 

Quadrilat… aire ?

Dans un triangle, la formule de Héron nous permet de calculer l’aire du triangle avec les mesures de ses côtés. Si \(K\) est l’aire du triangle dont les côtés mesurent \(a\), \(b\) et \(c\), alors \[K = \sqrt{s\, (s-a)(s-b)(s-c)}\]dans laquelle \(s\) est le demi-périmètre \[s = \frac{a+b+c}{2}\]Nous avons rencontré la formule de Héron ici et ici. Existe-t-il une telle formule pour quadrilatère ? En d’autres mots, connaissant les mesures des quatre côtés d’un quadrilatère, est-il possible de déterminer son aire et de l’exprimer à l’aide d’une formule ?

Hélas, non, ce n’est pas possible. À partir d’un quadrilatère ayant pour côtés \(a\), \(b\), \(c\) et \(d\) comme côtés, il est possible de former une infinité de quadrilatères avec les mêmes mesures de côtés. Il suffit d’imaginer des charnières aux sommets du quadrilatère. On obtient un quadrilatère articulé.

Bon. Si ce petit billet se terminait ici, cela serait un peu ordinaire. Est-il possible de trouver quelque chose d’intéressant quand même concernant le quadrilatère, ses côtés et son aire ? Bien sûr ! Il faut cependant fixer, par exemple, les diagonales ou certains angles. Supposons que l’angle entre \(a\) et \(b\) soit \(\theta\) et que l’angle entre \(c\) et \(d\) soit \(\phi\).

Les « charnières » ne sont plus libres de bouger. Le quadrilatère est fixe. Plus précisément, si on trace la diagonale \(x\) comme dans la figure ci-dessous, on obtient deux triangles.

Il n’est pas très difficile ni très spectaculaire d’obtenir une expression pour l’aire du quadrilatère. Il suffit de faire la somme des aires des deux triangles. La formule d’aire trigonométrique du triangle nous donne l’expression recherchée instantanément.

\[K = \frac{ab\sin(\theta)}{2} + \frac{cd\sin(\phi)}{2}\]

Cela reste tout de même un peu décevant et c’est pourquoi, chers lecteurs, je vous invite à suivre ces quelques étapes algébriques et trigonométriques supplémentaires.

La formule de Bretschneider

\[K = \frac{ab\sin(\theta)}{2} + \frac{cd\sin(\phi)}{2}\]

On double,

\[2K = ab\sin(\theta) + cd\sin(\phi)\]

et on élève au carré, \[4K^{2} = a^{2}b^{2}\sin^{2}(\theta) + 2abcd\sin(\theta)\sin(\phi) + c^{2}d^{2}\sin^{2}(\phi)\]

et on quadruple

\[16K^{2}= 4\left( a^{2}b^{2}\sin^{2}(\theta) + 2abcd\sin(\theta)\sin(\phi) + c^{2}d^{2}\sin^{2}(\phi)\right)\]

D’autre part, en utilisant la loi des cosinus deux fois, une fois dans chaque triangle,

\[x^{2} = a^{2} + b^{2}-2ab\cos(\theta)\] et \[x^{2} = c^{2}+d^{2}-2cd\cos(\phi)\]

on obtient \[a^{2} + b^{2}-2ab\cos(\theta) = c^{2}+d^{2}-2cd\cos(\phi)\]

On peut réarranger les termes pour obtenir \begin{align*}a^{2}+b^{2}-c^{2}-d^{2}&=2ab\cos(\theta)-2cd\cos(\phi) \\ \\ a^{2}+b^{2}-c^{2}-d^{2}&=2\left(ab\cos(\theta)-cd\cos(\phi)\right)\end{align*}

En élevant au carré et en regroupant les termes du même côté de l’égalité, on obtient

\begin{align*}\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}&=\left(2\left(ab\cos(\theta)-cd\cos(\phi)\right)\right)^{2} \\ \\ \left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}&=4\left(ab\cos(\theta)-cd\cos(\phi)\right)^{2} \\ \\ 0&=4\left(ab\cos(\theta)-cd\cos(\phi)\right)^{2}-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &= 4\left(a^{2}b^{2}\cos^{2}(\theta)-2abcd\cos(\theta)\cos(\phi)+c^{2}d^{2}\cos^{2}(\phi)\right)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}\end{align*}

Il est ensuite possible d’additionner ces deux équations. On obtient

\begin{align*}16K^{2}&=4\left( a^{2}b^{2}\sin^{2}(\theta) + 2abcd\sin(\theta)\sin(\phi) + c^{2}d^{2}\sin^{2}(\phi)\right) \\ &\phantom{=}+4\left(a^{2}b^{2}\cos^{2}(\theta)-2abcd\cos(\theta)\cos(\phi)+c^{2}d^{2}\cos^{2}(\phi)\right)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &=4\big(a^{2}b^{2}\sin^{2}(\theta) + 2abcd\sin(\theta)\sin(\phi) + c^{2}d^{2}\sin^{2}(\phi)+a^{2}b^{2}\cos^{2}(\theta) \\ &\phantom{=}-2abcd\cos(\theta)\cos(\phi)+c^{2}d^{2}\cos^{2}(\phi)\big)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &= 4\big(a^{2}b^{2}(\sin^{2}(\theta) +\cos^{2}(\theta))+c^{2}d^{2}(\sin^{2}(\phi) + \cos^{2}(\phi))\\ &\phantom{=}-2abcd(\cos(\theta)\cos(\phi)-\sin(\theta)\sin(\phi)\big)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &=4\left(a^{2}b^{2}+c^{2}d^{2}-2abcd\cos(\theta+\phi)\right)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &=4\left(a^{2}b^{2}+c^{2}d^{2}\right)-8abcd\cos(\theta+\phi)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}\end{align*}en utilisant, à l’avant-dernière étape \[\sin^{2}(x) + \cos^{2}(x) = 1\]et \[\cos(x + y) = \cos(x)\cos(y)-\sin(x)\sin(y)\]On note au passage que le choix de \(\theta\) et \(\phi\) comme étant les angles entre \(a\) et \(b\) et entre \(c\) et \(d\) plutôt que les angles entre \(a\) et \(d\) et \(b\) et \(c\) n’a pas d’importance. Puisque la somme des mesures des angles intérieurs d’un quadrilatère est \(360^{\circ}\), si les angles entre les côtés \(a\) et \(d\) et \(b\) et \(c\) sont respectivement \(\alpha\) et \(\beta\), on a \[(\theta + \phi) + (\alpha + \beta) = 360^{\circ}\] ce qui implique \[\cos(\theta + \phi)  = \cos(\alpha + \beta)\]Cela étant dit, à \[16K^{2}=4\left(a^{2}b^{2}+c^{2}d^{2}\right)-8abcd\cos(\theta+\phi)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}\]on ajoute et retranche ensuite \(8abcd\) à droite et on factorise le trinôme carré parfait. On obtient \begin{align*}16K^{2} &= 4\left(a^{2}b^{2}+c^{2}d^{2}\right)+8abcd-8abcd-8abcd\cos(\theta+\phi)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &=4\left(a^{2}b^{2}+2abcd+c^{2}d^{2}\right)-8abcd\left(1+\cos(\theta+\phi)\right)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &=4\left(ab+cd\right)^{2}-8abcd\left(1+\cos(\theta+\phi)\right)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}\end{align*}Sachant que \begin{align*}1+\cos(x)&=1+\cos\left(\frac{x}{2}+\frac{x}{2}\right) \\ \\ &=1 + \cos\left(\frac{x}{2}\right)\cos\left(\frac{x}{2}\right)-\sin\left(\frac{x}{2}\right)\sin\left(\frac{x}{2}\right) \\ \\ &=1+\cos^{2}\left(\frac{x}{2}\right)-\sin^{2}\left(\frac{x}{2}\right) \\ \\ &=1+\cos^{2}\left(\frac{x}{2}\right)-\left(1-\cos^{2}\left(\frac{x}{2}\right)\right) \\ \\ &=2\cos^{2}\left(\frac{x}{2}\right)\end{align*}on peut poser \(x = \theta + \phi\) et l’équation devient \begin{align*}16K^{2}&= 4\left(ab+cd\right)^{2}-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2} \\ \\ &= 4\left(ab+cd\right)^{2}-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)\end{align*}Suivent quelques étapes algébriques élémentaires dans lesquelles on factorise des différences de carrés et des trinômes carrés parfaits.\begin{align*}16K^{2}&=4\left(ab+cd\right)^{2}-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &=\left(2ab+2cd\right)^{2}-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)^{2}-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)\\ \\ &=\left(2ab+2cd-\left(a^{2}+b^{2}-c^{2}-d^{2}\right)\right)\left(2ab+2cd+\left(a^{2}+b^{2}-c^{2}-d^{2}\right)\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &=\left(c^{2}+2cd+d^{2}-\left(a^{2}-2ab+b^{2}\right)\right)\left(a^{2}+2ab+b^{2}-\left(c^{2}-2cd+d^{2}\right)\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &=\left(\left(c+d\right)^{2}-\left(a-b\right)^{2}\right)\left(\left(a+b\right)^{2}-\left(c-d\right)^{2}\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &=\left(c+d-(a-b)\right)\left(c+d+(a-b)\right)\left(a+b-(c-d)\right)\left(a+b+(c-d)\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)\\ \\ &=\left(b+c+d-a\right)\left(a+c+d-b\right)\left(a+b+d-c\right)\left(a+b+c-d\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)\end{align*}

Il ne reste que quelques étapes simples. On introduit le demi périmètre \(s\) : \[s = \frac{a+b+c+d}{2}\]On obtient enfin \begin{align*}16K^{2}&=\left(b+c+d-a\right)\left(a+c+d-b\right)\left(a+b+d-c\right)\left(a+b+c-d\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &= (2)\left(\frac{b+c+d-a}{2}\right)(2)\left(\frac{a+c+d-b}{2}\right)(2)\left(\frac{a+b+d-c}{2}\right)(2)\left(\frac{a+b+c-d}{2}\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &=16\left(\frac{a+b+c+d}{2}-\frac{2a}{2}\right)\left(\frac{a+b+c+d}{2}-\frac{2b}{2}\right)\left(\frac{a+b+c+d}{2}-\frac{2c}{2}\right)\left(\frac{a+b+c+d}{2}-\frac{2d}{2}\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \\ \\ &=16\left(s-a\right)\left(s-b\right)\left(s-c\right)\left(s-d\right)-16abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right) \end{align*}En divisant par \(16\) et en retenant la racine positive, on obtient \[K = \sqrt{\left(s-a\right)\left(s-b\right)\left(s-c\right)\left(s-d\right)-abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)}\]La formule de Bretschneider. Voilà !

Formule de Brahmagupta

En bonus ! La formule de Bretschneider se simplifie dans le cas d’un quadrilatère inscriptible (voir la toute fin de ce billet). En effet, dans un tel quadrilatère, les angles opposés sont supplémentaires. On a \(\theta + \phi = 180^{\circ}\) ce qui fait \begin{align*}\cos^{2}\left(\frac{\theta + \phi}{2}\right) &= \cos^{2}\left(\frac{180^{\circ}}{2}\right)  \\ \\ &=\cos^{2}\left(90^{\circ}\right) \\ \\ &=0\end{align*}et \[K = \sqrt{\left(s-a\right)\left(s-b\right)\left(s-c\right)\left(s-d\right)-abcd\cos^{2}\left(\frac{\theta+\phi}{2}\right)}\]devient \[K = \sqrt{(s-a)(s-b)(s-c)(s-d)}\]toujours avec \[s = \frac{a+b+c+d}{2}\]

Formule de Héron

En bonus (bis) ! On peut obtenir facilement la formule de Héron à partir de celle de Brahmagupta. En effet, à partir d’un quadrilatère inscriptible, on fait tendre \(d\) vers \(0\). Le demi-périmètre \[s = \frac{a+b+c+d}{2}\] tend vers \[s = \frac{a+b+c}{2}\] et l’aire du quadrilatère \[K = \sqrt{(s-a)(s-b)(s-c)(s-d)}\]tend bien sûr vers \[K = \sqrt{s(s-a)(s-b)(s-c)}\]Voilà !

 

Références :

Niven, Ivan, Maxima and Minima Without Calculus, MAA, 1981

Alsina, Claudia and Roger B. Nelsen, A Cornucopia of Quadrilaterals, AMS/MAA Press, 2020

Johnson, Roger A, Advanced Euclidean Geometry, Dover, 2007