function [ x ] = newton1( f, f1, x0, n )
%Newton iteration

x = x0;

for i = 1:n
    x = x - f(x) / f1(x);
end


end