function derivace
x0=pi/4;
zmena=1e-4;
for k=1:5
k
x=linspace(x0-zmena/2,x0+zmena/2,k+1);
y=sin(x);
dx=zmena/k;
dernum=diff(y,k)/dx^(k)
endfor
endfunction

