function [x] = my_new_alg_eq(c,x0,tol, maxit)
%metodo delle tangenti o di Newton per la risoluzione
%di equazioni algebriche
%x0=-max(1, sum(abs(c(2:end)/c(1))));
disp('punto iniziale') 
x0
err=inf;
it=0;
n=length(c);
c1=[n-1:-1:1].*c(1:n-1);
while(err>tol & it<=maxit)
    x=x0-my_rh(c,x0)/my_rh(c1,x0);
    err=abs(x-x0);
    disp('**********************');
%    abs(x-appv)
    it=it+1;
    it
    x0=x;
end
end