function [X,it]=newton(F,JF,X0,tol,maxit)
%NEWTON resuje sistem nelinearnih enacb z Newtonovo metodo
%[X,it]=NEWTON(F,JF,X0,tol,maxit)
%X je priblizek za resitev, it stevilo potrebnih iteracij
%F je siste, JF Jacobian sistema
%X0 zacetni priblizek, tol toleranca za konec
%maxit maksimalno stevilo ponovitev

it=0;
X=X0;
Y=feval(F,X0);;
napaka=norm(Y,'inf');
while (napaka>tol) & (it<maxit)
	X=X-feval(JF,X)\Y;
	Y=feval(F,X);
	it=it+1;
	napaka=norm(Y,'inf');
end
