function y=eros(f,x,h)
% Erotusosamääräfunktio (fd_deriv Hig-Hig s. 145)

if nargin < 3, h=sqrt(eps); end 
      % Numeerinen derivointi on häiriöaltis tehtävä, askeleena
      % ei voi käyttää pienempää kuin sqrt(eps). (eps on aivan
      % liian pieni)
 y=(f(x+h)-f(x))/h;     