% 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=newton2d(fun,jac,Z0,k,tol);
	for i=1:k
		Zk=feval(jac,Z0)\(-feval(fun,Z0))+Z0;
		Z0=Zk;
	end
