func     = inline('x^2-2'); % Definerer f(x)
derivert = inline('2*x');   % Definerer f'(x)

NullTol = 1e-10;
x = 3;
Nmax = 100;   % Maksimalt antall iterasjoner
err  = 1.0;   % Dummy verdi for første iterasjon
Nint = 0;     % Antall iterasjoner utført
Tol  = 1e-12; % Toleranse for f'(x)

fprintf('Nit              x                f(x)\n',Nint,x(end),func(x(end)));

while ((Nint < Nmax) & (err > Tol)),
  fprintf('%2i\t %5.15f \t%0.5g\n',Nint,x(end),func(x(end)));
  x(end+1) = NewtonsMetode(x(end),func,derivert,NullTol);
  err = abs(x(end) - x(end-1));
  Nint = Nint + 1;
end