x=pi/4;

% x = 2;

f=@sin;             % funktionshandtag
fprimx=cos(x);      % f'(x)

% f=@sqrt;                % funktionshandtag
% fprimx=1/(2*sqrt(x));   % f'(x)

for i=1:16
    h(i)=10^(-i);
    Df(i)=derivata1(f,x,h(i));
  % Df(i)=derivata2(f,x,h(i));
    e(i)=abs( Df(i) - fprimx );
end

a=[h', Df', e']
[emin,i]=min(e)
