    function [x,it]=newton(f,fprim,x0,tol)
    % Lös f(x)=0 med noggrannhet tol med Newton-Raphsons metod 
    % och startvärdet x0
    % Iterera tills korrektionen är mindre än tol
    it=0; 
    korr=feval(f,x0)/feval(fprim,x0);
    while abs(korr)>tol,
      x=x0-korr
      x0=x;
      it=it+1;
      korr=feval(f,x0)/feval(fprim,x0);
    end
