function [x,it] = newton_rov(fun,dfun,x0,tol)

% [x,it] = newton_rov(fun,dfun,x0,tol)
%
% funkce nalezne koren rovnice 
%      f(x) = 0
% Newtonovou metodou
%
% x  ... koren
% it ... pocet iteraci
%
% fun ... f(x) zadana ve stringu 
% dfun .. derivace f(x) zadana ve stringu 
% x0  ... pocatecni aproximace
% tol ... pozadovana absolutni presnost

x=x0;
xn = x - feval(fun,x)/feval(dfun,x);
it = 1;

while norm(xn-x) >= tol 

  it = it + 1;
  x = xn ;
  xn = x - feval(fun,x)/feval(dfun,x);

end

x = xn;

% end