first lab / NewtonMeth (1)

function x = NewtonMeth(func, dfunc, x0,n);
x=x0;
for i=1:n
y=func(x);
dy = dfunc(x);
xnew= x-y/dy;
if abs(xnew-x)<1e-6;

disp(xnew);
error_newton = abs (func(x));
disp(['Невязка метода Ньютона = ', num2str(error_newton)])
break
end
x=xnew;

end

end