%Aufgabe 3b Serie 5 marco Weber

function[x_n,x_hist,k]=my_damp_newton(F,Df,x0,tol,maxit,alpha)
x_o=x0
x_hist=x0;
for k=1:1:maxit
    
    J=Df(x_o);
    fun=F(x_o);
    x_n = x_o - alpha*(J\fun)
    
    x_hist = [x_hist,x_n];

    if(norm(x_n-x_o)<tol)
        return
    end
    
x_o = x_n;
end

%alles ähnlich gemach wie bei Aufgabe 2
