Iterative Berechnung der Quadratwurzel

program Wurzel(input, output);

    {Berechnet eine Wurzel}
    
    var
        x, a, eps: double;
        
    begin {Wurzel}
        writeln('Berechnet die Quadratwurzel einer positiven reellen Zahl.');
        writeln;
        write('Argument: '); readln(a);
        write('Absolute Genauigkeit: ');
        readln(eps);

        x := a;
        while abs(x * x - a) >= eps do begin
            x := (x + a/x) / 2
        end; {while}
        
        writeln;
        writeln('Die Wurzel aus ', a:7:3, ' ist: ', x:7:3)
    end. {Wurzel}