Voir la source de Méthode de Héron

{{Confusion|Formule de Héron}}
En [[mathématiques]], la '''méthode de Héron''' ou '''méthode babylonienne''' est une méthode efficace d'[[extraction de racine carrée]], c’est-à-dire de résolution de l'équation ''x''{{2}} = ''a'', avec ''a'' positif. Elle porte le nom du mathématicien [[Héron d'Alexandrie]], qui l'expose dans le tome I de son ouvrage ''Metrica'' (''Les métriques''), découvert seulement en 1896<ref>Cf. ''[[Encyclopædia Universalis]]'', édition 2008, ''Thésaurus - tome III'', {{p.|2517}}.</ref> mais certains calculs antérieurs, notamment égyptiens<ref>{{Ouvrage|langue = français|auteur1 = Marianne Michel|titre = Les mathématiques de l'Égypte ancienne. Numération, métrologie, arithmétique, géométrie et autres problèmes|lieu = |éditeur = Safran Bxl|année = 2014|pages totales = 608|isbn = 978-2-87457040-7|lire en ligne = |passage = }}.</ref>, semblent prouver que la méthode est plus ancienne.

Héron expose ainsi sa méthode dans le problème 8 du tome I des ''Métriques''. Il détaille initialement une méthode pour calculer l'aire d'un triangle en connaissant ses trois côtés ({{cf.}} [[formule de Héron]]), en prenant pour exemple un triangle de côtés 7, 8 et 9 unités. Il obtient alors le nombre 720 comme résultat intermédiaire, dont il doit calculer la racine carrée pour aboutir au résultat final. Il propose alors la méthode de calcul suivante<ref>Citation extraite de {{Chapitre|auteur=Bernard Vitrac|titre=Euclide et Héron : Deux approches de l'enseignement des mathématiques dans l'Antiquité ?|auteurs ouvrage=Gilbert Argoud|titre ouvrage= Science et vie intellectuelle à Alexandrie (Ier-IIIe siècle après J. C.)|lieu = |éditeur = Publications de l'Université de Saint-Étienne|année = 1995|pages totales = 226|isbn = 2-86272-058-5 |lire en ligne = https://hal.archives-ouvertes.fr/hal-00175155|page=121-145}}.</ref> :
{{citation bloc|Puisque alors les 720 n'ont pas le côté exprimable, nous prendrons le côté avec une très petite différence ainsi.
Puisque le carré le plus voisin de 720 est 729 et il a 27 comme côté, divise les 720 par
le 27 : il en résulte 26 et deux tiers.
Ajoute les 27 : il en résulte 53 et deux tiers.
De ceux-ci la moitié : il en résulte 26 2' 3'.
Le côté approché de 720 sera donc 26 2' 3'.
En effet 26 2' 3' par eux-mêmes : il en résulte 720 36',
de sorte que la différence est une {{36e|part}} d'unité.
Et si nous voulons que la différence se produise par une part plus petite que le 36', au
lieu de 729, nous placerons les 720 et 36' maintenant trouvés et, en faisant les mêmes choses, nous trouverons la différence qui en résulte inférieure, de beaucoup, au 36'.|[[Héron d'Alexandrie]]|''Metrica'', tome I, 8}}
== Exposé de la méthode ==
=== Approche géométrique ===
[[image:Heronsmethod.svg|thumb|right|Les rectangles ont même aire. Chaque rectangle a pour longueur la moyenne des dimensions du rectangle précédent.]]
Il est intéressant de mettre en évidence le principe géométrique sous-jacent à la méthode. Chez les mathématiciens grecs, extraire la racine carrée de {{math|''a''}} revient à trouver un carré dont l'aire soit {{math|''a''}}. En prenant un rectangle de côté arbitraire {{math|''x''}} et de même aire, il est nécessaire que l'autre côté ait pour longueur {{math|''a''/''x''}}. Pour le rendre « moins rectangle », il suffit de considérer un nouveau rectangle dont la longueur est la [[moyenne arithmétique]] des deux côtés précédents soit
{{Retrait|<math>\frac{x+a/x}2</math>}}
et dont l'aire reste {{math|''a''}}.

En réitérant infiniment le processus, le rectangle se transforme petit à petit en un carré de même aire. Cette constatation est à la base de la méthode de Héron.

=== Principe ===
Pour déterminer la racine carrée du nombre (positif) {{math|''a''}}, il convient dès lors de considérer la [[suite définie par récurrence]] de la façon suivante :
{{Retrait|<math>\forall n\in\N\quad x_{n+1}=\frac{x_n+\frac{a}{x_n}}2,</math>}}
de premier terme <math>x_0>0</math> choisi si possible « assez proche » de {{sqrt|''a''}}, en général la [[partie entière]] de {{sqrt|''a''}}.

La suite ainsi obtenue est une suite décroissante à partir du second terme, convergeant vers {{sqrt|''a''}}.
{{Démonstration|
En premier lieu, il est immédiat que, puisque <math>x_0 >0</math> et <math>a>0,</math> alors <math>x_n> 0.</math> De plus, pour tout <math>n\in\N</math> :
:<math>x_{n+1}^2-a=\left(\frac{x_n^2+a}{2x_n}\right)^2-a=\frac{x_n^4+a^2+2a x_n^2}{4x_n^2}-a=\frac{x_n^4+a^2-2a x_n^2}{4x_n^2}=\left(\frac{x_n^2-a}{2x_n}\right)^2 \ge0,</math>
ce qui implique :
:<math>\forall n\in\N^*\quad x_n\ge\sqrt a.</math>
Par ailleurs pour tout <math>n \in\N^*</math> on a :
:<math>x_{n+1}-x_n=\frac{x_n^2+a}{2x_n}-x_n=\frac{a-x_n^2}{2x_n},</math>
or puisque <math>x_n\ge\sqrt a</math>, on a <math>a-x_n^2 \leq 0</math> ce qui implique <math>\forall n\in\N^*\quad x_{n+1}\le x_n.</math>
La suite <math>(x_n)_{n\in\N^*}</math> est donc décroissante positive donc elle converge vers une limite <math>l\ge0</math> qui vérifie <math>l=\frac12\left(l+\frac al\right)  </math>, et donc <math> l= \sqrt a </math>.
}}
Si le premier terme de la suite est un [[nombre rationnel]], il est clair que tous les termes successifs seront des nombres rationnels, ce qui permet d'approcher un [[nombre irrationnel]] tel que la [[racine carrée de deux]] par une suite de rationnels.

=== Convergence ===
Il est par ailleurs facile de vérifier que la convergence en est [[Convergence quadratique|quadratique]] : l'écart entre chaque terme et la limite {{sqrt|''a''}} évolue comme le carré de l'écart précédent, en effet pour tout ''n'' > 0 :
: <math>x_{n+1}-\sqrt a=\frac{(x_n-\sqrt a)^2}{2x_n},</math>
soit puisque <math>x_n\ge \sqrt a</math> :
:<math>\left|x_{n+1}-\sqrt a\right|\le\frac{\left |x_n-\sqrt a\right|^2}{2\sqrt a},</math>
ce qui correspond bien à la définition de la convergence quadratique, c’est-à-dire que le nombre de décimales exactes double à chaque itération.

L'algorithme nécessite à chaque étape de faire une division, qui elle-même requiert une suite d'opérations d'autant plus longue que la précision demandée est importante. Néanmoins, l'algorithme est robuste, il supporte bien quelques approximations (et même quelques erreurs, dont l'effet sera de retarder l'obtention du résultat mais n'empêchera pas de l'obtenir), ce qui permet de se contenter de divisions (pas trop) fausses, au moins au début.

Du fait de sa convergence rapide, la méthode de Héron permet d'obtenir une bonne approximation de la valeur de {{sqrt|''a''}} même après peu d'étapes de calcul.

'''Exemple : calcul de {{sqrt|2}}'''
Soit <math>x_0 = 1 </math>, il vient successivement :
:<math>x_1=\frac{1+2}2=\frac32=1{,}5,</math>
:<math>x_2=\frac{3/2+4/3}2=\frac{17}{12}=1{,}4166\ldots,</math>
:<math>x_3=\frac{17/12+24/17}2=\frac{577}{408}=1{,}4142156\ldots,</math>
:<math>x_4=\frac{577/408+816/577}2=\frac{665857}{470832}=1{,}4142135623745\ldots,</math>
:<math>x_5=\frac{886731088897}{627013566048}=1{,}4142135623730950488016896\ldots,</math>
or en comparant avec la valeur exacte {{sqrt|2}} = 1,4142135623730950488016887…, on constate bien la convergence quadratique (2 décimales exactes au deuxième calcul, 5 au troisième, 11 au quatrième, 23 au cinquième…). En seulement trois étapes, la précision relative sur la valeur de {{sqrt|2}} est déjà de 10{{exp|–6}}, ce qui est excellent, et de moins de 10{{exp|–12}} en quatre étapes. De fait, une des principales problématiques est de choisir une « bonne » valeur pour <math>x_0</math>, idéalement l'entier dont le carré est le plus proche de {{math|''a''}}, ce que suggérait d'ailleurs Héron lui-même dans la partie des ''Metrica'' consacrée à cette question.

== Généralisation de la méthode ==
=== Racine ''n''-ième d'un nombre ===
{{Article détaillé|Algorithme de calcul de la racine n-ième}}
Une méthode analogue existe pour extraire la [[Racine d'un nombre|racine ''n''-ième d'un nombre]] ''A'', il convient de considérer alors la suite de terme général
<math>\forall k\in\N\quad x_{k+1} = \frac{1}{n} \left[{(n-1)x_k +\frac{A}{x_k^{n-1}}}\right]</math>.

L'idée géométrique sous-jacente est la même, puisque déterminer la racine ''n''-ième d'un nombre ''A'' consiste à trouver le côté d'un [[hypercube]] dont le « volume » est ''A''. La suite considérée revient à partir d'un (hyper)parallélépipède à ''n'' dimensions dont ''(n-1)'' côtés sont égaux, le dernier étant ajusté de façon à obtenir un volume égal à ''A''.

=== Lien avec la méthode de Newton ===
La méthode de Héron est un cas particulier de la [[méthode de Newton]]. En effet, dans la méthode de Newton, il s'agit de trouver un zéro d'une fonction ''f'' en utilisant la récurrence suivante :
: <math>x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)}.</math>
En prenant
: <math>f(x) = x^2 - a,</math>
la récurrence devient
: <math>x_{n+1}= x_n - \frac{x_n^2-a}{2x_n} = \frac{x_n^2+a}{2x_n} = \frac{x_n + \frac a{x_n}}2.</math>

== Notes et références ==
{{Références}}

== Voir aussi ==
* [[Racine carrée entière]]
* [[Extraction de racine carrée par la méthode du goutte à goutte|Méthode du goutte à goutte]]

{{Palette|Méthodes de résolution}}

{{Portail|analyse}}

[[Catégorie:Histoire de l'analyse]]
[[Catégorie:Algorithme de recherche d'un zéro d'une fonction|Héron]]
[[Catégorie:Racine carrée]]

[[es:Cálculo de la raíz cuadrada#Algoritmo babilónico]]
[[it:Metodi per il calcolo della radice quadrata#Metodo babilonese]]
[[pl:Metody obliczania pierwiastka kwadratowego#Metoda babilońska]]