function w = wurzel(x, n)

x0 = 1;
for i = 1:100
    x1 = x0 - (x0^n - x) / x0^(n-1) / n;
    if abs(x1 - x0) < 10^-15
        break
    end
    x0 = x1;
end
w = x1;
disp(i)

end  % endfunction