function x = newton(f, df, x0, epsilon)
% wyznacza miejsce zerowe f z tolerancja epsilon dla przyblizenia
% poczatkowego x0
	x = x0;
	while( abs(f(x)) > epsilon )
		x = x - f(x)/df(x);
		f(x) % sledzimy na biezaco postep iteracji
	end
end
