function list = newt(f,df,x0,n)


x=[x0];

for k=1:(n-1)
    y=x(k)-f(x(k))/df(x(k));
    x=[x,y];    
end


list = x';