function Newton
fprintf('\牛顿迭代法\n');
x(1)=1;
N0=100;
eps=0.00001;
fprintf('k x(k) f(x(k)) \n' )
for k=1:N0
x(k+1)=x(k)-f(x(k))/dif(x(k));
fprintf('%d %2.10f %2.10f \n',k ,x(k),f(x(k)));
if abs(x(k+1)-x(k))<eps
fprintf('success!');
break;
end
end
function y=f(x)
y=x^3-2;
function y=dif(x)
y=3*x^2;