function x = newtonfler(start,F,J)
  % Funksjonsprosedyre som utf?rer Newtons metode p? en 
  % vektorvaluert funksjon F med Jacobi-matrise J
  % Innparametre:
  % 
  % start: s?ylevektor som inneholder startpunktet for iterasjonen 
  % F: funksjon som regner ut funksjonsverdiene
  % J: Funksjon som regner ut Jacobideterminanten
  maxit=30; 
  tol=0.0000001; 
  diff=tol+1; 
  x=start; 
  for t=1:maxit 
    xny=x-J(x)\F(x); 
    diff=norm(xny-x); 
    x=xny; 
    fval=F(x);
    fprintf('itnr=%2d x=[%13.10f,%13.10f] F(x)=[%13.10f,%13.10f]\n',... 
             t,x(1),x(2),fval(1),fval(2)) 
    if diff<tol 
      break; 
    end %if 
  end %for 
  % Kalles f.eks med kommandoen: 
  % >> x=newtonfler([-7;1],F,J)