function y=newton(effeeeffeprimo, x0,k)
y=zeros(k,1);
z=x0;
for i=1:k
    [f fprimo]=effeeeffeprimo(z);
    z=z-f/fprimo;
    y(i)=z;
end
end