function dy= deriv_fce(f,x0,xN,n,h)
%deriv_fce - priblizny vypocet derivace funkce
%x0,xN...zacatek a konec intervalu, kde se bude derivovat
%n,,,pocet bodu pro derivaci
%h..male cislo pro aproximaci limity (implicitne h=1e-3)
%priklad:
%fce=@(x) x.^2-3*x+2    polynom
%d=deriv_fce(fce,-2,2,10,1e-5)
%x=linspace(-2,2,10)    kvuli grafu a dosazeni presne derivace
%syms t;dfce=subs(diff(fce(t)),x)   presna derivace pomoci symbolickeho
%toolboxu
%plot(x,d,'k.',x,dfce,'rx'), legend('num.derivace','presna derivace',4)

if nargin<5
    h=1e-3;
end
if nargin<4
    error('Zadejte alespon 4 vstupy!')
end
if h>1e-1 || h<1e-10
    error('Posledni vstup musi byt mezi 1e-10 a 1e-1!')    
end

x=linspace(x0,xN,n); %n hodnot v danem intervalu
dy=x;

for i=1:n
    dy(i)=(f(x(i)+h) - f(x(i)))/h;  %aproximace derivace    
end




