Heron-Verfahren
Das Heron-Verfahren, Heronsche Näherungsverfahren oder Babylonische Wurzelziehen ist ein Rechenverfahren zur Berechnung einer Näherung der Quadratwurzel einer reellen Zahl <math>a>0</math>. Hierbei wird die Zahl <math>a</math> als Flächeninhalt eines Rechtecks aufgefasst (z. B. mit Seitenlängen <math>a</math> und <math>1</math>). Dieses Rechteck wird dann schrittweise in ein flächengleiches Quadrat transformiert, indem man in jedem Rechenschritt die längere Seite des vorherigen Rechtecks verkürzt und seine kürzere Seite so verlängert, so dass der Flächeninhalt <math>a</math> erhalten bleibt. Die verkürzte neue längere Seite berechnet sich dabei als Mittelwert der beiden Seiten des vorherigen Rechtecks (siehe Grafik rechts). Das Verfahren ist nach dem griechischen Mathematiker Heron von Alexandria benannt, der es in seinem Werk Metrika beschrieb. Allerdings wurde es schon über 1000 Jahre früher von den Babyloniern benutzt.
Im Gegensatz zum schriftlichen Wurzelziehen benötigt man keinen festgelegten, also korrekten, Ausgangswert. Zudem ist das Verfahren relativ robust gegen Rundungsfehler und konvergiert in der Regel schneller. Jedoch können Wurzeln mit dem Heronverfahren prinzipiell nur näherungsweise berechnet werden.
Geometrische Veranschaulichung und Grundidee
Dem Heron-Verfahren liegt die Idee zu Grunde, dass ein Quadrat mit Flächeninhalt <math>A</math> eine Seitenlänge von <math>\sqrt A</math> hat. Ausgangspunkt des Verfahrens ist ein beliebiges Rechteck mit Flächeninhalt <math>A</math>. Schritt für Schritt wird das Seitenverhältnis des Rechtecks so geändert, dass sich seine Form immer mehr der eines Quadrats annähert, während der Flächeninhalt gleich bleibt. Die Seitenlängen des Rechtecks sind die Näherungswerte für <math>\sqrt A</math>.
Im ersten Schritt wird eine beliebige Seitenlänge <math>x_0</math> für das Rechteck gewählt. Damit dieses den gewünschten Flächeninhalt hat, wird die zweite Seitenlänge mit der Formel
- <math>y_0 = \frac A {x_0}</math>
berechnet. Als Beispiel soll die Wurzel aus 9 berechnet werden. Für die eine Seitenlänge wird der Wert 9 gewählt, sodass sich die andere Seitenlänge zu 1 berechnet. Das erste Rechteck hat deshalb die folgende Form.
Die Ähnlichkeit dieses Rechteckes mit einem Quadrat ist gering. Das kommt auch dadurch zum Ausdruck, dass die Seitenlängen 1 und 9 sehr schlechte Näherungen für die Wurzel aus 9 sind.
Um eine bessere Annäherung an ein Quadrat zu erhalten, muss die lange Seite gekürzt und die kurze Seite verlängert werden. Als neue Länge der langen Seite wird der Mittelwert
- <math>x_1 = \frac{x_0 + y_0}{2}</math>
der beiden bisherigen Seitenlängen genommen. Die Länge der anderen Seite berechnet sich wie oben zu
- <math>y_1 = \frac A {x_1}</math>.
Im Beispiel ergibt sich als Mittelwert die Seitenlänge 5. Die dazugehörige kurze Seite hat eine Länge von 1,8.
Auch hier ist die Ähnlichkeit zu einem Quadrat noch gering. Allerdings ist das neue Rechteck im Vergleich zum vorhergehenden kompakter.
Der beschriebene Ablauf wird in jedem weiteren Schritt des Heron-Verfahrens wiederholt. Der Mittelwert der Seitenlängen eines Rechtecks entspricht der Länge der langen Seite des neuen Rechtecks und die Länge der kurzen Seite lässt sich daraus jeweils wie oben beschrieben berechnen. Im Beispiel entstehen so in den nächsten zwei Schritten die folgenden beiden Rechtecke.
Datei:Part 3 of a geometric example of Herons method.svg Datei:Part 4 of a geometric example of Herons method.svg
Das letzte Rechteck ist schon annähernd quadratisch. Die Seitenlänge 3,024 liegt entsprechend nah bei 3, dem exakten Wert von <math>\sqrt 9</math>.
Iterationsverfahren
Entwicklung aus der Grundidee, Rekursionsgleichung
Aus der geometrischen Grundidee erhält man ein allgemeines Iterationsverfahren zur näherungsweisen Berechnung der Wurzel einer reellen Zahl <math>a>0</math>:
Man geht von einem beliebigen Startwert <math>x_0 \neq 0</math> (idealerweise in der Nähe von <math>\sqrt a</math>) aus und setzt <math>y_0 = a/x_0</math>. Da in jedem Iterationsschritt die eine Seite durch den Mittelwert der beiden Seiten ersetzt wird und die andere Seite so angepasst wird, dass der Flächeninhalt unverändert bleibt, lautet die Iterationsvorschrift
- <math>x_{n+1} = \frac{x_n + y_n}{2} \quad \text{und} \quad y_{n+1}=\frac{a}{x_{n+1}} </math>.
Häufig wird die Iteration in einer Form geschrieben, in der nur noch die Variable <math>x</math> vorkommt. Dazu setzt man <math display="inline">y_{n}=\frac{a}{x_{n}} </math> in die Gleichung für <math>x_{n+1}</math> ein und erhält die Rekursionsgleichung
- <math>x_{n+1} = \frac{1}{2}\cdot \left (x_n + \frac{a}{x_n}\right)</math>,
die eine Folge <math>x_0, x_1, x_2, \ldots</math> („Heron-Folge“) von immer besseren Näherungen von <math>\sqrt a</math> liefert.
Alternativ kann diese Rekursionsgleichung auch aus dem Newton-Verfahren für die Nullstelle der quadratischen Funktion <math>f(x)=x^2-a</math> hergeleitet werden.
Beispiel
Mithilfe des Heron-Verfahrens soll die Wurzel aus 2 angenähert werden. Als Startwert wird <math>x_0=2</math> gewählt. Die ersten Glieder der Heron-Folge lauten
- <math>x_1=\frac{1}{2}\cdot\left(2 + \frac{2}{2}\right) =\frac{3}{2} = 1{,}5</math>
- <math>x_2=\frac{1}{2}\cdot\left(\frac{3}{2} + \frac{2}{\frac{3}{2}}\right) =\frac{17}{12} =1{,}41666666\dots </math>
- <math>x_3=\frac{1}{2}\cdot\left(\frac{17}{12} + \frac{2}{\frac{17}{12}}\right)=\frac{577}{408} =1{,}41421568\dots .</math>
Nach drei Iterationen ist die Näherung bereits auf fünf Nachkommastellen genau und die Abweichung vom wahren Wert <math>\sqrt 2 = 1{,}41421356\ldots</math> beträgt somit weniger als 0,0001 %.
Konvergenz
Die Heron-Folge <math>(x_n)</math> mit <math display="inline">x_{n} = \frac{1}{2}\left (x_{n-1} + \frac{a}{x_{n-1}}\right) </math> konvergiert für jeden Startwert <math>x_0 > 0</math> gegen <math>\sqrt a</math>. Somit kann jede Wurzel durch das Iterationsverfahren beliebig angenähert werden.
| Beweis |
|
Alle Glieder der Heron-Folge sind positiv. Man kann nun zeigen, dass sie auch alle wenigstens so groß wie <math>\sqrt a</math> sind und die Heron-Folge somit nach unten beschränkt ist. Dazu zeigt man für beliebiges <math>n \geq 1</math> die Ungleichung <math> x_n^2 - a \ge 0 </math>:
Weiter zeigt man, dass die Heron-Folge monoton fallend ist:
Aufgrund der Beschränktheit und Monotonie muss die Folge nach dem Monotoniekriterium gegen einen Grenzwert <math>x</math> konvergieren. Es bleibt noch zu zeigen, dass <math>x=\sqrt a</math>. Hierzu ist es zweckmäßig, die Folge <math>x_0^2, x_1^2,x_2^2, \ldots </math> zu betrachten, die gegen <math>x^2</math> konvergiert. Aus <math>x_n^2 \geq a > 0 </math> folgt <math>x^2\geq a > 0</math>, also ist <math>x \neq 0</math>. Nun lässt sich der Grenzwert <math>x</math> berechnen:
Hieraus erhält man durch elementare Termumformungen <math>x^2 = a</math>, woraus schließlich <math>x=\sqrt a</math> folgt. |
Der Startwert <math>x_0</math> der Iteration kann sogar, solange er von Null verschieden ist, beliebig festgesetzt werden, die Iteration konvergiert immer. Zu beachten ist, dass bei negativen Startwerten die Iteration gegen die negative Quadratwurzel konvergiert. Da sich das Heron-Verfahren aus dem Newtonschen Näherungsverfahren ableiten lässt und die zu berechnende Nullstelle einfach ist, ist die Konvergenzordnung 2. Die Zahl der richtigen Stellen wird mit jedem Schritt etwa verdoppelt.
Bei einem Startwert <math>x_0</math> in der Nähe von <math>\sqrt a</math> erhält man mithilfe des Heron-Verfahrens schnell gute Näherungswerte, wie das Beispiel verdeutlicht. Wenn die Anfangsnäherung jedoch schlecht ist, sind viele Schritte für eine gute Näherung nötig. Wenn zum Beispiel die Wurzel einer ganzen Zahl <math>a</math> mit 200 Binärstellen berechnet werden soll und man <math>x_0 = a</math> als Startwert verwendet, dann wird die Näherung mit jedem Schritt um etwa eine Binärstelle kürzer, d. h. erst nach etwa 100 Schritten hat die Näherung die richtige Länge von 100 Stellen. Danach reichen sechs bis sieben weitere Schritte (<math>\log_2 (100)</math>), um alle 100 Stellen vor dem Komma richtig zu berechnen. Es empfiehlt sich somit, einen möglichst genauen Startwert <math>x_0</math> zu bestimmen. Im Beispiel sollte man zuerst die Bitlänge <math>\lfloor \log_2(a) \rfloor + 1</math> von <math>a</math> ermitteln und einen Startwert mit der halben Länge verwenden.<ref group="A" name="numerisch">Startwert: Sofern der Ausgangswert bereits als Binärzahl (im Stellenwertsystem) vorliegt, kann einfach gezählt werden, an welcher Stelle <math>i</math> seine höchstwertige '1' steht; Startwert wird dann <math>1\cdot 2^{i/2}</math>. Sofern der Ausgangswert in (Binär-)Exponentialdarstellung vorliegt, kann als Startwert einfach der Exponent halbiert werden (um 1 Bit nach rechts schieben). Siehe auch Abschnitt Implementierung in Software</ref> Eine qualifizierte Schätzung für den Startwert erhält man aus der Taylorreihen-Entwicklung der binomischen Reihe um 1, deren zwei erste Glieder die Gleichung <math>x_0 = \tfrac{a+1}{2}</math> liefern.
Das Heron-Verfahren gehört zu den Fixpunktverfahren. Setzt man <math display="inline">\varphi(x) = \frac{1}{2}\cdot\left(x + \frac{a}{x}\right)</math>, so gilt für den Fixpunkt (der die Bedingung <math display="inline">\varphi(x) = x</math> erfüllt) <math display="inline">x^2 = a</math> mit der (positiven) Lösung <math display="inline">x = \sqrt{a}</math>.
Fehlerabschätzung
Für die Heron-Folge <math>(x_n)_ {n \ge 1}</math> gilt die Abschätzung
- <math>\frac{a}{x_n} \le \sqrt{a} \le x_n</math>,
und für den Fehler die Abschätzung
- <math>x_n-\sqrt{a} = \frac {1}{2x_{n-1}} \left( x_{n-1}-\sqrt{a} \right)^2 \le \frac {1}{2x_{n-1}} \left( x_{n-1}-\frac{a}{x_{n}} \right)^2 = \frac{1}{2x_{n-1} \cdot x_n^2} \left(x_{n-1} \cdot x_n-a \right)^2</math>.
Angewandt auf obiges Beispiel erhält man:
- <math>x_3-\sqrt 2 = \frac {1}{2x_2} \left( x_2-\sqrt 2 \right)^2 \le \frac{1}{2x_2 \cdot x_3^2} \left(x_2 \cdot x_3-2 \right)^2=0{,}00000212\dots</math>
Für den relativen Fehler
- <math>\varepsilon_n = \frac{x_n - \sqrt{a}}{\sqrt{a}}</math>
gilt die Rekursion
- <math>\varepsilon_{n+1} = \frac{\varepsilon_n^2}{2(1+\varepsilon_n)}</math>.
Die Folge der <math>\varepsilon_n</math> ist also bei gegebenem relativen Fehler <math>\varepsilon_0</math> der Startnäherung unabhängig von <math>a</math>.
Implementierung in Software
Das Verfahren eignet sich besonders gut zur Implementierung in Software, da nur Grundrechenarten benötigt werden, s. o. Es wird heute angesichts der breiten Verfügbarkeit numerischer Prozessorhardware aber nur noch selten benötigt.
Wenn dazu noch eine Gleitkommadarstellung mit einem Zweier-Exponenten benutzt wird, wird der Ansatz relativ einfach, als Beispiel wird die Wurzel aus 5 betrachtet und der relative Fehler zum Endwert <math>\frac{|x_i - x|}{x}</math> verfolgt:
- Zunächst wird von diesem Zweier-Exponenten eine gerade Anzahl abgespaltet, so dass als Exponent entweder eine 0 oder 1 übrig bleibt, die Zahl also auf das Intervall <math>[\tfrac{1}{2}, 2]</math> normalisiert wird. In diesem Intervall ist die Wurzelfunktion eine nur schwach gekrümmte Kurve, lässt sich also numerisch gut behandeln. Beispiel:<math>\sqrt{5} = \sqrt{4 \cdot 1{,}25} = 2 \cdot \sqrt{1{,}25} \approx 2 \cdot 1{,}118034 = 2{,}236068</math>, es wird also vorerst nur noch <math>a = 1{,}25</math> mit dem Ziel <math>x = 1{,}118</math> behandelt.
- Als Startwert für die eigentliche Iteration approximiert man diese Kurve durch eine noch einfachere, die sich direkt ohne Iteration berechnen lässt. Mit dieser Anfangsberechnung wird der Startwert ermittelt, mit dem die folgende Iteration begonnen wird. Man kann diese Kurve mehr oder weniger aufwendig ansetzen, mit den steigend komplizierteren Ansätzen unten lässt sich gegebenenfalls ein Iterationsschritt einsparen:
- eine einfache Konstante (beispielsweise 1),
Beispiel: <math>x_0 = 1</math>, relativer Fehler <math>1{,}1 \cdot 10^{-1} </math> - eine Gerade mit Steigung <math>\tfrac{1}{2}</math> und einer additiven Konstante von <math>\tfrac{1}{2}</math> als Vereinfachung des nachfolgenden Falls
Beispiel: <math>x_0 =\tfrac{1}{2} + \tfrac{1{,}25}{2} = 1{,}125</math>, relativer Fehler <math>6{,}2 \cdot 10^{-3} </math> - eine Gerade mit Steigung <math>\tfrac{1}{2}</math> und einer additiven, optimierten Konstante von <math>\left( 2\cdot
- eine einfache Konstante (beispielsweise 1),
\sqrt[4]{2} - \sqrt{2} \right)^2 / 2 \approx 0{,}4648415</math>,
Beispiel: <math>x_0 = \tfrac{0{,}929683}{2} + \tfrac{1{,}25}{2} \approx 1{,}089841</math>, relativer Fehler <math>2{,}5 \cdot 10^{-2}</math>.
- eine Gerade mit optimierter Steigung und einer additiven Konstante hier nicht näher betrachtet.
- Ausgehend von dem so ermittelten Startwert <math>x_0</math> führt man eine feste Anzahl von Iterationsschritten durch. Die nötige Anzahl, um die gewünschte Genauigkeit zu erreichen, lässt sich dank der obigen Fehlerabschätzung als Worst Case innerhalb des Startintervalls direkt ausrechnen. Bei 32 Bits Mantisse und dem mittleren Startansatz braucht man beispielsweise nur drei Schritte. Diese fest gewählte Anzahl erspart wesentlich aufwendigere Abfragen auf Erreichung der Genauigkeit. Der Ersatz der genannten Konstanten durch die Zahl 1,0 ändert daran nichts. Auch der noch kompliziertere Ansatz brächte zumindest bei dieser Genauigkeit keine Einsparung eines weiteren Iterationsschritts. Bei höheren Genauigkeitsanforderungen kann das anders aussehen.
Beispiel mit drei Schritten nach Ansatz 1 (Konstante 1, mit den anderen Ansätzen konvergiert es noch einen Schritt schneller):
<math> x_1 = \tfrac{x_0 + \tfrac{a}{x_0}}{2}x_1 = \tfrac{x_0 + \tfrac{1{,}25}{x_0}}{2} = \tfrac{1 + \tfrac{1{,}25}{1}}{2} = 1{,}125 </math>, relativer Fehler <math>6{,}2 \cdot 10^{-3}</math>
<math> x_2 = \tfrac{x_1 + \tfrac{a}{x_1}}{2} = \tfrac{1{,}125 + \tfrac{1{,}25}{1{,}125}}{2} \approx 1{,}118056 </math>, relativer Fehler <math>2{,}0 \cdot 10^{-5}</math>
<math> x_3 = \tfrac{x_2 + \tfrac{a}{x_2}}{2} = \tfrac{1{,}118056 + \tfrac{1{,}25}{1{,}118056}}{2} \approx 1{,}118034 </math>, relativer Fehler kleiner als <math>10^{-6}</math>
Man sieht die Wirkung der quadratischen Konvergenz, dass sich der relative Fehler von Schritt zu Schritt jeweils quadriert oder die Anzahl gültiger Stellen bzw. der negative Fehlerexponent etwa verdoppelt. - Zum Schluss wird der Exponent restauriert, indem man die Hälfte des im ersten Schritt abgespalteten Werts wieder hinzufügt.
Beispiel: <math> 2 \cdot x_3 = x_2 + \tfrac{a}{x_2} = 1{,}118056 + \tfrac{1{,}25}{1{,}118056} \approx 2{,}236068 </math>.
Verallgemeinerung des Verfahrens
Näherung höherer Wurzeln
Dieses Verfahren lässt sich verallgemeinern, um die k-te Wurzel einer Zahl <math>a>0</math> näherungsweise zu berechnen. Dabei wird das Newton-Verfahren zur Bestimmung der positiven Nullstelle <math>\sqrt[k]{a}</math> der Funktion <math>f(x)=x^k-a</math> angewandt. Mit <math>f'(x)=kx^{k-1}</math> folgt aus der Rekursionsformel <math display="inline">x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)} </math> des Newton-Verfahrens die Iterationsvorschrift
- <math>x_{n+1}=\frac1k\left((k-1)x_n+\frac{a}{x_n^{k-1}}\right).</math>
Beispielsweise lautet die Rekursionsformel zur Berechnung der Kubikwurzel
- <math>x_{n+1}=\frac13\left(2x_n+\frac{a}{x_n^2}\right).</math>
Hier muss die Folge mit einem geeigneten Startwert <math>x_0</math> für den gesuchten Wert von <math>\sqrt[k]{a}</math> gestartet werden.
Je größer <math>k</math> ist, desto mehr Schritte werden benötigt, um die Wurzel genau zu berechnen. Für ganzzahliges positives <math>k</math> gelten die gleichen Konvergenzaussagen wie oben für <math>k=2.</math>
Näherung des Kehrwerts
Verfahren und Ersetzen der Division durch Verwendung des Kehrwerts
Für <math>k=-1</math> erhält man ein Verfahren, mit dem (ohne Verwendung der Division!) der Kehrwert <math>a^{-1}=1/a</math> näherungsweise errechnet werden kann:
- <math>x_{n+1}=\frac{(-1-1)x_n^{-1}+a}{(-1)x_n^{-1-1}}=2x_n-ax_n^2=(2-ax_n) \cdot x_n.</math>
Dieses Verfahren konvergiert für alle <math>x_0\in\left( 0,2/a \right) </math> quadratisch gegen <math>1/a.</math>
Die Iteration ermöglichte für erste Computer ohne eingebaute Division die Zurückführung dieser Operation auf Multiplikation und Subtraktion. Die Division von zwei Zahlen wurde so ausgeführt, dass der Kehrwert des Nenners bestimmt wurde und mit dem Zähler multipliziert wurde.
Beispiel
Es soll <math>1/3</math> näherungsweise berechnet werden mit dem Startwert <math>x_0=\tfrac12 < \tfrac23</math>:
- <math>x_1=\left(2-3\cdot \frac 12 \right) \cdot \frac 12=\frac14=0{,}25,</math>
- <math>x_2=\left(2-3\cdot \frac 14 \right) \cdot \frac 14=\frac5{16}=0{,}3125,</math>
- <math>x_3=\left(2-3\cdot \frac 5{16} \right) \cdot \frac 5{16}=\frac{85}{256}=0{,}33203125.</math>
Für den Startwert <math>x_0=\frac23</math> erhält man
- <math>x_1=\left(2-3\cdot \frac 23 \right) \cdot \frac 23=0,</math>
- <math>x_2=\left(2-3\cdot 0\right) \cdot 0=0,</math>
somit keine Konvergenz gegen den gesuchten Wert von <math>\tfrac13.</math>
Geschichte
Das Verfahren war in Mesopotamien bereits zur Zeit des babylonischen Königs Hammurapi I. (ca. 1750 v. Chr.) bekannt.<ref>Jochen Ziegenbalg, Oliver Ziegenbalg, Bernd Ziegenbalg: Algorithmen: Von Hammurapi bis Gödel. Springer Spektrum 2016, S. 55 ({{#if: VsdeCwAAQBAJ
| {{#if: {{#if: ||1}} {{#if: VsdeCwAAQBAJ ||1}}
| <0|&pg={{#if:|RA{{{Band}}}-}}PA55|&pg=55}}{{#if:|&q=}}#v=onepage|{{#if:|&pg=|}}{{#if:|&q=}}}}{{#if:|q=%7B%7B%7BSuchbegriff%7D%7D%7D}}|{{#if:|q=%7B%7B%7BSuchbegriff%7D%7D%7D}}}} {{#if:Auszug (Google)|{{#invoke:WLink|getEscapedTitle|Auszug (Google)}}|eingeschränkte Vorschau}}{{#if:ja|| in der Google-Buchsuche}}{{#ifeq:|US|-USA}}{{#if: VsdeCwAAQBAJ |{{#invoke: Vorlage:GoogleBook|fine |id=VsdeCwAAQBAJ |errN=Parameter „BuchID“ hat falsche Länge |errC=Parameter „BuchID“ enthält ungültige Zeichen |errH=# in der „BuchID“ |errP=Parameterzuweisungen in der „BuchID“ |class=editoronly |cat={{#ifeq: 0 | 0 | Wikipedia:Vorlagenfehler/Vorlage:Google Buch}} |template= Vorlage:Google Buch}}
}}
| Es darf nur genau einer der beiden Parameter „Suchbegriff“ oder „BuchID“ ausgefüllt werden. Bitte beachte die in der Vorlage:Google Buch befindliche Dokumentation und prüfe die verwendeten Parameter.{{#ifeq: 0 | 0 | }}}}
| Es muss mindestens einer der beiden Parameter „Suchbegriff“ oder „BuchID“ ausgefüllt werden. Bitte beachte die in der Vorlage:Google Buch befindliche Dokumentation und prüfe die verwendeten Parameter.{{#ifeq: 0 | 0 | }}}}{{#invoke:TemplatePar|check
|all=
|opt= Suchbegriff= BuchID= Seite= Band= SeitenID= Hervorhebung= Linktext= Land= KeinText=
|cat= {{#ifeq: 0 | 0 | Wikipedia:Vorlagenfehler/Vorlage:Google Buch}}
|template= Vorlage:Google Buch
|format=
}}{{#if:Auszug (Google)|{{#if:{{#invoke:WLink|isBracketedLink|Auszug (Google)}}|}}}}).</ref> Um 100 n. Chr. wurde es von Heron von Alexandria im ersten Buch seines Werkes Metrika beschrieben.<ref>{{#if: |{{{author}}}: |John J. O’Connor, Edmund F. Robertson: }}[https://mathshistory.st-andrews.ac.uk/{{#switch:
| bio = Biographies
| cur = Curves
| ex = Extras
| ht = HistTopics
| misc = Miscellaneous
| soc = Societies
| #default = Biographies
}}/Heron/ Heron-Verfahren.] In: {{#invoke:Vorlage:lang|flat}}{{#if: |, {{#invoke:DateTime|format|{{{datum}}}|T._Monat JJJJ}}}} (englisch).</ref>
Literatur
- Jochen Ziegenbalg, Oliver Ziegenbalg, Bernd Ziegenbalg: Algorithmen: Von Hammurapi bis Gödel. 4. Auflage. Springer Spektrum, Wiesbaden 2016, ISBN 978-3-658-12362-8, S. 55–60.
- David Fowler, Eleanor Robson: Square Root Approximations in Old Babylonian Mathematics. In: Historia Mathematica, Band 25, Nr. 4, 1998, S. 366–378 (sciencedirect.com).
- Hans Rudolf Schwarz, Norbert Köckler: Numerische Mathematik. 7. Auflage. Vieweg+Teubner, Wiesbaden 2009, ISBN 978-3-8348-0683-3, S. 192–194.
- Dietmar Herrmann: Die antike Mathematik. 3. Auflage. Springer Spektrum, Berlin 2024, ISBN 978-3-662-68477-1, S. 337–339.
- Konrad Königsberger: Analysis 1. 6. Auflage. Springer, Berlin / Heidelberg 2004, ISBN 978-3-540-40371-5, S. 48–49.
Weblinks
|X|x= |0|-= |S|s= – Sammlung von Bildern |1|= – Sammlung von Bildern{{#if:
| {{#switch: {{#invoke:TemplUtl|faculty|1}}/{{#invoke:TemplUtl|faculty|1}}
|1/= und Videos
|1/1=, Videos und Audiodateien
|/1= und Audiodateien}}
| , Videos und Audiodateien
}}
|#default= – }}{{#if: Heron's method
| {{#ifeq: {{#invoke:Str|left|heron's method|9}}
| category:
| FEHLER: Ohne Category: angeben!}}}}Vorlage:Wikidata-Registrierung
- Das Heron-Verfahren zur Wurzelberechnung auf arndt-bruenner.de (Erläuterung und Beispielrechner)
- Interaktive Illustration in GeoGebra (Papp)
- Interaktive Illustration in GeoGebra (Leo)
Anmerkungen
<references group="A" />
Einzelnachweise
<references />
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