function[x_new,x_hist,iter]=damped_newton(f,Df,x0,tol,maxit,alpha)
x_old=x0;
x_hist=x0;
for iter=1:maxit
    Delta=%!! Newtonkorrektur berechnen
    x_new =%!! Schritt mit Daempfungsfaktor alpha
    x_hist = [x_hist,x_new];
    if norm(x_new-x_old)<tol*norm(x_new)+tol
        return
    end
    x_old = x_new;
end

