function x = dnewton(x0, F, tol, maxit)
%x = dnewton(x0, F, tol, maxit) vrne resitev 
%enacbe F(x) = 0, ki jo poisce z diskretizirano Newtonovo 
%metodo. x0 je zacetni priblizek, F funkcija, tol 
%toleranca in maxit najvecje st. iteracij.

%velikost podatkov
n = length(x0);

%standardni bazni vektorji R^n
e = eye(n);

%izberemo korak
delta = sqrt(tol);

for k = 1:maxit
	F0 = feval(F, x0);
	%poracunamo priblizek za JF...
	for j = 1:n
		ajac(:,j) = (feval(F, x0 + delta*e(:,j)) - F0)/delta;
	end
	%... in izvedemo korak Newtonove iteracije.
	x = x0 - ajac\F0;
	if(norm(x - x0) < tol)
		break;
	end
	x0 = x;
end