% NEWTON2D je Newtonova iteracija v 2d
% Zk=newton2d(fun,jac,Z0,k,tol)
% fun je funkcija, jac je Jacobijeva matrika, Z0 zacetni priblizek,
% k je stevilo iteracij
% tol je toleranca med dvema zaporednima priblizkoma
function [Zk,it]=newton2d(fun,jac,Z0,k,tol);
	Zk=Z0;
        napaka=tol+1;
        it=0;
        while (it<k)&(napaka>tol) 
	  Zk=feval(jac,Z0)\(-feval(fun,Z0))+Z0;
          it=it+1;
          napaka=norm(Zk-Z0);
          Z0=Zk;
       end;
