a = input('Aus welcher positiven Zahl soll die Wurzel bestimmt werden :> ') ;
x = a ;
x_neu = 2*a ; % Damit im while die Bedingung erfüllt ist
epsilon = 0.0001 ; 
zaehler = 0 
while abs(x_neu - x)/abs(x_neu) > epsilon
    zaehler = zaehler + 1 ;
    x = x_neu
    x_neu = 1/2 * ( x + a/x ) 
end
fprintf( 'Die Wurzel aus %9.3f ist = %12.8f\n', a, x ) 
rel_fehler = abs( a - x*x )/ abs(a) ;
fprintf( ' Nach %2d Iterationen ist der relative Fehler : %12.8f\n', zaehler, rel_fehler) 