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