function y=f(x) y=x-cos(x) endfunction function y=fp(x) y=1+sin(x) endfunction x0=-1.0; n=10; x(1)=x0; for i=1:n x(i+1) = x(i) - f(x(i)) / fp(x(i)); end disp(x)