Das Heron-Verfahren
> restart;
Iterative Bestimmung der Quadratwurzel von a
Die Iterationsvorschrift:
> heron:= x -> 1/2*(x + a/x);
> a:= 13;
Wahl eines Startwertes:
> x[alt]:= 4;
![x[alt] := 4](images-heron/heron4.gif){kind=small}
> x[neu]:= heron(x[alt]);
![x[neu] := 29/8](images-heron/heron5.gif){kind=small}
> while (abs(x[neu]-x[alt]) > 1/1000000) do
x[alt]:= x[neu]:
x[neu]:= heron(x[alt]):
print(x[neu]);
od:
x[neu]:= heron(x[alt]):
print(x[neu]);
od:
