    function [x,it]=newtsys(f,J,x0,tol)
    % Lös ickelineära systemet f(x)=0 med Newton-Raphsons 
    % metod och startvärdet x0
    % Avbryt när normen av korrektionen är mindre än tol
    it=0; x=x0;
    korr=1e10; % Alltid minst en iteration 
    while korr>tol,
      h=feval(J,x)\feval(f,x);;
      korr=norm(h);
      x=x-h;
      it=it+1;
%      disp(x'),
    end

