function [x,it] = iter(fun,x0,tol)
% [x,it] = iter(fun,x0,tol)
%
% funkce nalezne pevny bod 
%      f(x) = x
% prostou iteracni metodou
%
% x  ... pevny bod
% it ... pocet iteraci
%
% fun ... f(x) zadana ve stringu (promenna x)
% x0  ... pocatecni aproximace
% tol ... pozadovana absolutni presnost


it = 0;
x=x0;
xn = feval(fun,x);

while norm(xn-x) >= tol 

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

end

x = xn;

% end