a = input('Aus welcher Zahl soll die Wurzel bestimmt werden :> ') ;
x = a ;
for zaehler = 1 : 16
    zaehler ;
    x_neu = 1/2 * ( x + a/x ) ;
    x = x_neu ;
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) 