Abt. Wurzelbehandlung (Heron)
====================
Das Heron-Verfahren zum Quadratwurzelziehen, hier auf XProfan-nahe Genauigkeit.
Gruss
P.S.: Daß auch die Lösungen von Polynomgleichungen als Wurzeln bezeichnet werden, ist auf eine kapitale Fehlübersetzung zurückzuführen. Details hier!
Code
Windowtitle "Heron-Verfahren für das Wurzelziehen"
Windowstyle 24:Window 800,600:declare v!,w!,f!
set("decimals",17)
whileloop 0,46:v!=2^(&Loop-16):w!=heronsqrt(v!)
f!=w!-sqrt(v!):print v!,tab(30);w!,tab(59+(f!>=0));f!
endwhile :waitinput :END
proc heronsqrt :parameters a!:declare xn!,x!:xn!=2
repeat :x!=xn!:xn!=0.5*(x!+a!/x!)
until abs(xn!-x!)=0:return xn!
endproc
Alles anzeigen