% coraz lepsze(?) przybliżenia pochodnej ilorazami różnicowymi
% skrypt w MATLABie/Octave
k = 150;
x = 3; h = 1;
for i = 1:k
	h = h/2;
	blad(i) = (sin(x+h) - sin(x))/h - cos(x);
end
semilogy(abs(blad));
