function x = newton(start,f,df) 
  % Funksjon som utf?rer Newtons metode p? en funksjon 
  % f med derivert df spesifisert via de anonyme funksjonene f og df
  % Innparametre: 
  % start: startpunkt for iterasjonen 
  % f: anonym funksjon som returnerer funksjonsverdiene
  % df: anonym funksjon som returnerer verdiene til den deriverte
  maxit=30; 
  tol=0.0000001; 
  diff=tol+1; 
  x=start; 
  for t=1:maxit 
    xny=x-f(x)/df(x); 
    diff=abs(xny-x); 
    x=xny; 
    fprintf('itnr=%2d x=%15.12f f(x)=%15.12f\n',t,x,f(x)) 
    if diff<tol 
      break; 
    end %if 
  end %for