function [x_koz,status,k] = NewtonModszer(f, df, x0, eps, status)
  delta   = 10^(-10);
  x_old   = x0;
  elteres = 2*eps;                    % hogy elsõre lefusson a while ciklus
  k = 0;
  
  while elteres > eps 
    k = k+1;
    dfxold = df(x_old); 
    if abs(dfxold) < delta            % a derivált (majdnem) nulla 
      warning('newton: a derivált (majdnem) 0');
      status = 'zerus derivalt';
      x_koz  = x_old;
      return;
    endif
    x_new   = x_old - f(x_old)/dfxold;
    elteres = abs(x_new - x_old);     % a két közelítés eltérése
    x_old   = x_new;   
  endwhile
  
  x_koz  = x_new;
  status = 'konvergált';
endfunction